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Chapter 1 
Introduction 



What happens to matter when you squeeze it further and further? And what if it will 
be heated more and more? Liquid water for example will turn at some point into 
a different phase called steam when it is heated. If instead the density is increased 
at room temperature by applying an external pressure, water will subsequently turn 
into different types of ice, called ice VI, ice VII, etc. Such phase transitions are not 
specific for water, but they can take place in any interacting substance, like for example 
hadronic matter. 

Hadronic matter is matter built out of quarks and gluons. The neutron for instance 
is a form of hadronic matter, it is a bound state of two down-quarks, one up-quark 
and gluons which keep the quarks together. Imagine a hypothetical situation in which 
one has lots of these neutrons in a box. Then increase the temperature. What will 
happen? At some point the kinetic energy of the quarks which build up the neutron 
will become larger than the energy that is gained by confining the quarks inside the 
neutron. At this point the neutrons will cease to exist. The matter in the box is now in 
a new phase, which is called the quark gluon plasma. An order of magnitude estimate 
of this transition temperature can easily be made. Classically the kinetic energy of 
a quark at temperature T is about 3T/2^. Since the masses of the quarks are much 
smaller than the neutron mass, the mass of the neutron is almost completely due to 
confining energy. So the confining energy per quark is about m/v/3 « 300 MeV. This 
implies that the transition temperature Tc to the quark gluon plasma should be around 
Tc = 200 MeV, which is very hot (actually it is about 10^ times the temperature of the 
solar core). In nature these temperatures were achieved in the early universe and are 
possible in relativistic heavy ion collisions for extremely short periods of time. 

Now start again from scratch at low temperatures and squeeze the box further 
and further. At some point the neutrons will start to overlap and at even higher den- 
sities they will cease to exist as separate entities. Around this point it is expected 
that quarks will form Cooper pairs. The matter will transform into a so-called color- 



In this thesis natural units are used, so c = 1, ^ = 1 and kg = I. 
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superconducting state. An order of magnitude estimate shows that this will occur 
around densities of about mn/Vn ~ 160 MeV • fm~^ w 3 X 10^' kg • m"-'. This density 
is quite large, and takes only place in extremely dense objects like for example neutron 
stars. 

The theory which describes the interactions between the quarks mediated by glu- 
ons is called quantum chromodynamics (QCD) and is treated in more detail in the fol- 
lowing section. Using QCD one could in principle predict the behavior of matter under 
these extreme circumstances by calculating its equation of state (that is the relation be- 
tween its pressure and energy density) and its phase diagram, which are discussed in 
Sees. 1.2 and 1.3 respectively. The situations in which these extreme circumstances 
are realized in nature are reviewed in Sec. 1.4. Since it turns out that QCD is very 
complicated at the energy scales around the phase transition, this thesis will deal with 
models inspired by QCD to describe matter at high temperatures and densities, as will 
be discussed in more detail in Sees. 1.5 and 1.6. A more extensive review on QCD at 
high temperatures and densities can be found in for example Meyer-Ortmanns (1996), 
Rajagopal and Wilczek (2000) and Rischke (2004). 

1.1 Quantum chromodynamics 

Quantum chromodynamics (QCD) is a non-Abelian SU(3) gauge field theory which 
describes the interactions between the quarks. It is a generalization of Maxwells theory 
of electromagnetism. Like the electrons, quarks carry a charge, called color. Unlike 
the photons in electromagnetism, the gluons, which are the force carriers of QCD carry 
a color charge as well. As a result the gluons interact with themselves and with quarks. 
Due to the SU(3) gauge symmetry QCD has three different color charges, named red, 
blue and green. Together with the electroweak theory, QCD is one of the building 
blocks of the Standard Model of elementary particle physics. QCD is defined by the 
following Lagrangian density 

X = «A ('^ £>^ - m) + yuro) -A - \K"F''^y , 

where g is the QCD coupling constant, Ta is a hermitian generator of SU(3) and /^^ 
denotes its corresponding structure constant. The matrices mq and fi are diagonal 
and contain the current quark masses and the quark chemical potentials respectively. 
There are six different quark flavors. The up, down and strange quark are relatively 
light, while the charm, bottom and top quark are heavy. Since the masses of the heavy 
quarks are so much larger than the estimated transition temperature of about 200 MeV, 
these quarks will play a minor role at these energies and will therefore be neglected 
in this thesis. The discussion in this thesis will mainly deal with two {Nf = 2) and 
three-flavor (Nf = 3) situations. 
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The chemical potentials are necessary to describe a system at finite density. The 
baryon chemical potential yUg - Hu+Hd+Hs for example, is basically the energy it takes 
to add one additional baryon to the system. Temperature is introduced by considering a 
Euclidean space in which the xq direction is made periodic (for bosons) or antiperiodic 
(for fermions) with periodicity 1/r. Field theory at finite temperature and densities is 
discussed in more detail in Chapter 2. 

QCD has different symmetries which are reflected in the hadron spectrum as a con- 
sequence. First of all it is invariant under local SU(3) transformations. This implies 
for example that red up quarks are as heavy as blue up ones. In addition, in absence of 
quark masses and chemical potentials, QCD has a global chiral S\]{Nf)L x S\J{Nf)R 
symmetry. Moreover it has a global U(1)b symmetry related to baryon number conser- 
vation and a global U(1)a (axial) symmetry. At low temperatures and chemical poten- 
tials it turns out that the chiral symmetry is spontaneously broken down to S\J{Nf)v 
giving rise to A'^^. - 1 massless pseudoscalar Goldstone modes. For Nf = 2 these 
are the three pions, for Nf = 3 also the four kaons and the 77 particle are among the 
pseudoscalar Goldstone modes. If chiral symmetry is broken the (ilu), (dd), and (ss) 
condensates obtain a vacuum expectation value. 

However, in reality the quarks have a small mass. Therefore chiral symmetry is 
only an approximate symmetry, as a result the pions, the kaons and the 77 particle be- 
come massive. This remaining (approximate) S\J{Nf)v symmetry is the reason why 
the constituent quark model of Gell-Mann works as well as it does. Particles are eigen- 
states of the QCD Hamiltonian which due to the symmetry commutes with S\J{Nf)v- 
Hence the particles can be classified by the representations of S\J{Nf)v- At high tem- 
peratures and/or densities the chiral symmetry is approximately restored, giving rise 
to a phase transition. Although the U(1)a symmetry is broken due to nonzero quark 
masses as well, it also has another reason of breakdown. The non-trivial topological 
vacuum structure of QCD due to instantons is causing U(l) axial symmetry breaking 
too, which explains the relatively high mass of the 77 meson ('t Hooft, 1976). 

One of the mysteries in QCD is confinement. It turns out experimentally that 
hadrons, which are bound states of quarks, carry no color. Colored objects, like freely 
moving quarks, do not occur in nature at low energies. In numerical computations 
(lattice QCD) it is confirmed that QCD has this confinement property. But a detailed 
understanding of the confinement mechanism is still lacking. At high temperatures 
and/or chemical potentials it is expected that matter will be in a deconfined phase, 
which means that in that situation quarks are liberated from the hadrons. Whether the 
deconfinement phase transition for light quarks coincides with the chiral symmetry 
restoration transition is an important issue which has not yet been resolved. 

QCD is asymptotically free, this implies that the effective coupling of quarks to 
gluons becomes smaller at high energies. So at high energy scales, larger than about 
I GeV, QCD is a theory of weakly interacting quarks and gluons. Due to the small 
coupling constant it is possible to perform calculations in this regime using perturba- 
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tion theory. But, at lower energies QCD becomes strongly coupled and perturbation 
theory breaks down. The order of magnitude estimate of the transition temperature 
Tc ~ 200 MeV from the first paragraph tell us that QCD is likely to be strongly cou- 
pled around the phase transition. Hence it is expected that perturbation theory will fail 
to describe QCD near T^. This will be illustrated next by comparing perturbative and 
lattice calculations of the QCD pressure. 

1.2 QCD equation of state 

Two important macroscopic thermodynamical quantities are the pressure P and the 
energy density £ of QCD. The relation between & and P is called the equation of 
state and determines for example the behavior of matter created in a relativistic heavy 
ion collision, the properties of a (neutron) star and the evolution of the early stage 
of the universe shortly after the big bang, see Sec. 1.4. Due to asymptotic freedom, 
QCD describes a gas of weakly interacting quarks and gluons in the limit of very high 
temperature. In that case the coupling constant is small, and hence perturbation theory 
should be applicable. However as will be shown next perturbation theory does not 
work for the temperatures of in the neighborhood of T^. 

Perturbative calculations 

The results of the perturbative calculation of the pure glue (which means no quarks) 
QCD pressure compared to numerical calculations (lattice QCD) are displayed in 
Fig. 1.1 up to order a^^^ = (g-^/An)^^^. It can be seen in the figure that for temper- 
atures in the neighborhood of the transition temperature T^ the results are varying a 
lot upon inclusion of higher order corrections. This implies that the perturbative ex- 
pansion breaks down, even for very high temperatures of around lOOOr^. However, 
by reorganizing the perturbative series using resummation methods (Andersen et al, 
1999, 2002) and self-consistent approaches (Blaizot et al, 1999) based on hard thermal 
loops (Braaten and Pisarski, 1990) it is possible to obtain reliable results for T > 3Tc. 
In perturbation theory it is impossible to calculate the order g^ contribution because an 
infinite number of diagrams are contributing in this order and cannot be resummed as 
was argued by Linde (1980). 

It is possible to perform perturbation theory for finite chemical potentials as well, 
but it again fails for densities where the phase transition occurs. Also this perturbation 
series can be improved by applying the so-called hard dense loop resummation method 
(Andersen and Strickland, 2002). 

Lattice calculations 

The best known method to obtain the QCD thermodynamical quantities from first 
principles near the phase transition temperature is by lattice calculations. In these 
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Figure 1.1: Perturbative results (indicated by their orders) and lattice results of (Boyd 
et al, 1996) (indicated with black diamonds) of the pure glue QCD pressure normal- 
ized to the ideal gas value, as a function of T/Tc. The shaded regions arise due to 
varying the MS renormalization scale between nT and 47tT and are hence an indi- 
cation of the error that is inherent to truncating the perturbative series. The order a 
calculation was performed by Shuryak (1978), a^^^ by Kapusta (1979), a^^^ logo- by 
Toimela (1983), a^ by Arnold and Zhai (1995), a^^^ by Zhai and Kastening (1995) and 
Braaten and Nieto (1996) and order a^ log a (not shown in figure) by Kajantie et al. 
(2003). This figure is adapted from Andersen et al. (2002). 



calculations, space-time is discretized and replaced by a lattice of a finite size. The 
fermion fields live on the vertices of this lattice, the gauge fields are replaced by links 
connecting the different vertices. All thermodynamical quantities can be obtained by 
numerically calculating the partition function (discussed in Chapter 2). This requires 
integration over the fermion fields and links. Since this results in a huge number of 
integrations, typical lattice sizes are taken to be rather small, in the order of ten points 
for each dimension. Due to this modest lattice sizes, the particles with low mass (like 
the pion), which can propagate over longer distances are not very well described. The 
integrations in lattice calculations are performed statistically, using importance sam- 
pling Monte-Carlo methods. This works fine for QCD at zero chemical potential. 
But at finite baryon chemical potential, the contribution from the fermions, the so- 
called fermionic determinant, becomes imaginary (see also Chapter 7). As a result, 
the integrand of the partition function becomes oscillatory, which hampers the impor- 
tance sampling methods. This complication is called the fermion sign problem. Small 
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Figure 1.2: Lattice calculation by Karsch et al. (2000) of the pure glue and full QCD 
pressure (for 2 light quark flavors, 2 light & 1 heavy flavor and 3 light flavors), nor- 
malized to T^, as a function of T . The arrows indicate the limit of the pressures when 
T ^ oo and are calculated in Sec. 2.4. This figure is adapted from Karsch (2002). 



baryon chemical potentials are, however, accessible by making a Taylor expansion 
around /ub = (AUton et al, 2002; Fodor and Katz, 2002; de Forcrand and Philipsen, 
2002; D'Eha and Lombardo, 2003). 

In Fig. 1.2 the lattice results of the pressure of QCD for diff"erent numbers of quark 
flavors are displayed as a function of T for /i = 0. Here T^- ~ 170 MeV. From the 
figure it can be seen that the QCD pressure rises quickly after passing the transition 
point. This may be an indication that many new degrees of freedom are formed, just 
what one would expect if the quarks become deconfined from the hadrons. The figure 
shows that even around 4Tc the pressure is still far away from that of a freely inter- 
acting gas of quarks and gluons P = !Psb- While the lattice gives reliable results for 
T > Tc, the data at T < Tc can not be trusted. This is because the pion, the lightest par- 
ticle of QCD, which is expected to dominate the pressure of QCD at low temperatures 
is still far too heavy on the lattice. 

For T > Tc the lattice data for the QCD pressure can be fitted to quasiparticle 
models (Peshier et al, 1996; Levai and Heinz, 1998; Schneider and Weise, 2001). 
The results of these fits can be extended to finite chemical potential (Rebhan and Ro- 
matschke, 2003; Thaler et al, 2004). In this way the equation of state of QCD at finite 
chemical potential for T > Tc can be predicted. 

It is generally believed that the lattice calculations of the pressure are reliable for 
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temperatures around T^ and higher. The results at high temperatures can be predicted 
by hard-thermal-loop resummation methods. At low temperatures lattice QCD is still 
unreliable due to the unrealistically large pion mass. In order to predict the pressure of 
QCD at low temperatures one can use a low-energy effective theory in order to describe 
the hadron gas phase as is done in this thesis in Chapter 6 and will be discussed in more 
detail in Sec 1.5. 



1.3 QCD phase diagram 

By considering the behavior of order parameters (parameters which vanish in one 
phase and are non-vanishing in another) one can determine a phase diagram. One 
typically distinguishes between two different types of phase transitions, a first order 
phase transition in which the order parameter changes discontinuously and a second 
order phase transition in which the derivative of the order parameter changes discon- 
tinuously. Higher order phase transitions are also possible, but are often called second 
order transitions as well. In addition, cross-over transitions in which the order param- 
eter changes smoothly can occur. The different possibilities are displayed in Fig. 7.1. 

In the limit of zero quark masses the chiral condensates, (mm), (dd) and (ss) are 
order parameters for the breaking of chiral symmetry. In Sec. 1. 1 it was mentioned that 
the chiral symmetry is already explicitly broken in the QCD Lagrangian density due 
to the non-zero quark masses. In that case the chiral condensates are only approximate 
order parameters. The order parameter for the confinement/deconfinement transition 
in the limit of infinitely heavy quarks is the trace of the so-called Polyakov loop. For 
finite quark masses no order parameter for this transition is known, see for example 
Weiss (1993). 

In Fig. 1.3 the current understanding of the QCD phase diagram is displayed 
schematically. As was discussed in the previous section, only results for zero and 
small baryon chemical potential can be obtained from lattice QCD. Lattice calcula- 
tions find a cross-over transition at T^ = 170 MeV. The rest of the phase diagram is 
not yet obtained from first principles QCD but can be estimated by means of effective 
models like the NJL model studied in Chapter 7. However, the phases with tempera- 
tures and densities much higher than the densities and temperatures where the phase 
transition takes place are accessible by hard-thermal loop and hard-dense loop resum- 
mation techniques as was discussed in the previous section. The NJL phase diagram 
as a function of jUg and T is displayed in Fig. 7.4. 

As mentioned before, finite chemical potentials are needed to describe a system 
at finite density. However, it is important to keep in mind that the relation between 
chemical potential and number density is not linear. Especially at a first order phase 
transition, a single value of a chemical potential can correspond to a whole range of 
densities, as is illustrated for the NJL model in Fig 7.2. In this case one also speaks of 
a mixed phase, two phases can occur together. The world we live in is an example of 
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Figure 1.3: Schematic structure of the current understanding of the phase diagram of 
QCD. First order phase transitions are indicated with a solid line, second order with a 
dashed line and a cross-over with a dotted line. In the diagram the cooling trajectory 
of the early universe and of matter produced in heavy ion collisions is sketched. It is 
also roughly indicated in which phase the matter inside neutron stars can be. 



a mixed phase of nuclear matter and vacuum. One should be aware that the real phase 
diagram of matter is not just the QCD phase diagram. To obtain the complete phase 
diagram of matter, one should also take into account the electromagnetic and weak 
interactions. 

The best known point in the QCD phase diagram, is the transition from the vacuum 
to the nuclear matter phase, there is a first order phase transition at fig - 922 MeV, 
see also Halasz et al. (1998). 

It is illustrated in Fig. 1.3 that matter at low chemical potentials and temperatures 
matter is in a confined phase in which chiral symmetry is broken as well. If the temper- 
ature is increased, matter goes according to the current understanding via a cross-over 
transition to the deconfined and chirally symmetric phase at low chemical potentials, 
and via a first-order transition at higher chemical potentials. This deconfined phase is 
called the quark-gluon plasma. The point in which the first-order transition goes over 
to a cross-over is called the critical endpoint. It still is uncertain where this critical 
endpoint lies exactly in the phase diagram. At low temperatures and high chemical 
potentials due to an attractive interaction, quarks can form Cooper pairs just like elec- 
trons in ordinary superconductivity. This phenomenon will be discussed in more detail 
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in Chapter 7. 

The phase diagram in Fig. 1.3 is displayed as a function of baryon chemical po- 
tential and temperature. Of course it is interesting to investigate other phase diagrams 
as well, for example as a function of the quark masses, as is discussed in Laermann 
and Philipsen (2003). In Chapter 7 phase diagrams of the NJL model will be investi- 
gated for unequal chemical potentials and temperature. In that case a new possibility 
appears which is not present in Fig. 1.3, namely quarks can form pseudoscalar con- 
densates, like the pion condensate (uiysd). 

In Fig. 1.3 it is also indicated which part of the phase diagram can be investigated 
using relativistic heavy ion collisions, which kind of matter neutrons stars are presum- 
ably made of, and through which phases the early universe went shortly after the big 
bang. In the following section these three situations will be examined in somewhat 
more detail. 

1.4 Matter under extreme conditions 

Typical situations in which temperatures and densities could be high enough for de- 
confinement to occur, are the universe just after the big bang, during heavy ion col- 
lisions and inside very compact neutron stars. In this section a short overview of 
these situations is given. More extensive discussions can be found in for example El- 
lis (2005) (the big bang in relation to heavy ion collisions), Gyulassy and McLerran 
(2005) (heavy ion collisions) and Weber (2005) (neutron stars). 

The big bang 

In one of the earliest stages of the universe, about 10"^ seconds after the big bang, the 
universe was still so hot that the matter inside was in the deconfined phase, i.e. the 
quark gluon plasma. Since at that time the particles and antiparticles had not annihi- 
lated yet, the baryon chemical potential was very small. As is indicated in Fig. 1.3 
when the universe cooled, it probably went through a cross-over transition to the con- 
fined phase. Since a cross-over transition is smooth, it is unlikely that the expansion 
of the universe was modified substantially during this transition. 

To describe the evolution of the universe one has to use an equation of state. Most 
often a simple equation of state is used, P = £/3 for the radiation dominated era 
at early times, and !P = for the matter dominated era which occured later. The 
description can be made more realistic by using the QCD equation of state. 

Relativistic iieavy ion collisions 

The behavior of matter under extreme circumstances can be studied experimentally us- 
ing heavy ion collisions. Such experiments have been performed at the Super-Proton- 
Synchrotron (SPS) at CERN and are being performed at the Relativistic Heavy Ion 
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Collider (RHIC) at BNL. Currently new accelerators which, among other things, will 
be used for relativistic heavy ion collision studies are being build at CERN (the Large 
Hadron CoUider (LHC)), and at GSI (Schwerionen-Synchrotron (SIS 200)). 

In a typical relativistic heavy ion collision two incoming heavy nuclei (for example 
gold with 197 nucleons) collide at relativistic energies. At RHIC these energies are up 
to 200 GeV per nucleon. During the collision a large fraction of the kinetic energy 
is converted into particles. Therefore statistical methods can be used to describe the 
system. High temperatures and energy densities are achieved during these collisions. 
The baryon chemical potential remains low in a heavy ion collision. The reason for this 
is that due to the large production of particle antiparticle pairs, the initial dominance of 
particles over antiparticles is washed out. Since SIS 200 will operate at a lower energy 
than RHIC and LHC, a higher baryon chemical potential can be achieved. However 
as a result the final temperature will be lower, which could make it more difficult to 
probe the phase transition. 

At RHIC it seems that the produced matter quickly achieves thermal equilibrium. 
After that moment, relativistic hydrodynamics can describe the evolution using the 
QCD equation of state. This indicates that the matter created during the collisions at 
RHIC has a very low viscosity, possibly the most perfect fluid ever made. 

That a statistical description to model heavy ion collisions works, is illustrated 
in Fig. 1 .4. In that figure, a prediction of the ratio of particle yields is compared to 
the experimental data of the four different experiments at RHIC (Braun-Munzinger 
et ai, 2001). The statistical model only has two fit parameters, T k 174 MeV and 
fiB ~ 46 MeV. It can also be seen from the figure that the baryon chemical potential 
is not that big compared to the temperature, because the ratio between particles and 
antiparticles of the same hadronic species are close to unity. 

The main objective of these heavy ion collision experiments is to produce the quark 
gluon plasma and measure its properties. One of the clearest signals for the creation 
of a very dense and hot state of matter comes from the suppression of back-to-back 
correlations between two high transverse-momentum jets. In a collision at low energy 
densities, such a jet should be correlated with a jet produced in the opposite direction 
due to momentum conservation. This is indeed observed in heavy ion collision experi- 
ments. However, at energies of 200 GeV per nucleon in central gold-gold collisions at 
RHIC, this correlation has suddenly disappeared. This indicates that the momentum of 
the opposite jet is absorbed in a very hot and dense medium, probably the quark-gluon 
plasma. 

Using heavy ion collisions it is difficult to reach high chemical potentials, so in 
order to investigate that situation one has to look to extremely compact objects like 
neutron stars. 
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Figure 1.4: Prediction of ratios of produced particles for gold-gold collisions at V^ = 
130 GeV by a statistical model compared to results of the four different experiments 
at RHIC. This figure is adapted from Braun-Munzinger et al. (2001). 



Neutron stars 



If the mass of a star is larger than about ten times the solar mass, the fusion process can 
continue until an iron-nickel core is formed. Then fusion will stop since iron is the nu- 
cleus which has the lowest binding energy per nucleon. As a result the temperature of 
the core will drop, hence the pressure will go down. Then the gravitational interactions 
cause the core to collapse until nuclear densities are reached. At this point the collapse 
stops because it takes a lot of energy to squeeze the core further. This creates a shock 
wave which as a result emits all of the matter from the shells surrounding the collapsed 
core. This event is called a supernova explosion. The remaining core cools down and 
becomes a neutron star or at even higher densities a black hole. Typical densities could 
be so high that the core of a neutron star is a color superconductor. Using the QCD 
equation of state at high-densities, the so-called Tolman-Oppenheimer-Volkov relation 
can be applied to calculate the mass as a function of the radius of a neutron star, see for 
example Fraga et al. (2001) and Andersen and Strickland (2002). These mass-radius 
relationships can be compared to observations. However, until now there has not yet 
been discovered a neutron star from which one is certain that the inner core is so dense 
that it must be in a color superconducting phase. 
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1.5 QCD inspired theories 

As was argued in the previous sections, due to the nonperturbative nature of QCD near 
the phase transition, it is not known how to obtain for example the order of the phase 
transition, the equation of state and the phase diagram using analytical methods from 
first principles QCD for all temperatures and chemical potentials. Also the confine- 
ment/deconfinement mechanism is not understood analytically. However, toy models 
(models which share features with QCD) and low-energy effective theories (theories 
which describe QCD in the non-perturbative regime) can be studied to learn about 
certain aspects to be addressed below. 

Toy models 

Toy models for QCD are models which have features in common with QCD. One can 
study these models in order to learn about nonperturbative methods and phenomena 
like for example the behavior of thermodynamical quantities, the generation of a mass 
gap, confinement and the importance of topological configurations like instantons. In 
this thesis two of these toy models are studied, the 0(A'^) nonlinear sigma model in 
1 + 1 dimensions in Chapter 3 and the CP'^^^ model in 1 + 1 dimensions in Chapter 4. 
The nonlinear sigma model is a real scalar field theory which like QCD is asymptot- 
ically free and has a dynamically generated mass gap. The CP^~^ model contains 
complex scalar fields and U(l) gauge fields. This model is also asymptotically free, 
has a dynamically generated mass gap, contains instanton configurations and has the 
confinement property. In these two models it is possible to expand in the number of 
fields N. This expansion is called the \aige-N approximation and is a method which 
can give insight in the nonperturbative behavior of these theories. This method will be 
used throughout this thesis and is discussed in more detail in Chapter 3. In this thesis 
the pressure for both models is calculated to next-to-leading order in 1 /N. In this way 
it is possible to check the validity of the 1 /N expansion and to investigate the effects 
of instantons on the thermodynamical quantities in the CP^~^ model. 

Low-energy effective theories 

Low-energy effective theories describe QCD in the non-perturbative regime where 
perturbation theory is no longer applicable. The simplest examples of such an effective 
theory are the 0(4) linear and nonlinear sigma model in 3 -H 1 dimensions. The 0(4) 
symmetry (which is locally isomorphic to the chiral SU(2)i x SU(2)« symmetry of the 
two-flavor QCD Lagrangian) of these models is spontaneously broken to 0(3) (which 
is isomorphic to the remaining chiral SU(2)y symmetry of the QCD vacuum state). 
Since these models have the same symmetry breaking pattern as two-flavor QCD, they 
serve as a low-energy effective theory for 2-flavor QCD in which only the pion and 
the sigma meson can occur. Moreover for the same reason, the phase transition of the 
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0(4) (non)linear sigma model falls in the same universality class as the 2-flavor QCD 
chiral phase transition. This allows one to study the order of the phase transition and 
critical exponents of QCD using the 0(4) (non)linear sigma model. In this thesis the 
pressure of these models is calculated to next-to-leading order in 1 IN in Chapter 6. In 
this way a prediction for the pressure of QCD at temperatures below Tc is made, where 
the lattice calculations are not reliable. 

Another effective theory studied in Chapter 7 of this thesis, is the Nambu-Jona- 
Lasinio (NJL) model. In this model the gluon exchange between the quarks is replaced 
by a 4-point quark interaction. As a result some features of QCD like confinement and 
asymptotic freedom are lost. However, low-energy properties like the meson masses 
are described very well using this model. Furthermore this model has the same pattern 
of chiral symmetry breaking as QCD. Therefore, it is expected that the NJL model 
gives a realistic qualitative description of the QCD phase diagram for low temperatures 
and chemical potentials. In Chapter 7 phase diagrams with different up, down and 
strange quark chemical potentials are calculated, in order to study the competition 
between phases in which color superconductivity is possible and in phases where the 
pions or kaons condense. 
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To summarize, the thermodynamics of QCD inspired theories is studied in this the- 
sis. In Chapter 2 a short introduction to finite density and temperature field theory is 
given. Particular emphasis is put on the analytic and numerical calculation of a com- 
bination of a sum and an integral, which will be required frequently throughout this 
thesis. In Chapter 3 the effective potential is derived from which by minimization one 
can derive the thermodynamical quantities like the pressure and determine the phase 
diagram. Almost all calculations in this thesis are based on evaluating this effective po- 
tential. Moreover, Chapter 3 discusses the 1 IN approximation which can give insight 
into non-perturbative physics. This is also a basic ingredient for all the calculations 
performed in this thesis. The thermodynamics of the 0(A'^) nonlinear sigma model in 
1-1-1 dimensions is studied in Chapter 4. In Chapter 5, the effect of quantum ins tan- 
tons on the thermodynamical quantities is investigated using the CP^~^ model in 1 -i- 1 
dimensions. Chapter 6 is devoted to the study of the thermodynamics of the O(A^) 
linear and nonlinear sigma model in 3 -i- 1 dimensions, which can be used to predict 
the pressure of QCD at low temperatures. The NJL model and its phase diagrams are 
discussed in Chapter 7. 
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1.7 Notations 

Several notations are used throughout this thesis. These will be summarized here. 

• Euclidean momentum vectors are denoted by a capital letter, that is P = (po, p). 
The length of a momentum vector /? is denoted us p -\p\. 

• The integral over Euclidean momenta, the bosonic sum-integral and the fermio- 
nic sum-integral are respectively defined as 

(1.2) 

• The difference of a sum-integral and an integral are for bosonic fields and fermi- 
onic fields respectively defined as 

i f{p)^i- r, i ^i - { ■ (1-3) 

Jp Jp Jp J{P] J{P] Jp 

• Integration over space(-time) will be depending on the context written as 

r= fdxQ fd^x or f=f dxo fd^x. (1.4) 



Chapter 2 

Finite temperature and density field theory 



In a system with a large number of particles like for example a gas, it is very cum- 
bersome, if not impossible, to calculate the trajectories of individual particles. On the 
other hand collective properties, like the pressure and number densities, characterize 
the system of particles as a whole and are therefore in many cases much more inter- 
esting than the behavior of individual particles. These collective properties can really 
be calculated using statistical methods. 

In this chapter first the basics of classical statistical physics will be summarized. 
Then, by using path-integrals in Euclidean space-time, classical statistical physics will 
be cast in a form suitable for quantum field theories, called finite temperature and den- 
sity field theory. As an illustration the pressure of a free scalar and of a free fermion 
field theory are obtained. Both expressions for the pressure are explicitly evaluated 
after a general explanation of how frequency sums, that naturally arise in finite tem- 
perature calculations, can be computed. At the end of this chapter further techniques 
for the evaluation of frequency sums are developed, which are useful for the numerical 
computation of more complicated sums that arise in interacting field theories. 

A more extensive introduction to finite temperature and density field theory can be 
found in the books by Kapusta (1989) and Le Bellac (2000) and the review article of 
Landsman and Van Weert (1987). 

2.1 Classical statistical physics 

Consider a box of volume V, having total energy E and filled with N particles. This 
box can be divided into regions 1 and 2. Clearly, it holds that V = V\ + V2,E - E1+E2 
and N - N\ + A^2- Let Q.{E, V,N) be the number of ways in which the total energy E 
can be distributed over N particles in a volume V. This quantity Q. is called the number 
of micro-states. Now the two postulates of statistical physics are, 

1 . All micro-states are equally likely to occur. 

2. In equilibrium the system will choose the state that is the most likely to occur. 
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Combining postulate 1 and 2 gives that the equilibrium state is the one with the highest 
number of micro-states. The number of micro-states of the complete box can be writ- 
ten as the product of the micro-states of the two different regions, Q - QiQ2» where 
O,- = 0.i{Ei, Vi,Ni). It is convenient to turn this equality into an additive relation by 
introducing a quantity called the entropy S which is defined as S = log Q. Then the 
total entropy of the box S is equal to the sum of the entropies of the different regions, 
S = S\ + 82- The two postulates can be translated into the condition that the total 
entropy is maximal in equilibrium. If the entropy is maximal it holds that 

{dEijy^N \5^i/y,,iv, \^^1/V2,/V2 \^^i/y,,M \^^2/v,,;v, ' 

where it was used that the total energy E = Ei + E2 is constant. It follows that in 
equilibrium 

ffi] ^(^] . (2.2) 

In equilibrium the temperatures T, of the two regions 1 and 2 should be equal, that 
is Ti = T2- Hence dS\ldE\ should be some function of temperature. The correct 
definition of temperature turns out to be 

Ti \dEjy^j^^ 

because in that way it is possible to derive the experimentally verified ideal gas law, 
PV = NT. Similarly since N = N\ + N2 is constant it holds that in equilibrium 

(2^) =(^) . (2.4) 

In equilibrium it take as much energy to transfer one particle from region 1 to region 
2 as to do the opposite. This energy is called chemical potential, so in equilibrium the 
chemical potentials should be equal, that is jui = //2- As a result dS i/dNi should be 
some function of chemical potential. It turns out that the correct definition is 

^^-Mw) ■ ^2.5) 

because it gives rise to the correct distribution function for fermions, Eq. (2.22). 

Now consider a box of fixed volume V placed in a very large heat bath of constant 
temperature T and constant chemical potential jj. The box is allowed to exchange 
energy and particles with the heat bath. The total system of heat bath and box together 
has energy Eq and contains No particles. The probability p,- that the box has energy 
Er and contains A^^ particles is equal to the probability that the heat bath has energy 
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£"0 - Er and contains A^o - A^r particles. So it follows that pr is proportional to 0.{Eq - 
Er, No - Nr), the number of micro-states of the heat bath. In terms of entropy one has 
that 

Pr = C-^ exp [S (Eo -Er,No- Nr)] , (2.6) 

where C is a normalization factor. Assuming the heat bath is large implies that Eq » 
Er and A^o ^ Nr. Hence it is possible to expand the entropy of the heat bath around 
Eq and A'^o. 



S{Eo-Er,No-Nr) = S(Eo,No)-E, 



dS{E,N) 



dE 



EoM 



-Nr 



dS{E,N) 



dN 



Eq,No 



(2.7) 



In higher orders of the expansion one gets terms like d^S /dE'^lEg^No * d{l/T)/dE, 
which reflects the change of the heat bath temperature when energy is transferred into 
the box. Because it is assumed that the heat bath has constant temperature d{l/T)/dE = 
0, so this higher order term can be neglected. Another higher order term of the expan- 
sion is d'^S /SN'^Ieq^No ~ -d(j^/T)/dN which reflects the change of chemical poten- 
tial divided by temperature when particles are transferred into the box. Again be- 
cause it is assumed the heat bath has constant temperature and chemical potential, this 
higher order term can be neglected too. For the same reason one can assume that 
d^S /dEdN\Eo,NQ vanishes as well. 

As a result, the probability pr that the box has energy Er and contains A^^ particles 
is equal to 



Pr^Z~'exp[-j3{Er-fiNr)] , 



(2.8) 



where Z is a normalization factor different from C andyS = l/T. 

The normalization factor Z which is also called the partition function, is equal to 
the sum of all probabilities. 



Z = Y,^xp[-^(Er-iuNr)] 



(2.9) 



The partition function contains all information of the collective or macroscopic behav- 
ior of the thermodynamic system. Strictly speaking this is already a quantum mechan- 
ical equation, since the energies are assumed to be discrete. In the classical case, one 
has to replace the sum over states by an integral. Using the partition function one can 
calculate thermodynamical quantities, like the energy density of the box. 



^=vT,p'-^'- = 



]_d\ogZ 
V d/3 



(2.10) 



the number density of particles in the box, 

1 ^ 1 51ogZ 



/3V dp 



(2.11) 
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and the entropy density, 

S=yl.PAogpr^--^^. (2.12) 

Using the definition of the pressure, AEr = -PrdV which is valid at constant tempera- 
ture and number density, it follows that the pressure is given by 



„ v'' oEr 1 d log Z 



Typically, the width L of a system is much larger than the inverse temperature, (i.e. 
L » 27r/3), such that one can use the infinite volume limit to describe the thermo- 
dynamics of a finite volume to good approximation. The advantage of the infinite 
volume limit is that field theoretic calculations simplify. In all calculations performed 
in this thesis, this infinite volume limit is taken. Then it turns out that logZ becomes 
proportional to V, such that the pressure becomes 

P^ — logZ. (2.14) 

Instead of calculating log Z//iV directly, in this thesis the pressure will be calculated 
via the effective potential (see Chapter 3) which in its minimum equals log Z//iV. 

2.2 Quantum statistical physics 

Of many physical systems one does not know the energies Er and the number densities 
Nr exactly. Most often only the Hamiltonian H and a corresponding number operator 
Q which commutes with H is known. Denoting the eigenstates of H and 2 by |r) , the 
partition function expressed in terms of the Hamiltonian and number operator becomes 

Z = 2](r|e-^(^-^*)|r) = tr e'^i^-'^^^ . (2.15) 

r 

The thermal expectation value of an operator A can also be expressed in terms of a 
trace, 

(A)-itr[Ae-^(^-^^)] , (2.16) 

An expectation value is independent of the choice of basis due to the cyclic property 
of the trace. 

As an example of the formalism, it will be shown how to derive the Bose-Einstein 
distribution function. This distribution function gives the number of states as a function 
of energy and temperature for a non-interacting bosonic system. The Bose-Einstein 
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distribution function is the expectation value of the number operator. Non-interacting 
bosons obey the following harmonic oscillator Hamiltonian 

H = ^^ojk [aj^al + alcij^) , (2.17) 

k 

where cok is the energy of the state with momentum k and a^ and dj^ are respectively 
annihilation and creation operators which satisfy the usual commutation relation for 
bosonic operators, [a^, aj] - 6ki and [a^, d^] = [aj, aj] = 0. The bosonic number 
operator is given by Nj^ = S^a^. The expectation value of the number operator is 

(N,) = i tr (e~^(^~^^^hld^) = - tr {a^e-/'(^'^^hlA"-^^^e'''(^-^'^^) , (2. 1 8) 
With use of the following equality 

e'^Ae''^ =A + [B, A] + \[B, [B, A]] + ... , (2.19) 
and the fact that [H - ^N, aj^] = {a>k - yu) a^ it can be shown that 

<A^,) = -tr (a,a^e-^("*-^V-^(^-^^)) - <1 - A^,)^-^^'^*-'^) . (2.20) 

The Bose-Einstein distribution function n(a)k) = (Nt) follows from the last equation, 

n{ci)k) - -jT ^ ■ (2.21) 

In a similar way it is possible to derive the Fermi-Dirac distribution, which is the 
expectation value of the fermionic number density. The fermion creation and annihila- 
tion operators satisfy anti-commutation relations. As a result one picks up a minus sign 
in Eq. (2.20) when swapping the annihilation and creation operators. The Fermi-Dirac 
distribution is given by 

n{a)k) = —;t, ; ■ (2.22) 

2.3 Statistical field theory 

Using the path integral formalism it is possible to obtain the partition function of a field 
theory. Consider a bosonic field ^{t, x). The thermal expectation value of a product of 
two bosonic fields in equilibrium with a heat bath of temperature T is given by 

{4>{tuxi)kt2, X2)) = -tr [0(?2, ^2)e-'^^khJi)e'^^e-P^] . (2.23) 
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The dynamics of a field ^{t, x) is entirely described by its Hamiltonian H, which can 
be used to determine the time evolution of the fields, 



^{t, x) - e''^^(0, x)e-''" 



(2.24) 



By identifying t = ijS a connection between inverse temperature and imaginary time is 
found. As a result 



{4>(tuXi)4>{t2,X2)) ^ {hh,X2)4>(ti + iP,xi)) . 



(2.25) 



Using the relation between inverse temperature and imaginary time the partition func- 
tion can be written in terms of a path integral. For this, consider a transition matrix 
element between an initial bosonic state cpi and a final state cpf in ordinary field the- 
ory. Such a transition element in terms of a path integral is given by the following 
expression 



{0/|exp[-/(f/-f,)^]|0,) - I D> exp i ( dt ( d'^x £[(/>] 



(2.26) 



where the prime (') on the measure indicates that the path integral is taken over fields 
which satisfy the following boundary condition (p{ti, x) - {4>i\4>{ti, x)\4>i) and (p{tf, x) - 
{(pf\^{tf, x)\(pf). If one makes the identification t = -ir and if one chooses tj = and 
tf = -ifi one finds 



(0/|exp[-y6i/]|0,) - I £)> exp - I dr I d'^x£[(p] 



(2.27) 



The last equation enables one to write the partition function in terms of a path integral. 



Z = 2] (0„ |exp[-/3i/]| 0„) = Td^ exp - f dr fd^xXL^] 



(2.28) 



where the integration is implicitly over all fields which obey the condition 0(t = 
0,x) - ±(p{T - /3,x). Since |0) and -|0> describe the same physical state the sign 
of the boundary conditions on the bosonic fields cannot be determined in this way. 
However, this can be done by considering a two-point function. 

For bosonic fields the two-point function evaluated at t = and t, where t is 
between and /3, is given by 



{Trko)kr)} = (mko)} = {mm) = (Trmm) , 



(2.29) 



where Eq. (2.25) was used. Here Tt indicates time ordering in imaginary time. By 
choosing t = y6 it follows from Eq. (2.29) that the boundary condition on bosonic 
fields has a -i- sign, 

0(r - 0, x) = 4>{t =I5,x) . (2.30) 
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For fermions similar arguments can be used, but since time ordering for fermions 
requires an additional minus sign, one finds an anti-periodicity condition for fermionic 
fields 

)/,(t - 0, x) - -ilf(T =/3,x) . (2.31) 

So thermal field theory is in essence a Euclidean field theory where one dimension 
(t) is compactified to a circle. As a consequence of this, the Fourier transform of a 
bosonic field becomes a sum over so-called Matsubara frequencies, 

^^^' ^^ = ^T,j l^^"'""''''''^^^) ^ £ e"^"^^'^-^^{K) , (2.32) 

where the Matsubara frequencies are 6J„ = InnT. The capital ^ is a momentum 
vector in Euclidean space, K - {oJn,k)- The symbol £^ denotes a so-called sum- 
integral, where the sum is over bosonic modes, and will arise often in finite temperature 
calculations. The momentum representation for fermions is given by 

^{t, x) = ]-Y{ ^e'^-^'^-'^iK) ^ { e'^'''^'^MK) , (2.33) 

where the Matsubara frequencies for fermions are a)n - (2« -i- l)7rr. The symbol £. 
denotes a sum-integral, where the sum is over fermionic Matsubara modes. 

Two different integration contours are often used in equilibrium finite temperature 
field theory. In the derivation above, the Matsubara contour was used. This is a contour 
starting at f = straight down the imaginary axis to t = -ip, which gives rise to the 
so-called imaginary-time formulation of thermal field theory. Another possibility is 
the Keldysh contour which starts at ti = -oo, goes along the real axis to ?i = oo, down 
to t2 = ti - ie, back under the real axis to ti, = tj - ie and finally to tf - ti - j/6. 
This Keldysh contour gives rise to the so-called real-time formalism. The real-time 
formalism is favored over the imaginary-time formalism when quantities have to be 
obtained in Minkowskian space-time at finite temperature as is for example the case for 
spectral densities. To calculate such a spectral density in the imaginary-time formalism 
one has to make an analytic continuation, which can be avoided by using the real-time 
formalism. All calculations in this thesis will be performed using the imaginary-time 
formalism. 

As an example of finite temperature field theory, the pressure of a free scalar theory 
will be calculated. The Lagrangian density of this theory in EucUdean space is given 
by 

£ = ^d^(t>d^(P + ^m202 , (2.34) 

where m is the mass of the scalar field. The action of a free field theory is quadratic in 
the fields, hence the Gaussian path integral can be computed exactly (see for example 
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Weinberg (1995), Chapter 9). One finds using the fact that logDet A = TrlogA 

logZ = Q = -ixr \og[-d' + n?) , (2.35) 

where the single closed loop denotes the corresponding Feynman diagram of this con- 
tribution to log Z. In contrast to for example a two-point function, log Z contains no 
external vertices, so all its diagrams are necessarily closed. In an interacting field 
theory more complicated closed loop diagrams contribute to the pressure next to the 
single closed loop, for examples see Figs. 3. 1 and 3.2. The Feynman rules needed to 
evaluate these kind of loop diagrams at finite temperature can be found in for example 
Kapusta (1989). The functional trace in Eq. (2.35) is over a complete set of functions 
that satisfy the periodic boundary conditions in imaginary time for scalar fields. The 
trace can be evaluated by going to momentum space. As a result 

log Z = -^£ log(p2 + m^) , (2.36) 

were the sum-integral is defined in Eq. (2.32). The pressure can now be calculated 
by applying Eq. (2.14). Since it is only possible to measure pressure differences, it is 
convenient to normalize the pressure at zero temperature to zero. In order to do this 
the contribution at zero temperature which is 



lip 



log(p2 + m^) , (2.37) 



will be subtracted from the contribution at finite temperature which is 



\i 



log(p2 + m^) . (2.38) 



The zero temperature contribution, Eq. (2.37) is clearly ultraviolet divergent. It can be 
evaluated by applying an ultraviolet momentum cut-off or using dimensional regular- 
ization. The finite temperature contribution, Eq. (2.38) is ultraviolet divergent as well. 
Because high-momentum modes at finite temperature are exponentially suppressed by 
a Bose-Einstein distribution function (as will be shown in the next section) and since 
Eq. (2.38) becomes equal to Eq. (2.37) in the limit of zero temperature, the divergences 
of Eq. (2.37) and Eq. (2.38) are the same. Hence the difference between those equa- 
tions, which is the normalized pressure is ultraviolet finite. One then finds that the 
pressure of a free scalar field md + \ dimensions is given by 

'P = -\\£ ^og{P'^ + m^) - riog(p2 + m2) = -}.[ \og{P'^ + m^) . (2.39) 

In the following section it will be explained how this expression can be computed. 
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To obtain the pressure of a free fermion field tlieory consider its Lagrangian density 
in Minkowskian space 

£ - (A (i'fd^ -m + nyo) ^ , (2.40) 

where // is a chemical potential for the fermion particle minus antiparticle number 
J if/'' if'. In Euclidean space this Lagrangian density becomes 

£- iff {-yodo + ijidi -m + jijo) if/ . (2.41) 

Performing the Gaussian path integral gives that 

log Z - ^y = log Det (-/5o + iyA -m+ /jyo) , (2.42) 

where the determinant is over the Dirac indices and a complete set of functions that 
satisfy anti-periodic boundary conditions in imaginary time. After going to momentum 
space it follows that 

log Z = 13V X log det {iyocbn + yiPi -m + fiyo) . (2.43) 

J IP] 

Evaluating the determinant over the Dirac indices and subtracting the divergent zero 
temperature contribution one finds that the pressure of a free fermion field in 4 dimen- 
sions is given by 

P^lf log( P^ + m^) , (2.44) 

JlP] 

where pa = &>„ + i/u. Like in the bosonic case discussed in the previous paragraph, this 
pressure is finite. In the following section it will be explained how this pressure can 
be calculated. 

2.4 Analytic calculation of sum-integrals 

As was discussed in the previous section one often has to evaluate sum-integrals in 
finite temperature field theory. In this section a method to perform these sum-integrals 
analytically will be discussed. The following section is devoted to the numerical eval- 
uation of sum-integrals. 

In order to calculate a sum-integral, one has to perform an infinite sum over Mat- 
subara modes, after which the integration over momenta has to be done. Such a sum 
over Matsubara modes can be obtained by using contour integration. Consider a par- 
ticular sum 

1 "' 

- 2] /(z - iojn) , (2.45) 

'"^ n=—co 

where a»„ = 2nnT as for bosonic fields. This expression can viewed as a sum over 
residues of some function which has simple poles located at z = /a>„. A sum over 
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Figure 2.1: Contour for summation formula. The black dots on the imaginary axis 
denote the poles of coth(jSz/2), while the other dots are possible poles of /(z) (in the 
left figure) and of /(z) + /(-z) (in the right figure). 



residues is equivalent to an integration around all the poles, which is sketched in the 
left-hand part of Fig. 2.1. Consider the function coth(j6z/2). It has only simple poles at 
z = ioJn, which all have residue 2/yS. So assuming /(z) has no poles on the imaginary 
axis, g(j) - co\h{]3zl2)f{z) has simple poles at z = iiOn with residue 2f{ia)n)/j3- This 
allows one to write the sum as the following integral 

^ !;/(. = -.) 4^ /dz/(z)coth(f), (2.46) 

'^ n=—co ^'~- 

where the contour C is depicted in the left part of Fig. 2.1. 

Now one can use that coth(y6z/2) = 1 + 2?i(z) = -1 - 2«(-z), where n(z) = 
l/[exp(/3z) - 1] is the Bose-Einstein distribution function. If lim^^oo z^(z) = 0, one 
can split C in two pieces along the imaginary axis and bring them together 



/-• r-ioo+£ 

I dzg(z) = I dz [g{z)-g{-z)] . 
The last equation can be used to write the sum as 

- y f{z = ico„) = -—. I dz [/(z) + /(-z)] [1 + 2n{z)] . 



(2.47) 



(2.48) 



If /(z) falls off rapidly enough at z ^ +oo it is possible to close the contour as is done 
in the right-hand part of Fig. 2. 1. The integral can now be calculated straightforwardly 
by summing over the residues. One should be aware that the contour in the right part 
Fig. 2. 1 goes clockwise, so one picks up an additional minus sign when applying the 
residue theorem to calculate the integral. 
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A frequently arising sum (see for example the gap equations calculated in Chapters 



4, 5 and 6) is the one over the propagator which using f{z) = l/(-z^ + w?,) with 



ap-p - p^ + m^ results in, 



- > ' —^ " T^ [l + ^"('^p)] ' (2-49) 



-y- 

' n 



where P = (po, p) with po = co^. After integrating over spatial momenta one finds the 
following important result 

By integrating Eq. (2.49) over ojp the sum-integral of a logarithmic function can be 
obtained. This sum-integral arises in the calculation of the pressure of a bosonic field 
theory (see for example Eq. (2.38) of the previous section and Chapters 4, 5 and 6). 
One finds 

^ log (P^ + m2) - J ^ [cOp + IT log (l - e-P'^")] + C , (2.51) 

where C is an infinite constant which is independent of ojp and temperature. Equa- 
tion (2.51) is ultraviolet divergent; all divergences arise from the integral over cop and 
the constant C. The high-momentum modes which depend on temperature are ex- 
ponentially suppressed so they do not give rise to divergences. In the limit of zero 
temperature a sum over Matsubara modes changes into an integration over pQ, hence 
Eq. (2.51) becomes in the limit of zero temperature 

j^log(p2 + m2).|i^^, + C, (2.52) 

which shows that all ultraviolet divergences of Eq. (2.51) are contained in the zero- 
temperature contribution. The pressure of a free bosonic field as is defined in Eq. (2.39) 
is hence ultraviolet finite and given by 

'p = -\£ iog(^' ^ ""'^ " "^ / ?l;7 ^°^ (^ ~ ''"^"") • ^^-^^^ 

In Chapter 7 a theory with fermions in the presence of a chemical potential will be 
discussed. In that case p^ = 2n{n + 1/2)7 + ijx - cOn + 7iT + i/u. Using the formalism 
presented above it can be shown that 
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where h+{(ji)p) = l/(expyS(a»p ±//) + 1). After integrating Eq. (2.54) over cOp the 
fermionic sum-integral of a logarithmic function can be obtained. This sum-integral 
arises in the calculation of the pressure of a fermionic field theory (see for example 
Eq. (2.44) in the previous section and Chapter 7). One finds 



|^_ log (P^ + „^) = J jgj L„ + T 2 log ( 1 + .-*"'*') 



+ C , (2.55) 



where C is a divergent constant which is independent of temperature and cOp. Like in 
the bosonic case, the high-momentum modes of the term which depends on temper- 
ature is exponentially suppressed an hence does not give rise to an ultraviolet diver- 
gence. All divergences are contained in the zero-temperature contribution, as a result 
the pressure of a free fermionic field in 3-1-1 dimensions (defined in Eq. (2.44)) is finite 
and given by 

If m = and ju = the temperature-dependent parts of the sum-integrals can be ob- 
tained exactly. The integrals which appear after summing over Matsubara frequencies 
can be evaluated with use of the following identities with rj = ±1 and n > 1 

n-l r°° °° °° „n 



noo n-l noo ^ ^ n 



s=0 .v=l 

r{nX{n) ifj] = +\ 

{\-2'-")Y{n)an) if 77 = -1 



(2.57) 



where r(«) is the gamma function which obeys: r(« -1- 1) = nr{n) and ^(n) is the 
Riemann zeta function. Some useful values of the Riemann zeta function are: ^(2) = 
;r2/6, ^(3) « 1.202 and ^(4) - ;rV90. Combining Eqs. (2.39), (2.51) and (2.57) gives 
the following result for the pressure of a massless bosonic spin-0 degree of freedom in 
3-1-1 dimensions 



7r2 
d2 _ ^ rr4 



n = -5;^;ogF' = -rV (2.58) 

By combining Eqs. (2.44), (2.55) and (2.57) the pressure of a massless spin- 1/2 fermionic 
field in the absence of a chemical potential in 3-f1 dimensions is obtained, which reads 

The last two equations can be used to obtain the pressure of QCD in the limit of 
infinite temperature. Due to asymptotic freedom the quarks and the gluons are effec- 
tively non-interacting in this limit. Moreover since the temperature is then much larger 
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than the quark masses, the quarks are effectively massless. So in the limit of infinite 
temperature the pressure of QCD is the sum of the pressure of eight massless gluons 
and Nf massless quarks. Since all eight gluons have two transverse polarizations their 
contribution to the QCD pressure in the infinite temperature limit is l&Ph- The Nf 
quarks all carry three colors, hence their contribution to the QCD pressure in the infi- 
nite temperature limit is 3NfPf. Adding the gluonic and quark contribution gives the 
QCD pressure in the infinite temperature limit which is !P = 37 n'^T^/ 90 for Nf = 2 
and P = l9n'^T^/36 for Nf = 3. In the confined phase at low temperatures one ex- 
pects the QCD pressure to be dominated by a gas of particles with lowest mass, which 
are the pions. Since there are three pions which are spin-0 particles, the pressure of 
QCD at low temperatures is something in the order of !P = n^T^ 130 (if one takes into 
account that the pions have a mass this pressure even becomes smaller, see Chapter 6). 
Hence the pressure divided by T^ at low temperatures is much smaller than in the in- 
finite temperature limit as also can be seen from Fig. 1.2, which display the results 
of lattice calculations of the pressure of QCD. In this figure, the infinite temperature 
limits calculated in this paragraph are drawn as well. 



2.5 Numerical computation of sum-integrals 

As was discussed in the previous sections, partition functions can be obtained by cal- 
culating sum-integrals. It was shown that these sum-integrals are typically ultravio- 
let divergent, however the difference between a sum-integral and an integral is finite 
because of the exponential suppression of high momentum modes. In the examples 
discussed in the previous section, the sum-integrals were obtained relatively straight- 
forwardly because the sum over Matsubara modes could be performed analytically. 
However in more complicated cases, which for example arise in Chapters 4, 5 and 6, 
an analytic result for this sum can no longer be obtained. In this section a method will 
be developed, which can be used to calculate the sum-integrals numerically. 

The sum-integrals which arise in calculating partition functions are typically of the 
following form (see for example Eq. (2.35)) 

^f{P^). (2.60) 

As was argued in the previous section, such sum-integrals are most often ultraviolet 
divergent. Hence, an immediate problem which arises when calculating a sum-integral 
numerically in the brute-force way (that is to perform the sum over Matsubara modes 
numerically and then to do the integration over spatial momenta numerically as well) 
are the ultraviolet divergences. To treat these ultraviolet divergences in a consistent 
way it is useful to split a sum-integral into a finite part containing the difference be- 
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tween a sum and an integral and a divergent part in the following way 
^ f{P)=^ f{P) + ^ f{P) 

In the following two subsection it is subsequently discussed how the integral and the 
difference between the sum-integral and the integral can be computed numerically. 

Computation of the integral 

The term j f{P) contains all possible ultraviolet divergences of £, f{P) and can be 
calculated for example by applying a cut-off or dimensional regularization. Numeri- 
cally a cut-off regularization is the easiest, however dimensional regularization is also 
possible numerically as was shown by Caravaglios (2000). In the cases considered in 
this thesis, the ultraviolet divergences D of the term proportional to j f{P) can al- 
ways be extracted analytically by considering the high-momentum behavior g{P) of 
f(P), i.e. D - jp g{P) (see also Blaizot et al. (2003)). The term which contains the 
divergences can now be written as 

J f{P) = D + j [f{P) - g{P)] , (2.62) 

where the last integral is finite and can easily be evaluated numerically using standard 
techniques like Gauss-Legendre integration. The rewriting of the divergences in terms 
of an integral prevents subtraction of two large quantities which, due to the finite ma- 
chine precision, can give rise to huge numerical errors. If f{P) depends explicitly on 
temperature (see Chapter 4, 5 and 6 for an example) D can have a temperature de- 
pendence as well, giving rise to temperature-dependent ultraviolet divergences. Such 
a divergence gives rise to renormalization problems, therefore a careful analysis is 
required. 

Computation of the difference between a sum-integral and an integral 

The difference between a sum-integral and an integral is most often finite because it 
turns out that its high momentum modes are exponentially suppressed (as was shown 
explicitly for two examples in the previous section). Since both terms in the differ- 
ence are divergent, it is not easy to obtain this difference numerically by using the 
expression in the second line of Eq. (2.61). However, using contour integration it is 
possible to derive an expression which does not contain divergent parts and hence will 
be suitable for numerical evaluation. After that, this expression has to be integrated 
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Figure 2.2: Contour C for summation formula. The black dots on the real axis denote 
the poles of cot(z/2), while the other dots indicate possible poles and cuts of /(z). 



over the spatial momenta which can be done using standard numerical techniques like 
Gauss-Legendre integration. 

For simplicity T = 1 is taken in this section. Results can however be easily 
generalized to any T by rescaling. Consider a function /(z) which is analytic for 
2n{N - 1/2) < Rez < 2n{M + 1/2), where A^ and M are integers. Using that cot(z/2) 
has poles at z = 2nn with residue 2, it holds that 



M 



I/*.' 



n=N 



2nn) = ^^i) dz/(z)cot(z/2) . 



(2.63) 



Because of the analyticity requirement on /(z) one is free to choose the contour C as 
long as it is closed, the points z - 2nn with N < n < M are included and other possible 
cuts and poles of /(z) are excluded. The contour C used in this section is displayed 
in Fig. 2.2, where C\ goes from a + r to a + re"^, C2 from a + re"^ to a, C3 from a to 
a + re""^ and C4 from a + re~"^ to a + r. Here a - 2n{N - 1/2) and r = 2t:{M -N +1). 
Now it can be used that 



cot(z/2) --/ 1 + 



e-'z - 1 



-/ 1 + 



e'^- 1 



(2.64) 



to obtain 



y f{z^2nn) = -— dz/(z) + — dz/(z) 



2njc 



1 



C1UC2 e '^-\ 



+ 



2njc 



1 



dzfiz)——-. (2.65) 
C3UC4 « ^ - 1 
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Since it is assumed that /(z) has no poles or cuts within the integration contour, 



r dz/(z)+ r 'dz/(z) = 0, 



r dzfiz)- r dz/(z) = 0. 



(2.66) 
(2.67) 



Combining the three equations above gives for the difference A/vm of a sum and an 
integral 



^N,M ^ y f{z = 2nn) - — dzfiz) = 

±Xi ^^ J2n{N-ll2) 

r dz/(z)-i-- r 



1 

2n 



1 



C1UC2 e '^ - 1 



. (2.68) 



Now consider the limit M ^ 00. For functions f{z) which grow slower than an expo- 
nential in the limit \z\ -^ 00, the contribution coming from integration along Ci and C4 
can be neglected. In that case 



A/v,oo - -:r- 
2n 



f 

Jo 



Ape-'^fia+pe''"^) 



1 



Jo 



exp(/pe '^) + 1 



+ dpe'V(a+pe'^)- 



(2.69) 



exp(-/pe'^) + 1 
A convenient choice is (^ = njl, which gives 

/ r"" 1 

Ayv.co = — dp [f{a - ip) - f{a + ip)] —— , (2.70) 

In Jo eP + 1 

where a = 2n(N - 111) and f{z) should be analytic for Rez > a. This formula is 
similar to the Abel-Plana formula (where a is taken to be 2nN, see for example Barton 
(1981)). Eq. (2.70) is, however, more convenient for numerical purposes due to the 
larger suppression factor of l/[exp(p) + 1] as compared to l/[exp(p)- 1] in the original 
Abel-Plana formula. If /(z) e R for z e R, Eq. (2.70) can be simplified to 



^N,oo 



-\ r 

^ Jo 



dp Im f{a + ip) 



1 



eP + 1' 



(2.71) 



The difference between a sum-integral and an integral can now be calculated using 
the following expression 



//<-Jl^- 



(2.72) 
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The term A_oo,oo can be computed using Eq. (2.70) although not necessarily with 
A'^ - -oo, since to apply that equation f{P) should satisfy the analyticity requirements 
discussed below Eq. (2.70). To illustrate the computation of Eq. (2.72) two examples 
will be given below. 

In order to obtain the pressure of a free scalar field theory (see Eqs. (2.39) and 
(2.51)) the term X,log(/'^ + m^) has to be computed. To evaluate this expression nu- 
merically one takes f{z) = log(z^ + ojI). Clearly the sum over all Matsubara modes 
(z = Inn) diverges as does the integral over z- However, the difference between the 
sum and integral is finite. To calculate this difference it is not allowed to use Eq. (2.70) 
with N = -oo, because /(z) has a cut at Rez = 0. One therefore has to split the 
difference into parts, for example in the following way 

A_oo,oo = A_oo,-i + A-1,1 + Ai,oo . (2.73) 

Because /(z) is even in z, A_oo,-i = Aioo- The term Aioo can be calculated numeri- 
cally using Eq. (2.71). The term A_ij can be obtained numerically using the explicit 
expression for the difference between a sum and integral. As a result 

1 r" 

A^eo.oo = 2Ai,eo + log(a;^) -J^\<^z log(z^ + ^p) • (2-74) 

Now X, log{P^ + m^) can be obtained by integrating Eq. (2.74) over momenta using 
Eq. (2.72). 

In Chapters 4, 5, and 6 of this thesis, next-to-leading 1/A'^ corrections to the pres- 
sure are investigated. To obtain these corrections one has to calculate a sum-integral 
of the following form 

F=^f{P), (2.75) 

where f(P) = log [c -i- I(P, m)], c is some constant and 

liP, m) = { -^^ ^ 2 • ^2.76) 

Jg Q^ +m^{P+ QY + m2 

The function f{P) is even in po, has a cut for Re(/7o) - and depends explicitly on 
temperature. In order to use Eq. (2.71) to calculate the difference between a sum and 
an integral one has to split A-oo,oo as follows 



o,oo ^ 2Aa?,oo + V fipo = Inn) - — \ dpofipo) , (2.77) 

.._^, , 27r J-2MN-1/2) 



where N > I because of the cut in f{P). The term Aa?_oo can be obtained numerically 
using Eq. (2.71). It was checked numerically for several examples that changing A'^ has 
no effects on A_oo,oo as expected. After calculating A-oo,oo the integration over p can 
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be done straightforwardly to obtain F using Eq. (2.72). It was observed numerically 
that the difference of a sum and an integral is dominated by the low-momentum modes 
which can be understood as being due to the exponential suppression of the high- 
momentum modes at finite temperature. Hence after integration over momenta p one 
obtains a finite result for F as expected. This conclusion applies to difference of sum- 
integrals and integrals of more general f{P) as well. 



Chapter 3 

The effective potential and the IjN expansion 



As was discussed in the previous chapter, thermodynamical quantities can be derived 
by calculating the partition function. In this chapter, it will be made clear how such 
a partition function can be obtained by locating the extremum of a so-called effective 
potential. The derivation of this effective potential will be reviewed in the first section 
of this chapter. Then the \wge.-N approximation, which can be used to investigate non- 
perturbative phenomena, will be discussed. Thereafter, the auxiliary field method will 
be studied. Using the auxiliary field method the large-A'^ approximation of the models 
examined in this thesis can be obtained systematically. Finally it will be explained why 
some effective potentials expressed in terms of auxiliary fields can have temperature- 
dependent ultraviolet divergences outside the minimum. 



3.1 The IPI effective action and the effective potential 

The one particle irreducible (IPI) effective action is a useful tool for calculating par- 
tition functions. It is especially preferred if a certain field can obtain a vacuum ex- 
pectation value. This happens for example in a theory with spontaneous symmetry 
breaking. In ordinary perturbation theory one expands around the trivial vacuum, for 
which the vacuum expectation values of fields vanish, say {(p) - 0. In the IPI effective 
action approach however, one expands around the true vacuum. This means that the 
IPI method allows for a nonzero vacuum expectation value of 0, (0) - 4>- This vac- 
uum expectation value can be found by minimizing the IPI effective action. Because 
in the IPI method one perturbs around the true vacuum it is more advantageous to 
use this method over ordinary perturbation theory. Moreover it allows one to resum 
whole classes of Feynman diagrams and to investigate non-perturbative physics. More 
details concerning the derivation of the effective action given below can be found in 
quantum field theory textbooks, for example Weinberg (1995), Peskin and Schroeder 
(1995) and Zinn- Justin (1996). 
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Formal definition of the efl'ective action 

Consider a general scalar field theory with Lagrangian density X[0] in Euclidean 
space-time. The partition function of this theory in the presence of a source term J 
is given by 



Z[/] = \ Dcp exp -S [cp] - fd''xJ{x)(l>{x) 



(3.1) 



here S[(l>] - j d'^xJl[(p] denotes the classical action. From this partition function one 
can compute the generating functional for connected Green's functions which is 

W[/] = -logZ[7] . (3.2) 

These Green's functions are obtained by functional diff"erentiating W[J] with respect 
to /. For example the vacuum expectation value of (p of the theory in the presence of a 
source term is given by 

^^ = {<f>(x)}j ^ 4>(x) . (3.3) 

Similarly, the connected two-point correlation function can be found as follows 



6^W[J] 
5J{x)6J(y) 



= {Tcp{x)<p{y))j - {<P{x))j{<p{y))j ^ D{x, y) . (3.4) 



The IPI effective action T[4>] is the generating functional of IPI diagrams. An IPI 
diagram is a diagram that is still connected after cutting one internal line. The IPI 
effective action is the Legendre transform of W[J], 



-/d^ 



r[0] = W[J] - d^x/(x)0(x) . (3.5) 

The effective action is a function of only and not of /. The source term / in Eq. (3.5) 
is to be chosen in such a way that the vacuum expectation value of (p in the theory with 
a source term will become equal to (/>. The IPI effective action is, unlike the classical 
action, an action which contains all contributions arising from quantum fluctuations. 
Extremizing the effective action with respect to cp gives 



X) 



J 6J(y) 6(p{x) J 6(p{x) 

Hence at an extremum of the effective action the source term has to vanish, in other 
words / = 0. The extremum will be denoted by ^(x), 

-^^ = . (3.7) 
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From the definition of the effective action, Eq. (3.5), it follows that the effective action 
at the extremal value is equal to 

rM = -logZ[0]. (3.8) 

This important equation shows that the partition function can be obtained by calculat- 
ing the minimal value of the effective action. 

It often happens that the vacuum solution of Eq. (3.7) is translational invariant, 
which implies that (fi{x) is space(-time) independent. However, translational non- 
invariant solutions for (fi{x) are sometimes possible. Instantons which for example 
arise in the CP^^^ model discussed in Chapter 5, are translational non-invariant solu- 
tions. But due to their dependence on space(-time), these translational non-invariant 
solutions often have a larger action than the translational invariant solution. If ip{x) is 
translational invariant it is useful to define the effective potential 'V which is minus the 
effective action divided by the volume V of the space 

T(0)--^r[0], (3.9) 

where the sign of the effective potential is chosen in such a way that the extremal 
value of the finite temperature effective potential becomes equal to the pressure. The 
extremal value of the effective potential is given by 

^(^) = -ilogZ[0]. (3.10) 

At finite temperature the volume of the space 5 x R'^ is PV. Therefore at finite tem- 
perature the extremal value of the effective potential is given by log Z[0]/fiV which is 
equal to the pressure. In the rest of the thesis this equation will be used to calculate the 
thermodynamical quantities. 

Next to the IPI effective action, a 2PI effective action (Cornwall et at, 1974) and 
even more general nPI effective actions exist as well. The 2PI effective action is the 
Legendre transform of the generating functional for connected Green's functions in 
the presence of a source term / for the field (p and a source term K for the two-point 
function. The 2PI formalism is very useful for out-of-equilibrium quantum field theory 
calculations (see for example Berges (2005)). This is because unlike the IPI method, 
the 2PI formalism lacks the so-called secularity problem, which causes the pertur- 
bation series of a time-dependent quantity to diverge at late times (see for example 
Arrizabalaga (2004)). Moreover, it is very natural to introduce Gaussian initial density 
matrices in the 2PI formalism. Since in this thesis all results are obtained for equi- 
librium situations, the 2PI formalism will not be discussed further. Although it is of 
course also possible to perform equilibrium calculations using the 2PI formalism. 
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Perturbative calculation of the effective action 

In practice, the effective action has to be calculated in perturbation theory. Following 
the method of Jackiw (1974), the field (^(x) is replaced by the sum of its vacuum 
expectation value 0(x) and a quantum fluctuating field ;^'(x). In this way it is possible 
to perturb around the true vacuum ^(x). Taylor expanding the action and the source 
term around ^(x) gives. 



S{^+X\ + ^ d^^ J{x) \x{x) + 4>{x)] = S[4>] + Jd^x J{ 

^jd'xd 



x)0(x) 



- fd^wd'^xd^ 
3! j 



+ /(x) 



Xix) + 



''yxix)G(x,yr'x(y) 



S^S[(f>] 



6(p{w)64>{x)S(p(y) 



xiw)x(x)x(y) + ■ ■ • ■ 



where the bare inverse propagator is given by 

s2 



G(x,y) = 

6<Pix)64>(y) 



(3.11) 



(3.12) 



Using the expansion of the action, Eq. (3.11), and the definition of the effective action, 
Eq. (3.5), it follows that the effective action obeys 



rm^sm+r2m 



(3.13) 



So the effective action r[0] is equal to the sum of the classical action S [cp] and the 
quantum corrections r2[0] which are given by 



r2[^] - -log fox exp -S'lx]- fd'^xf 



'{x)xix) 



(3.14) 



In the expression for Y2 the action 5" |;kf] is equal to 



S'M=^ Jd''xd''yx{x)G{x, 

+ — fd^'wd^xd^v - 
3! J ^ i 

and the current /' is equal to 



yrxiy) 



{w)6cp{x)6(P(y) 



(3.15) 



x{w)x{x)xiy) + ■■■ 



J'{x) = 



ss{cp^ 



(X) 



(5r[0] _ 6Y2m 
^=4, Scpix) 6(pix) 



(3.16) 
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The quantum corrections to the effective action can be obtained by calculating 
r2[0] in perturbation theory by summing Feynman diagrams. In order to find out 
which kind of diagrams contribute to r2[0], consider Z[0]. In perturbation theory, Z[0] 
is equal to 1 plus the sum of all closed loop Feynman diagrams times the contribution 
of the single loop (see also Sec. 2.3). The quantity W[0] - -logZ[0] is equal to the 
sum all connected closed loop Feynman diagrams. And r[0 = 0] = -W[/'] is equal to 
the sum of all one particle irreducible Feynman diagrams, for a proof see for example 
Zinn-Justin (1996). Here /' is a source term that is chosen in such a way that the 
corresponding vacuum expectation value vanishes. As one might expect the source 
term /' of Eq. (3.16) is tuned automatically in such a way that the vacuum expectation 
value X of the fluctuating x field vanishes, for a proof see Jackiw (1974). Hence Y2 is 
equal to the sum of all closed loop IPI diagrams with bare propagator G and vertices 
as given by the shifted action S' . It is important to realize that although a tadpole term 
arises in the shifted action, tadpole diagrams do not contribute to T2- The tadpoles 
form part of the current /', see Eq. (3.16), which forces the vacuum expectation value 
of the X field to vanish. 

If the contribution of higher order IPI diagrams to Y2 is suppressed, for instance 
due to some small coupling constant, the main contribution to Y2 will arise from the 
single closed loop. In such a case the effective action can be approximated by 

r[0]^5[0] + ^logDetG, (3.17) 

where the log Det term arises from the Gaussian integration over the quantum fluctua- 
tions X- As will be argued in the following sections, Eq. (3.17) can be used to obtain 
thermodynamical quantities to next-to-leading order in 1 jN, where A'^ is the number of 
fields in the theory. 

An explicit calculation of the effective potential in a scalar A(p'^ field theory has 
been performed by Coleman et al. (1974) in an expansion in small X. This work was 
generalized to finite temperature by Dolan and Jackiw (1974). Although the discussion 
of the effective potential in this section was limited to scalar fields only, the effective 
potential for a theory with gauge fields and fermions is obtained analogously. In this 
thesis finite temperature effective potentials will be calculated not in a perturbative 
expansion in the coupling constant, but rather in an expansion in 1 IN. This 1 IN ex- 
pansion will be explained next. 

3.2 The 1 IN expansion 

To calculate a quantity in an interacting field theory it is often necessary to make some 
kind of approximation. Expanding in the coupling constant i is a widely used method. 
This however has the drawback that it only works for small couplings. Moreover, 
contributions which are non-analytic in A, like exp(-l//i), will not be found in an 
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expansion around A -0.\i the coupling constant is large, which happens for example 
to be the case in QCD at low energies, perturbation theory is not reliable. Hence in 
this non-perturbative regime one has to use other approaches like for example the 1 jN 
expansion, lattice discretization, Dyson-Schwinger or renormalization group methods. 
In the large-A'^ method the expansion parameter is the number of fields N. In order 
to perform this expansion one has to change the symmetry group of the theory under 
consideration. To illustrate this point, consider the 0(4) linear sigma model (which 
is discussed in more detail in Chapter 6). This model has the following Lagrangian 
density 

£=\ (^^^'p)^ + ^ [Mi-^fltf ' / = 1 . . . 4 , (3.18) 

where Aj, is a bare 4-point coupling constant and f„^h is the bare "pion decay" constant. 
To apply the 1 /N expansion in the 0(4) linear sigma model it has to be generalized to 
the 0(A'^) linear sigma model. The Lagrangian density of the 0(A'^) linear sigma model 
is given by 

X = ^(VO' + ^(0,0,-A^/,y' , /=l...A^. (3.19) 

As one can see from Eq. (3.19), the coupling constants Ah and fj^j, have been rescaled 
by factors of N. This redefinition of the coupling constants is important when applying 
al/N approximation, since it assures that the relative strength of the interactions is not 
changed when varying N. The coupling constants are rescaled in such a way that the 
action naturally scales with A'^. In the 0(A'^) linear sigma model the diagrams which are 
contributing to the pressure at leading order are the bubble diagrams and all possible 
insertions of bubbles, see for example Jackiw (1974). These diagrams are called daisy 
and superdaisy diagrams as well and are displayed in Fig. 3.1. The chain diagrams 
which are displayed in Fig. 3.2 contribute to next-to-leading order in l/N. All these 
diagrams can be resummed using the auxiliary field method which will be discussed 
in the next section. In this way it is possible to obtain an expansion which is really in 
1 /N and not for example in A/N. This allows one to investigate the non-perturbative 
large A behavior of the model. 

Unfortunately the situation in QCD is not that simple. The obvious generaliza- 
tion is to consider an SU(A^c) gauge theory and to take the number of colors, Nc, as 
an expansion parameter. The coupling constant ^^ has then to be rescaled to g^/Nc- 
As was shown by 't Hooft (1974), only planar diagrams contribute to leading order in 
the 1 /Nc expansion. Non-planar diagrams with only one handle contribute at next-to- 
leading order. Just like in the 0(A'^) linear sigma model an infinite number of diagrams 
contribute to the QCD pressure at leading and next-to-leading order in Nc- But unfor- 
tunately until now no one has found a way to resum all these diagrams, the auxiliary 
field method discussed in the next section for example does not work for QCD. 

Another possibility is to give all quark flavors the same mass and expand in the 
number of flavors Nf. Using this method Moore (2002) calculated for zero quark 
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Figure 3.1: Some examples of bubble diagrams in the OCA'^) linear sigma model. An 
infinite number of these types of diagrams with all possible insertions of bubbles con- 
tribute at leading order in l/N. The lines in these diagrams denote the bare propaga- 
tors, which can be read off directly from the Lagrangian density of the 0(A'^) linear 
sigma model, Eq. (3.19). 







Figure 3.2: Some examples of chain diagrams in the 0(AO linear sigma model. An 
infinite number of these types of diagrams with all possible insertions of bubbles con- 
tribute at next-to-leading order in 1 /N. 



mass the full non-perturbative large Nf contribution to the QCD pressure up to next- 
to-leading order. Ipp and Rebhan (2003) extended this study to non-zero chemical 
potential. Unfortunately this approach gives only limited insight in QCD, since large 
Nf QCD is not asymptotically free and hence behaves completely diff'erently from real 
QCD. 

A hint that in the future methods may become available which can be used to 
perform analytic calculations in the non-perturbative regime of QCD comes from the 
anti-de-Sitter conformal field theory (AdS/CFT) correspondence. Maldacena (1998) 
conjectured that N - A supersymmetric U(AO Yang-Mills theory in 4 dimensions 
(which is a conformal gauge field theory) is equivalent to a 10-dimensional type JIB 
string theory on AdSs x 5 ^ . It turns out that if the supersymmetric Yang-Mills theory is 
strongly coupled, the corresponding string theory is weakly interacting. So by apply- 
ing perturbation theory to this 10-dimensional string theory one can obtain insight in 
the non-perturbative behavior of the corresponding 4-dimensional gauge theory. One 
is of course eager to find an analogous result for real QCD in 4 dimensions. 
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3.3 The auxiliary field method 

The 1 IN expansion can in some theories be systematically performed by introduction 
of auxiliary fields. These auxiliary fields will have no dynamics so they will not influ- 
ence the physical content of the theory. They will be shifted by a constant in order to 
obtain a Lagrangian density which is quadratic in the A'^ original fields. This allows for 
Gaussian integration, after which the \IN expansion follows naturally. 

As an example consider the O(A^) model defined in Eq. (3.19). Its partition function 
is given by Eq. (3.1). An auxiliary field a can be added to this model by changing the 
Lagrangian density as follows 

Z^L+ ^a^ . (3.20) 

ZAh 

The original partition function times an infinite constant (which drops out when calcu- 
lating physical quantities) is recovered after integration over the a field. Since a field 
is nothing else than an integration variable of the path-integral, it is always possible to 
shift this field by a constant. If the a field is shifted according to 

a^a-!^{0icl>i-Nfl,), (3.21) 

the Lagrangian density becomes quadratic in the fields, 

X - ^ {d,<l>if - '-a [Mi - Nfl,) + ^a2 . (3.22) 

It is now possible to perform the integration over the scalar fields, resulting in the 
following eff"ective action (not to be confused with the IPI effective action, here a is 
still a field which should be integrated over) 

5eff M = ^Tr log (-52 - ia) + f r dr fd^x iiflj,a + ^ J . (3.23) 

The auxiliary field can obtain a vacuum expectation value a which is assumed to be 
translational invariant. Hence to obtain for example the pressure, it is necessary to 
calculate the effective potential and minimize it with respect to a. The leading order 
term of this effective potential is just the classical effective action, Sss{a) and is pro- 
portional to A'^. By expanding ^eff around its vacuum expectation value the quantum 
corrections to the effective potential can be obtained. As will be shown in detail in 
Chapters 4, 5 and 6, the leading term of these corrections is the logarithm of a deter- 
minant of the propagator of the quantum fluctuations which is proportional to A'^", see 
also Eq. (3.17). The contributions which arise from the IPI diagrams of the quantum 
fluctuations are of order 1 IN and higher because the vertices of these diagrams turn 
out to have 1 / VA^ suppressing factors. 
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3.4 Temperature-dependent ultraviolet divergences^ 

The finite temperature effective potential expressed in terms of the vacuum expecta- 
tion value of the auxiliary field discussed in the previous section is calculated explicitly 
for different models in Chapters 4, 5 and 6. In all these models it turns out that the 
effective potential contains temperature-dependent ultraviolet divergences. However, 
as follows from the calculations, these divergences become temperature-independent 
at the minimum. These temperature-dependent divergences make renormalization of 
the effective potential outside the minimum impossible because it will require a renor- 
malization prescription which is different for every temperature. Such a prescription 
would lead to arbitrary temperature dependence and hence is meaningless. 

At first sight it might seem surprising that these temperature-dependent ultraviolet 
divergences arise at all, since one expects that temperature cannot influence the ul- 
traviolet behavior of a theory due to the exponential suppression of high-momentum 
modes of finite temperature contributions. To see why it is nevertheless possible that 
an effective potential can have temperature-dependent divergences outside its mini- 
mum, consider as an example the 0{N) linear sigma model defined in Eq. (3.19). For 
other models the argument is similar. The effective action of the 0(A'^) linear sigma 
model expressed in terms of the auxiliary field, 5'eff[Q'] is given in Eq. (3.23). The 
effective action in the presence of a space-time independent source // for a is defined 
by 

5eff[/,af] ^ ^effM + // I dr | d^xa{x) . (3.24) 

The quantity Z[/]r is defined here as the partition function of the theory with action 
SssVJ,ct\ at temperature T. Now logZ[7 - Olr/ySV is equal to the pressure of the 
linear sigma model. One therefore expects that this quantity should be well-defined 
and that it contains no temperature-dependent divergences. So one can assume that 
logZ[0]r does not contain temperature-dependent divergences. As a result, logZ[7]7-, 
contains also no temperature-dependent divergences as long as J is independent of 
temperature. The reason for this is that since both / and /^^ are proportional to a in 
the Lagrangian density, any source / can be incorporated in /^^ by the redefinition 
NflJl + / -. NflJl. 

The effective potential in the minimum is equal to logZ[0]r. On dimensional 
grounds, the divergences of the minimum of the effective potential should depend on 
fn^b- For example in Chapter 6 it is found that the effective potential contains terms 
proportional to the divergent factor A^/^^, where A is an ultraviolet momentum cut- 
off. Since any / can be incorporated in /^^, varying / will change the divergences of 
logZ[/]r. Hence the divergences of logZ[/]7- will depend on /. 



'This section is based on: Thermodynamics of the 0(N) non-linear sigma model through next-to- 
leading order in 1/N, H.J. Warringa, proceedings of tlie SEWM 2004 meeting, World Scientific (2005). 
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The effective potential is given by 'V{m^) - log Z[J]t I P^ - Jm^, where a = 
im^ is the vacuum expectation value of the a field. That this vacuum expectation 
value is purely imaginary is proved in Chapter 4. In the expression for the effective 
potential, / is defined implicitly, since J is the source term which gives the a field 
the vacuum expectation value ini^. If the temperature is varied, the implicitly defined 
J has to change in order to keep the vacuum expectation value of the a field fixed to 
im^. Therefore the implicitly defined / has to depend on temperature. In the previous 
paragraph it was argued that the divergences of log Z[J]t depend on /. Since log Z[J]t 
is part of the effective potential, the divergences of the effective potential will depend 
on temperature for any non-zero /. 

As is shown in Eq. (3.10), the source term / vanishes at an extremum of the effec- 
tive potential. This holds for any temperature. So using the arguments in the previous 
paragraphs temperature-dependent divergences do not arise at the minimum of the ef- 
fective potential, but only outside the minimum. Since physical quantities like the 
pressure, entropy density and energy density are to be calculated at the minimum of 
the effective potential they do not suffer from the renormalization problems caused by 
temperature-dependent divergences. 



Chapter 4 

Thermodynamics of the nonhnear sigma model in 
d=l+l 



The thermodynamics of the 0(N) nonUnear sigma model is studied in this chapter in 
the l/N approximation. The effective potential, from which one can derive thermo- 
dynamical quantities, will be calculated to next-to-leading order in l/N. It is found 
that the effective potential contains temperature-dependent divergences which become 
independent of temperature at the minimum. As a result it turns out that the effective 
potential can only be renormalized meaningfully at the minimum, hence thermody- 
namical quantities which are obtained at the minimum, can be rendered finite consis- 
tently. The renormalized pressure will be calculated to next-to-leading order in 1 /N. 
Furthermore in intermediate steps of the calculation, thermal infrared renormalons 
are encountered which potentially can give rise to ambiguities in the end result. It 
will be shown that the correctly renormalized pressure is free from these ambigui- 
ties. This chapter is based on: Thermodynamics of the 0(N) nonlinear sigma model in 
(l+l)-dimensions, J.O. Andersen, D. Boer and H.J. Warringa, Phys. Rev. D69 076006, 
(2004). 

4.1 Introduction 

The 0(A'^) nonlinear sigma model in 1 -i- 1 dimensions has been studied extensively 
at zero temperature as a toy model for QCD. It is a scalar field theory which classi- 
cally is scale invariant. Due to renormalization of the quantum corrections a scale is 
introduced. As was shown by Polyakov (1975) the nonlinear sigma model is asymptot- 
ically free like QCD. It is renormalizable both perturbatively and in the l/N expansion. 
Moreover, as is discussed in more detail in Chapter 5, it contains instanton solutions 
ior N = 3 only. In that chapter it is also shown that the 0(3) nonlinear sigma model is 
equivalent to the CP^ model. 

Unlike the nonlinear sigma model in more than two dimensions, where the theory 
is no longer renormalizable, there is no spontaneous symmetry breaking of the global 
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0{N) symmetry for any value of the coupling constant. This reflects the theorem of 
Mermin and Wagner (1966) and Coleman (1973), which forbids spontaneous break- 
down of a continuous symmetry in a homogeneous system in one spatial dimension at 
any temperature. The nonlinear sigma model suffers from infrared divergences in per- 
turbation theory, because the scalar fields remain massless in that case. It was conjec- 
tured by Elitzur (1983)' and shown by David (1981) that the infrared divergences can- 
cel in 0(A^)-invariant correlation functions. However, like in QCD, a non-perturbative 
mass-gap is generated dynamically due to the interactions. In contrast to QCD where 
the spontaneous breakdown of chiral symmetry breakdown is responsible for non-zero 
masses, the non- vanishing mass of the scalar fields in the nonlinear sigma model is not 
related to symmetry breaking. In the \wge.-N limit, which is equivalent to summing 
all so-called daisy and super-daisy graphs (see Chapter 3), it turns out that the mass 
at zero temperature is given hy m - ^ exp {-2n/g^), where g is the coupling constant 
and fi is the renormalization scale. The mass is non-analytic in g which explains that 
it vanishes in perturbation theory. 

Dine and Fischler (1981) have investigated the nonlinear sigma model in 1 -i- 1 di- 
mensions at finite temperature. They calculated the free energy in perturbation theory 
and to leading order in the laige-N limit. In the weak-coupling expansion, they showed 
that the two-loop contribution to the ideal gas vanishes and that the three-loop contri- 
bution is infrared finite. In this chapter the analysis of Dine and Fischler (1981) will 
be extended to next-to-leading-order (NLO) in the 1 /N expansion. To obtain the ther- 
modynamical quantities the effective potential will be calculated. At zero temperature, 
the effective potential (or equivalently the Gibbs free energy) has been investigated at 
this order by Root (1974) and Biscari et al. (1990). 

There are several reasons to calculate the 1/A'^ corrections to the thermodynam- 
ics. Since the nonlinear sigma model is one of the few field theories in which non- 
perturbative calculations can be done (partly) analytically, it is very useful to check 
whether the 1 /N expansion is converging. Furthermore renormalization at finite tem- 
perature is non-trivial, since next to the renormalization scale an additional scale is 
introduced. Normally one expects that the ultraviolet behavior of a theory will not de- 
pend on temperature because the high-momentum modes of thermal contributions are 
exponentially suppressed due to the Bose-Einstein distribution factor. However, it is 
found that the next-to-leading order effective potential contains temperature dependent 
divergences. Fortunately it turns out that the effective potential can still be renormal- 
ized at the minimum, but it then contains an ambiguity which is related to a so-called 
infrared renormalon (David, 1982). It will be shown that this ambiguity vanishes when 
calculating the pressure. 

It was shown by David (1982) that the 1 4- 1 dimensional nonlinear sigma model 
contains infrared renormalons. Renormalons (see Beneke (1999) for a review) give 
rise to ambiguities in resummming a perturbative series in the coupling constant g. 



'This article was submitted in 1978. 
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The size of these ambiguities may give insight in the importance of non-perturbative 
corrections. In QCD these infrared renormalons appear as well, which in principle 
might teach us about the size of non-perturbative physics in QCD. Because the OCA'^) 
nonlinear sigma model can be accessed both non-perturbatively (via the 1 IN expan- 
sion) and perturbatively (in an expansion in g) the nonlinear sigma model is widely 
used (David, 1984, 1986; Beneke et al, 1998) to study these infrared renormalons. In 
this chapter such ambiguities are encountered when calculating the effective potential. 
However it turns out that they vanish when calculating the pressure. Furthermore it is 
investigated how these renormalons depend on temperature. It is found that thermal 
infrared renormalons in the nonlinear sigma model are independent of temperature at 
the minimum of the effective potential, however outside the minimum the renormalon 
residue become temperature-dependent while the position of the renormalon pole is 
not affected by temperature. Thermal renormalons have been studied before by Loewe 
and Valenzuela (2000) in 0"* theory in 3-1-1 dimensions. In this theory, one deals with 
ultraviolet renormalons only and thus it resembles QED rather than QCD. They found 
that the residues of the ultraviolet renormalon poles in the Borel plane are temperature- 
dependent while the position of the poles are not. 

The Gross-Neveu model, which is a 1 -i- 1 dimensional toy model which includes 
also fermions was recently studied by Blaizot et al. (2003) at finite temperature to next- 
to-leading order in the l/N expansion as well. While there are similarities between 
this model and the nonlinear sigma model, such as dynamical mass generation and 
asymptotic freedom, no problems related to renormalization of the effective potential 
was encountered by Blaizot et al. (2003). 

This chapter is organized as follows. In Sec. 4.2 the Lagrangian density of the 
nonlinear sigma model will be given. Thereafter, in Sec. 4.3 it will be discussed how 
the pressure could in principle be obtained by a weak-coupling constant expansion. 
Sec. 4.4 is devoted to the calculation of the effective potential to next-to-leading order 
in I /N. The zero-temperature case will be discussed first, after which the derivation 
of the finite temperature effective potential will follow. In Sec. 4.5 the pressure for the 
0(A'^) nonlinear sigma model will be presented. Some approximations for the pressure 
at high temperatures will be discussed in Sec. 4.6. Thermal infrared renormalons are 
investigated in Sec. 4.7. A summary and conclusions are given in Sec. 4.8. 

4.2 The nonlinear sigma model 

The nonlinear sigma model is a scalar field theory with an 0(A'^) symmetry. It is de- 
scribed by a Lagrangian density which only consists of a kinetic term and a constraint 

£=^df,cf>id^cl>i, cl)i{x)Ux)=Nlgl, i=\...N, (4.1) 

where gh is the bare coupling constant. The constraint forces the fields to live on 
a N - 1 -dimensional hypersphere, which causes the interactions between the fields. 
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Since gy has no dimension md - 1 + 1, the nonlinear sigma model is renormalizable. 
As one can see from Eq. (4.1), the classical 0(A'^) nonlinear sigma model is scale 
invariant. Due to renormalization of the quantum corrections a scale n is introduced in 
the theory. As will be shown explicitly in Sec. 4.4, the 0(A'^) nonlinear sigma model is 
asymptotically free. For N = 2 the nonlinear sigma model is a free field theory. This 
can be seen by choosing 0i - y2cos{-&) /gi, and 02 - y2.sm{-&) /gi,. The constraint is 
satisfied automatically by this choice, and the Lagrangian density turns into 

X = ^{^^^&f , (4.2) 

Sh 

which is the Lagrangian density of a free field theory. 

Thermodynamical quantities of the nonlinear sigma model can be derived from the 
partition function which is given by 



r r r*® r 

Y\^4>iY\^{N/8l-(f>icf>i)exp - dr dxX 
J .^j ^ \_ Jo J . 



(4.3) 



In the following section thermodynamics in the weak-coupling expansion will be dis- 
cussed. The remaining sections of this chapter will deal with the non-perturbative 1 /N 
approach. 

4.3 Thermodynamics in the weak-coupling expansion 

The weak-coupling expansion can be performed by writing the scalar fields as cp = 
(tti, 7r2, . . . , ttat-i, o"). Due to the constraint it holds that o"^ = N/g^ - ntJii with / = 
1 ... A'^ - 1 . Integrating over the cr field results in the following partition function 



Y\ !Dni Y\ & {N/gl - nin) exp - dr dx £,s{n) 
i=\ X V Jo J 



(4.4) 



where 7?(x) is the step function. This step function reflects the fact that due the original 
constraint on the (f) fields, :7r,7r,- never can be larger than N Ig^. The effective Lagrangian 
density Xeff can be found by inserting the expression for cr in terms of t: in the original 
Lagrangian density. As a result one finds 

Xeff (^) - \d^7:iifni + ^ ^""^"f - l6\0) log (N/gl - mn) , (4.5) 
2 2 N- gljTim 2 V o ; 

where the term proportional to 5^(0) arises from integration over the delta function 

6{N/gl - cf>i4>i). 
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For small values of g^ the 7?(x) function is only vanishing when ;r(x) is large. Since 
large values of tt give a small contribution to the partition function one can approximate 

{^{N Ig^ - niTTi) K 1 which gives 

.N-l r ^/3 

(4.6) 



I ]~[D7r,exp - I dr I dxXeffW 



li gf - it can be seen from Xeff that the nonlinear sigma model contains N-l non- 
interacting TT fields. Hence the leading contribution to the pressure of the nonlinear 
sigma model is that of a gas oi N - I non-interacting scalar fields which is 

9 = —j log(^2) ^ (^ _ i)_7^2 ^4 7^ 

Using the effective Lagrangian density Eq. (4.5), Dine and Fischler (1981) showed 
that the order g^ correction to the pressure vanishes. In their calculation, effects of the 
step-function were not taken into account. However, using the the 1 /N expansion, one 
does find corrections to the free pressure as will be discussed in detail in the following 
sections. 

4.4 The effective potential 

By introducing an auxiliary scalar field a, the constraint on the fields can be ex- 
pressed in terms of a path integral over the a fields in the following way 



Y\ 5 (cPicPi - N/gl) = C JDaexp J^a (cja^Pi - Nig]) 



(4.8) 



where C is a temperature and gb independent normalization constant. As a result one 
finds the following Lagrangian density which is equivalent to the original nonlinear 
sigma model Lagrangian density 



-df,4>id^'(f>i - -I 



X = ^df,4>i&'4>i - ^ai4>i<Pi - N/gl) + Ji4>i , (4.9) 



where Jt is a source term for the scalar field. It is useful to keep this source term 
for a moment in order to determine the scalar field propagator. By calculating this 
propagator to next-to-leading order in 1 /N one can determine the physical mass of the 
scalar field to next-to-leading order. Any renormalization prescription should give rise 
to a finite physical mass. The expression for the physical mass will be used later in 
order to renormalize the so-called gap equation to next-to-leading order in 1 /N. 

Using the auxiliary field method the Lagrangian density becomes quadratic in the 
scalar fields. Therefore one can perform the integral over these fields. This yields the 
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following exact expression for the partition function of the nonlinear sigma model 

Z= fDaexpC-^eff) , (4.10) 

where the effective action is given by (Novikov et at, 1984) 

Ses - ^Tr log (-52 - ia) + /-^ f a + ^ f f Ji{x){x\ J . |y)7,(j) . (4.11) 



IX ^y 



In the last expression the trace is over a complete set of functions that satisfy the 
boundary conditions for scalar fields at finite temperature. 

The vacuum expectation value of the a field, {a) is purely imaginary since from 
5'eff(Q')* - 5'eff(-Q') and changing the integration over a to -a it follows that 



<«>* 



r£)aaexp [-5eff(-a)] = - f^Jaaexp [-5eff(a)] = -(a) . (4.12) 



It is very likely that (in absence of a source term) the vacuum expectation value of a is 
space-time independent because the Euclidean space is homogeneous. The auxiliary 
field a can be written as the sum of its complex vacuum expectation value in? and a 
quantum fluctuating field a 

a = in? ^aj^ . (4.13) 

Putting 7 = and expanding the effective action around the vacuum expectation value 
of a gives 

N . o Nm^ i^/N ( 1 

5eff - :rTrlog(-52 + m^) - —^pv - -^Tr — ^ ^a 

2 2gl 2 \-d^ + m^ 



+ 



jTrf J 2 ^] +0{ll^). (4.14) 
4 \ -d^ + m^ I 



The last equation shows that a n-point vertices are suppressed by powers of 1 / VA^. 
This is of course due to the definition of a in Eq. (4.13). This definition is convenient 
since now one immediately sees from the vertices that higher loop diagrams are sup- 
pressed. One could of course choose another definition, for example, a = im^ + a. 
Then the same conclusions would follow, since in that case the a propagator has an 
additional suppressing factor of IjN. Inserting a plane wave basis inEq. (4.14) results 
in 



= —X log (P + m ) — + yN X terms hnear m a 

PV 2jp ^ I 2gl 



+ jX &{P)Tl{P, m)a{Py + 0{\/ Va^) , (4.15) 
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where the terms Hnear in a were not written down since these terms (tadpoles) do not 
contribute to the effective potential as explained in Chapter 3. The inverse a propagator 
is equal to 

n{P,m) = -i — -. (4.16) 

2 JQ Q^ +m^ (P + QY + m2 

As discussed in more detail in Chapter 3, the effective potential 'V{np-) can be obtained 
by integrating over the quantum fluctuations a. The leading term of this effective 
potential proportional to N stems from the classical action. The next-to-leading order 
term is originating from the Gaussian integral over the quantum fluctuations a. The 
remaining contributions are of order 1/A^, since one needs at least two three-point 
vertices (which are proportional to 1 / Va^) to form a closed diagram. Choosing the 
sign of the effective potential in such a way that the pressure is equal to the minimum 
of the effective potential one obtains 



2§; 



2 -|/log(p2 + m2)-l/logn(P,m) + 0(l/A^). (4.17) 



The ^ field propagator can be found by differentiating Eq. (4.11) twice with re- 
spect to the sources / and expanding a around its vacuum expectation value. If one 
then keeps terms to next-to-leading order in 1 IN and inserts the wave function renor- 
malization factor Z^ one finds that the scalar propagator in momentum space equals 



D^{P,m) = \ , (4.18) 



where the self-energy is equal to 



i:{P m) = ^ _ ^t. . . 7T77T-T • (4-19) 



1 1 

fg {P + QY + m^UiQ,m) ' 



It is important to keep in mind that m^ is the vacuum expectation value of the 
a field. This vacuum expectation value can be found by minimizing the effective 
potential with respect to m^. Since the effective potential is ultraviolet divergent and 
the renormalization of the nonlinear sigma model is non-trivial, the zero temperature 
case will be discussed separately in the following subsection. Thereafter, the finite 
temperature effective potential will be investigated. 

The effective potential at zero temperature 

At zero temperature one can obtain analytic results for the effective potential. The 
sum over Matsubara modes changes into an integration over momentum. The effective 
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potential 'V to next-to-leading order in the 1 IN expansion is at zero temperature given 



<y = 



It follows after integration over momenta that Ii{P, m), which is the inverse a propa- 
gator, at zero temperature equals (see for example Novikov et al. (1984)) 



1 { VP2 + 4m2 + VP2 ) 

n{P, m) = _ . log — 

47rVp2Vp2+4m2 [^fp^T4^- Vp^ 



(4.21) 



The integrals integrals in Eq. (4.20) will be evaluated using an ultraviolet momentum 
cutoff A. For that one needs the following integral 

r <fp ( Vp2 + 4m2 + Vp2 ) m^ r^ I 1\ 

-log log ^ ^=\^— df 1-- loglog(?) 

- ^ ([rloglog(0 - U(0]|[ + [?loglog(l/0 - U(0]|;^'') , (4.22) 

where 

VA2 -I- 4m2 + Va2 
X - ^ — , (4.23) 

VA2 + 4m2 - Va2 

and li(x) is the logarithmic integral 

\i(x) = P f dt-^. (4.24) 

Jo log(0 

Here P indicates a principal-value prescription for the integral. Since the integrand of 
the logarithmic integral has a pole at f = 1 , a prescription to evaluate the integral is 
necessary. This prescription is the source of the renormalon ambiguity discussed in 
Sec. 4.7. For example the logarithmic integral in the ±ie prescription 

hi(x)= fd? ] . , (4.25) 

Jo log(0 ± le 

differs from the logarithmic integral using the principal value prescription, that is 
h(x) - li±(x) ± in. 

Using that lim;t^i[li(x) - xloglog(x)] = je, where y^ is the Euler-Mascheroni 
constant, and after dropping divergences that are m-independent in the limit A2 ^ oo, 
one obtains 



^, Nm^ Nm^l^ , A2\ 1 



(A2.2.2)ioglog(^)-.2H(^) 



,2 



-I- 2m y^ - 1 - lo^ 



A^ 
4m2 



(4.26) 
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As one can see from Eq. (4.26), the effective potential is highly ultraviolet divergent. 
Its divergences can be removed by renormalization. One way to do this is by adding 
or subtracting infinite, m-independent constants. This merely shifts the whole effec- 
tive potential, but since relative energy differences will not change this does not have 
physical effects. This subtraction is already performed in Eq. (4.26) but does not 
remove all divergences. Another possibility to remove divergences is to absorb m- 
independent infinite constants in a redefinition of the coupling constant gh. However, 
this still leaves two divergences in Eq. (4.26). The only possibility to renormalize the 
term A^loglog(A^/m^) would be to subtract an m-dependent infinite constant since 
limA^oo A^loglog(A^/m^) - A^loglog(A^//i^) is still divergent if yu^ i^ m^. The 
m\i{h} Irn^) term has the same problem. This divergence should be absorbed into 
the coupling constant since it is proportional to m^. But there does not exist a single 
m-independent function which can be subtracted from li(A^/m^) to obtain a finite re- 
sult in the limit A ^ oo. So it seems that the effective potential is not renormalizable, 
unless m-dependent and potentially T-dependent counter terms are introduced. 

Moreover, subtracting the A^loglog(A^/m^) -m^li(A^/m^) term as is done for 
instance by Biscari et al. (1990) for T -0, gives rise to another problem. This term is 
called the perturbative tail and does not arise in dimensional regularization suggesting 
that this subtraction is the correct thing to do. But if one would calculate the effec- 
tive potential using dimensional regularization one would find a ±in ambiguity instead 
(David, 1984) depending on whether the limit 6 J, or e t is taken. This ambi- 
guity is the same as the renormalon ambiguity mentioned above due to the choice of 
prescription in the logarithmic integral. In Sec. 4.5 it is shown that the perturbative 
tail and its accompanying ambiguity cancel when calculating the pressure which is the 
minimum of the effective potential at finite temperature minus the minimum of the 
effective potential at zero temperature. 

To obtain physical quantities like the mass of the scalar fields and the pressure, one 
has to evaluate the effective potential at its minimum. The condition for the minimum 
is given by the equation 

(4.27) 



drrfi 



0. 



Equation (4.27) is often referred to as a gap equation, since solving this equation de- 
termines in leading order the gap in the excitation spectrum at zero spatial momentum 
or equivalently the mass of the scalar particles. Differentiating the effective potential 
with respect to m^, Eq. (4.26), one obtains the gap equation which reads 
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(4.28) 



This gap equation still contains the problematic divergence proportional to li(A^/m^). 
But unlike the effective potential, the gap equation can be renormalized. For that the 
gap equation has to be expressed in terms of the finite physical mass m^. 
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The physical mass m^ of the scalar fields follows from the pole of the (p propagator, 
Eq. (4.18), at zero spatial momentum in Minkowski space. To calculate this propaga- 
tor to next-to-leading order the self-energy, Eq. (4.19), has to be evaluated. At zero 
temperature this self-energy can be simplified by integrating over angles, resulting in 



!(/', m) 



_ 1 r'^ 1 1 

"2^ Jo [{P^ + e^ + m2)2 - 4^232] 1/2 n{Q,m) 



(4.29) 



Since the expression for the propagator, Eq. (4.18), is in Euclidean space-time, one 
has to analytically continue the propagator to Minkowski space to obtain the physical 
mass. Because the self-energy only depends on P^ as one can see from Eq. (4.29), 
continuation to Minkowski space boils down to replacing the Euclidean two- vector P^ 
by the Minkowskian two- vector -p^. Neglecting terms of order l/A'^2^ the physical 
mass m^ yields (Flyvbjerg, 1990), 



ml = m^ - T^^iP, m)\p2. 



nl=m + — ll — 



(4.30) 



Solving Eq. (4.30) for m-^ results in 
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(4.31) 



Since the physical mass m^ should be finite, it follows from Eq. (4.31) that m-^, which 
is the vacuum expectation value of the a field, receives divergent contributions due to 
the quantum fluctuations at next-to-leading order in 1 /N. 

Expressing the gap equation (4.28) in terms of the finite physical mass m^, using 
Eq. (4.31) gives 
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(4.32) 



Because UmA^oo[loglog(A2/m2) - loglog(A2/ju2)] = q, the gap equation can now be 
renormalized. Making the substitution g^ —^ ^g^S^ *i'^h g^ = g^ijj) and 
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1 + -:- 1 - T7 log— 7 + T7::^loglog— r , 



one obtains the renormalized gap equation which reads 



^---1-T7 log — + -(?£ + log 4) . 
g^ip) \ Nj mi N 



(4.33) 
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The expression Eq. (4.33) for Z 2 is exact in g^iji) up to order l/N'^ corrections and 
results in the known next-to-leading order /3-functions as calculated by Rim and Weis- 
berger (1984), Orloff and Brout (1986) and Biscari et al. (1990) 

Pisl) = A-7^--|l-^l^|l + 77^h (4-35) 

The negative sign in front of the y6-f unctions shows that the nonlinear sigma model is 
indeed asymptotically free. The next-to-leading order beta function gives already the 
correct result for A/^ = 2 because for a free theory the beta function vanishes. This is an 
indication that the 1 /N expansion in the nonlinear sigma model might work well down 
to low values of A'^. Solving the renormalized gap equation (4.34) for the physical mass 
m^ gives 




ml - M^ exp < - 



An 2 

N 



— --(T£ + log4) 



/|l-|ir (4-37) 



This equation illustrates the need for the non-perturbative approach, since the physical 
mass turns out to be non-analytic in the coupling constant. For N = 2 this gives a zero 
physical mass, reflecting the fact that the 0(2) nonlinear sigma model is a free field 
theory. 

Using the gap equation, one can obtain the value of the effective potential 'V{m^) at 
the minimum which will be denoted by 'V^^^^{m^^, where mo = m^ given by Eq. (4.37). 
In terms of bare quantities, one finds 

2 

T^=0(m2) ^ - (A^ - 2) ^ - ^a2 log % . (4.38) 

87r 87r g2 

In contrast to the eff'ective potential as a function of m^, this equation shows that the 
efl^ective potential at the minimum can be renormalized. It is expressed in terms of 
a finite mass which can be found by solving the renormalized gap equations. The 
quadratic divergence that is left in the minimum of the effective potential can be re- 
moved by a subtraction. However, as explained before this term is ambiguous, it 
would differ by a constant in if a different prescription is used. The expression for 
the minimum of the effective potential at zero temperature will be subtracted from the 
minimum of the effective potential at finite temperature in order to obtain a finite pres- 
sure which is moreover unambiguous. As is shown in Sec. 4.5 this ambiguity can be 
removed without renormalizing the effective potential as well. 

For completeness, the wave-function renormalization factor of the scalar fields will 
be derived in this paragraph. The wave-function renormalization constant Z^ follows 
from requiring that the pole of the propagator has a residue equal to unity. This implies 
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that 



1 dI(P,m) 
Zih — \ — —- 



N dP^ 



= l + iloglog(^| . (4.39) 



The wave-function renormahzation Eq. (4.39) is in accordance with that obtained 
by Flyvbjerg (1990). Rim and Weisberger (1984) have calculated the wave-function 
renormahzation constant in dimensional regularization. Their result also agrees with 
Eq. (4.39) upon identifying log(A^///^) -^ 2/e where d = 2 - e. 

The effective potential at finite temperature 

The results at zero temperature are obtained analytically, but at finite temperature this 
is in general not possible. Therefore, numerical methods will be applied to calculate 
thermodynamical quantities. However, it is still possible to isolate the ultraviolet di- 
vergences analytically and discuss the renormahzation at finite temperature without 
the need for numerical evaluation. 

The eff"ective potential through next-to-leading order in l/N can be written as 
n^(m2) = N'VLoim^) + ^nlo('«^) where 



^LO - ^ - \i log {P^ + n?) , ^NLO = -^i log n{P, m) . (4.40) 



After summing over Matsubara modes and subtracting divergent m- and T-independent 
terms one obtains in the limit A ^ oo for the leading order effective potential 



m^ m^ I , , A^\ 1 _2 



^LO = -T - ^ 1 + log ^ + —T'JoiJim) , (4.41) 

where 

hifim) - 4 r ^^«K) ■ (4-42) 
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Here cjOq - ^[cf^+n? and n{x) - [exp(/3x) - 1]~' is the Bose-Einstein distribution 
function. 

To calculate the next-to-leading order contribution to the eff"ective potential one 
needs to know the inverse a propagator Il{P, m) at finite temperature which is de- 
fined in Eq. (4.16). Summing over Matsubara frequencies and averaging over angles, 
IV{P, m) reduces to 
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Figure 4.1: The function F\/T^ as a function of ySm. This function contributes to the 
next-to-leading order finite temperature effective potential. 



This expression cannot be simplified further. It shows that the inverse a propaga- 
tor Il{P, m) is temperature-dependent and ultraviolet finite. The next-to-leading order 
contribution to the effective potential is not finite though. 

The next-to-leading order contribution to the effective potential can be written as 



'y^i.o = -^{Di+D2 + F^+F2) , 



(4.45) 



where D\ and D2 contain divergences and F\ and F2 are finite terms. As explained in 
Sec. 2.5 and as verified numerically, the quantity 



^£ 



\ogIV{P,m), 



(4.46) 



is finite since its integrand is suppressed by a Bose-Einstein like distribution factor. 
To calculate F\ a modified Abel-Plana formula (see Sec. 2.5) was used. The function 
Fi/T-^ depends on/Jm only and is displayed in Fig. 4.1. As is discussed in more detail 
in Sec. 4.6 limm^o F\{m) = 7t/3T^ which can also be seen from Fig. 4.1. The limit of 
/3m ^ 00 is equivalent to fixing m and taking T ^ 0. In this limit Fi should vanish 
because it is the difference of a finite temperature and zero temperature contribution. 

Because Fi is finite all divergences of the next-to-leading order effective potential 
arise from 

;^logn(P,m). (4.47) 
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The nonlinear sigma model only contains ultraviolet divergences. It is infrared finite 
since the mass m is non-vanishing. Therefore, in order to isolate the divergences in 
Eq. (4.47) one has to consider the limit of momenta large compared to T . Due to 
the suppression by the Bose-Einstein distribution function of high momentum modes, 
only internal momenta ^'s up to order T contribute to YVt{P, ni). Hence, if the external 
momenta P are much larger than T one can neglect factors of q compared to factors 
of P in the high-momentum limit. In this way one can approximate the inverse a 
propagator as follows 

1 P2 

YlT{P,m) ~ —rJi , (4.48) 
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(4.49) 



In order to split ofi" the prefactor in front of the logarithm of the zero temperature part 
of IV{P, m), it is useful to write n(P, m) as 



n(P, m) = 
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4^ yp2 (p2 + 4^2) 
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(4.50) 



(4.51) 



The factor ^P^{P^ + 4m^) gives the following divergent contribution to the free en- 
ergy 

1 f d^P / , ,x m^\ / A2 

The divergent term Di is temperature independent. Hence it is also contained in the 
eff"ective potential at zero temperature as can be verified from Eq. (4.26). The other 
divergent term D2 and the finite term F2 can be found by calculating 
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log U{P, m) = D2+F2 . 



(4.52) 



In order to isolate the divergent term D2, one needs the large-P behavior of n(P, m) 
which by using Eq. (4.48) is found to be 
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where HM stands for high-momentum approximation. By integrating log IIhmC^. w) 
over momentum it is found that the divergent term D2 is given by 
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here m^ = m^ exp[-/i(/3m)]. In the hmit of zero temperature it can be seen that the 
divergence D2 coincides with the divergences found in the effective potential at zero 
temperature, Eq. (4.26). Since fh^ depends expUcitly on temperature, Eq. (4.54) shows 
that the next-to-leading order effective potential contains temperature dependent diver- 
gences. This seems surprising since one expects that scales at which finite temperature 
are important are much smaller than the ultraviolet cutoff. As a result, ultraviolet diver- 
gences should be independent of temperature. In Sec. 3.4 a general argument is given 
why these T-dependent divergences arise. Hence apart from problems at zero temper- 
ature discussed in the previous subsection, it even becomes impossible to renormalize 
the effective potential in a temperature independent way. But, as will be discussed in 
the next section, these problems disappear at the minimum of the effective potential. 

The finite function F2 defined in Eq. (4.52) is the difference of two divergent quan- 
tities. Due to the finite machine precision the maximal reachable accuracy of the dif- 
ference would be rather small if the divergent quantities are calculated separately. It is 
therefore better to write D2 partly in terms of an integral in the following way 



,2 -fc ^"'^ \ ^ + ^ log log ^ 
r)^ \m^/ In \m^ 



£»2 = ^J^loglog| — | + — loglogl^l . (4.55) 



The last equation can be used to write F2 in terms of an integral and a very slowly 
diverging term. 
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(4.56) 



To calculate F2 numerically, the integral in Eq. (4.56) was evaluated up to a cutoff 
A. Then for increasing values of A, F2 was calculated until convergence to a finite 
number was found. This convergence indicates that by the method discussed above 
indeed all divergences were isolated. In the limit of zero temperature it can be seen 
from Eq. (4.26) that F2 = m^jElln. It is convenient to subtract this zero temperature 
contribution from F2 in order to define F2 as 

2 
F2 = F2- ^7E . (4.57) 

The function F2/T'^ is a function of Pm and is displayed in Fig. 4.2. In Sec. 4.6 it is 
shown that for small ySm, F2 — > 0. Consequently, F2 ^ in agreement with Fig. 4.2. 
For large values of /3m, which implies low temperatures, F2 -^ 0. This is because in 
the definition of F2 the zero temperature contribution is subtracted. 

Putting everything together, the finite temperature effective potential becomes 
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Figure 4.2: The function F2lT'^ as a function of ySm. This function contributes to the 
next-to-leading order finite temperature efi'ective potential. 

The gap equation (4.27) at nonzero temperature is 
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(4.59) 



From the fact that ^^ is temperature independent, one can conclude that m^ is also tem- 
perature independent at leading order in the 1 I'M expansion, when it is a solution to 
the gap equation. This will be used later on to conclude that the pressure can be renor- 
malized in a temperature-independent way. The renormalization of the gap-equation at 
finite temperature proceeds analogously to the zero-temperature case discussed in the 
previous subsection. To remove the problematic li divergence the gap equation will, 
like was done in the zero temperature case, be expressed in terms of the physical mass 
m^ at finite temperature. This mass can be obtained by finding the pole of the finite 
temperature (p field propagator, Eq. (4.18), analytically continued to Minkowski space. 
Here the finite temperature physical mass will be defined as the pole of the propagator 
at po = and p^ = -m^,. It is found that 
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Figure 4.3: The function logifim^) as a function of T, for different values of N. 



where m^ = m^exp[-J\(Jim^)] and Ft, is a finite function that depends on the temper- 
ature as well as m^. Since m^ is merely used as a way to express the renormalized gap 
equation in terms of finite quantities, any choice of Ft, will do and we choose Ft, = 0. 
It was checked numerically that other choices indeed do not alter the final result for the 
pressure. Using Eq. (4.60), the gap equation expressed in the physical mass becomes 
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(4.61) 



The gap equation can be rendered finite by the substitution g^ -^ Zgig^Qi), where Z'2 
is the zero temperature renormalization constant given by Eq. (4.33). Hence unlike 
the effective potential, the gap equation can be renormalized T -independently without 
problems. The renormalized gap equation becomes 
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(4.62) 



In Fig. 4.3 the solution to this renormalized gap equation for the arbitrary choice 
g^ifi = 500) - 10 is displayed as a function of temperature and for different values of 



68 



Thermodynamics of the nonlinear sigma model in d=l+l 



V 




Figure 4.4: Pressure of the 0{N) nonlinear sigma model to next-to-leading order in 
I IN, normalized to NT^, as a function of temperature, for different values of A'^. 

4.5 Pressure 



By solving the renormalized gap equation (4.62) one can find the minimum value of the 
effective potential which is equal to the pressure. The physical mass at this minimum 
will be denoted as m? - m^AT) and m^(0) - mi{T = 0). By expressing the effective 
potential in terms of the physical mass using Eq. (4.60) the pressure can be calculated. 
For this one has to expand the /q and J\ functions in 1 /N. Furthermore, one should 
also take into account the zero temperature results since the pressure is the difference 
of the minimum of the effective potential at finite temperature and at zero temperature 
(given in Eq. (4.38)). As a result due to this difference the problematic divergences 
drop out and one obtains 
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(4.63) 



This finite pressure was evaluated numerically, after solving Eq. (4.62) for m^(T). 
The result for different values of A'^ is shown in Fig. 4.4, for the arbitrary choice 
g^iM = 500) = 10, hence T is given in the same units as /i. As can be seen in the 
figure, P/NT^ approaches an A/^-dependent constant (to be evaluated in the next sec- 
tion) at large temperatures. This constant is exactly the pressure of a non-interacting 
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Figure 4.5: Pressure of the 0(A'^) nonlinear sigma model normalized its value at T 
oo, as a function of temperature for different values of N. 



gas of A'^- 1 bosons divided by N, which is "PsqIN - {N-l)nT^/6N. This is an imme- 
diate consequence of the asymptotic freedom and the fact that one degree of freedom 
is effectively removed due to the constraint. This expression coincides with the per- 
turbative result, Eq. (4.7), as expected. In the limit of zero temperature the pressure 
divided by T^ approaches zero exponentially as it effectively becomes a gas of heavy 
quasi-particles, see Fig. 4.3. The figure shows that corrections are really of order l/N, 
which indicates that the 1 /A'^ is a good expansion for the nonlinear sigma model. 

If the pressure, for a given value of N, is normalized to its value at T = oo, the 
normalized pressure has a very small dependence on A'^ as is displayed in Fig. 4.5. 
Similar behavior is found in lattice calculations of the normalized pressure of the pure 
SU(A^c) gauge theory in 4-dimensions (Bringoltz and Teper, 2005), see Fig. 4.6. Also 
lattice calculations performed by Karsch et al. (2000) show that for full QCD the 
normalized pressure turns out to be rather insensitive to the number of flavors. 

To obtain the pressure, the next-to-leading order gap equation was used. But as 
was argued by Root (1974) one in principle only needs the leading order gap equation 
in order to obtain the value of the effective potential at the minimum. Writing the 
solution to the gap equations as 
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(4.64) 
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Figure 4.6: Results of lattice calculations by Bringoltz and Teper (2005) of the pressure 
of pure S\]{N) Yang-Mills theory normalized to the pressure of a free gas iox N - 3 
(triangles), N - A (circles) and N = S (squares). The solid line is the N = 3 result of 
Boyd et al. (1996). 



and Taylor expanding the effective potential, one obtains (up to 0{1/N) corrections) 
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(4.65) 



The last term of Eq. (4.65) vanishes by using the leading-order gap equation. The 
pressure f can now be written as 



P = NPlo + :Pnlo 
From the discussion above it follows that 
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where m| is the solution of the leading-order gap equation, at temperature T. By 
using the leading order gap equation it can be shown that at the minimum, part of D2 
becomes 
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(4.68) 
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Since D2 contains the logarithmic integral it is dependent on the prescription and has 
an ambiguity. From Eq. (4.68) it follows that the problematic T-dependent divergences 
become T-independent at the minimum of the effective potential, and vanish in the 
calculation of the pressure due to the subtraction of the zero temperature contribution. 
As a result the pressure is unambiguous, i.e. independent of the choice of prescription. 

4.6 High-temperature approximations 

Bochkarev and Kapusta (1996) consider the nonlinear sigma model in 3+1 dimen- 
sions, at next-to-leading order in the 1/A'^ expansion. Since their result for the pressure 
cannot be obtained analytically, they resort to a "high-energy approximation". In this 
section the same approximation for the 1 -1- 1 -dimensional case will be made to compare 
it with the exact numerical results obtained in Sec. 4.5. 

The essence of the high-energy approximation of Bochkarev and Kapusta (1996) 
is that one neglects the zero temperature part of IV{P, m) and approximates it by its 
temperature-dependent part nr(/', m). According to Bochkarev and Kapusta (1996) 
one can then approximate n7-(P, m) by its high-energy behavior. In the present case of 
the nonlinear sigma model in 1 -1- 1 dimensions, these approximations amount to 

1 P^ 

I{{P, m) « IIt{P m) « —— ^-— -h , (4.69) 

4n (pI + col){pl + ojI) 

where aj+ = ^^fp^^+n? ± m. This expression is identical to the high-momentum ap- 
proximation to IVT{P,m), Eq. (4.48). In the high-energy approximation the effective 
potential is given by 

Nm^ NC^ / 2 2\ IX. 2 
^,,,.____|^log(p^..^)--|^logP^ 

+ \^ log (pI + ojI) + 1^ log (pI + oj^_) + 0{\/N) . (4.70) 
The resulting expression for the gap equation is 
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Again, the gap equation requires coupling constant renormalization. In this approxi- 
mation, it is found that the renormalization constant is 
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which differs slightly from the exact renormalization constant to next-to-leading order 
in 1/A^, Eq. (4.33). It is in fact equal to the perturbative one as for instance is discussed 
in Peskin and Schroeder (1995). Making the substitution g^ -^ '^g^S^^ ^"^^ obtains 
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(4.73) 



where the functions Kf are defined as 
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By substituting the gap equation (4.71) into the expression for the effective potential 
in the high-energy approximation, Eq. (4.70), the pressure becomes 



^HEA - ^ [joT^ + (Ji - l)m^] + \t^ 



^ [{K^ + K-)T^ + {K^ +K-- T)m'\ , (4.75) 



where the functions K^ are given by 

8 rdpp\ 

From Fig. 4.7, one can see that the high-energy approximation underestimates 
the pressure compared to the exact result. The advantage of an approximation like 
the high-energy approximation is that the analytic calculations are simpler and that it 
is easier to implement numerically. However, it is not a satisfactory approximation 
because one cannot really justify that the temperature-independent part oiIl{P,m) can 
be left out. It may even be considered surprising that this approximation works. In 
Chapter 6 it will be shown that the high-energy approximation in the nonlinear sigma 
model in 3 -I- 1 dimensions actually leads to wrong results. 

It is possible to derive a different approximation, which is better than the high- 
energy approximation. The inverse a-propagator IT at finite temperature can be ob- 
tained by first integrating over the momentum which gives 



qo=2nnT Jm'^ + ^^ 



U\P,m) = — > — ^^ ^-^ . (4.77) 



Since it follows from the leading order gap equation that for high temperatures and for 
all values of the coupling constant, m «; T see Fig. 4.3, Il{P, m) can at high tempera- 
tures (HT) be approximated by keeping only the ^o = mode in the sum Eq. (4.77) as 
follows 

n(P, m) « UniiP, m) = \^ ^^ ^ , i i ■ (4-78) 
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Figure 4.7: Exact pressure as a function of temperature at next-to-leading order for 
N = S compared with the high-energy and high-temperature approximations. 



By using that P^ + Am^p^ can be written as [/j^ + {p + im)^ + m'^][p^ + (p - ifn)^ + m^] 
and shifting p ^ p ±im after taking the logarithm, it is found that 



i \og{P'^ + Am^p^) = iL \og{P^ + m^) . 



(4.79) 



Hence the function F\ can be approximated by 



i' 



Fi « 4 logUm{P,m) - l-T^J^,(J3m) - ^T^ ^ ^T^ 



(4.80) 



The function F2 « in the high-temperature approximation because j log Uht does 
not contain a finite term that depends on temperature. The part proportional to logjSw 
gives rise to a divergence which by definition is not part of F2. These high temperature 
limits agree with the numerical calculations displayed in Figs. 4.1 and 4.2. 

The numerical calculation of Fi + F2 shows that for m/T < 0.1 this approximation 
has an error smaller than 10 percent. To obtain the pressure in the high-temperature 
approximation, Eq. (4.63) was applied using the approximations to Fi and F2. The 
result for the pressure in the high-temperature approximation is shown for comparison 
in Fig. 4.7 (again for g'^iji = 500) = 10). It can be seen that the high-temperature 
approximation is much better than the high-energy approximation. 

One can approximate the pressure even further by expanding the functions Jq and 
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J I in the limit fim -^ 0: 

/o - ^ - 47r/3m - 2 flog ^ + y^ - i j (fimf + O (oSm)^) , (4.81) 

7i = -^ + 2 (log ^ + ys) + O (0Sm)2) (4 32) 

Inserting the approximations given in Eqs. (4.80), (4.81) and (4.82) into the gap equa- 
tion (4.62), one obtains 



ySm w n 



In log4\/ 2\ n^r' 



(4.83) 



which indicates that fim ~ l/logT for large T. In the limit m/T -^ 0, the high- 
temperature approximation of the pressure is 

P n I 1\ / 2\ m 

7 1-77 - 1-77 b;;' (4-84) 



NT^ 6\ Nj \ NjAT 

where the first term is the pressure of a gas of free massless particles with A'^- 1 degrees 
of freedom. 

4.7 Thermal infrared renormalons 

An observable O in a quantum field theory like QCD is often only known in terms of 
an expansion in the coupling constant a. Because negative a gives rise to unphysical 
behavior as was argued by Dyson (1952) for QED, it is expected that this series expan- 
sion has zero radius of convergence. Hence one should view it at best as an asymptotic 
series. One might wonder whether it is possible by resumming the series expansion to 
learn about the non-perturbative physics. In specific cases this is possible, but more 
often one finds so-called renormalons, which give rise to ambiguities in the relation 
between O and its resummed version. These ambiguities may give an insight in the 
size of the non-perturbative corrections. 

Since one can apply the 1 IN expansion to the nonlinear sigma model, observables 
can be calculated exactly in terms of the coupling constant. These exact expressions 
can be used to study the renormalons in the nonlinear sigma model as was done at zero 
temperature by David (1982, 1984, 1986) and Beneke et al. (1998). In this section this 
analysis is extended to finite temperature. 

This section consists of two parts. In the first part a short introduction into renor- 
malons will be given. In the second part the nonlinear sigma model will be used to 
investigate the temperature dependence of the renormalons and how they affect the 
possible definitions of the renormalized effective potential, pressure and gap equation. 
For a more extensive review on renormalons see Altarelli (1995) and Beneke (1999). 
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Renormalons 

An observable O can in general be expressed in terms of a perturbation series ina - g^ 
in the following way 

0{a) ^Y^ana'' . (4.85) 

It often happens that this series is an asymptotic series which is not convergent. Such 
an asymptotic series can sometimes be resummed by the Borel transform method, 
although in general this resummation is not unique. The Borel transform of 0{a) is 
defined as 



5(0 =y,-/ ■ (4.86) 

The observable O can now be expressed in terms of an integral over its Borel transform 



1 C^ 

0(a) = - dte-""B(t), (4.87) 

« Jo 

which has the series expansion of Eq. (4.85). For example if one takes a„ = {-l)"n\, 
the series expansion for O is clearly divergent. The Borel transform however is con- 
vergent for \t\ < 1 and equal to B{t) = 1/(1 + t). To resum 0(a) it is assumed that 
B(t) = 1/(1 + for any value of t. The resummed expression for O now yields 



1 r 

{a) = - 

a Jo 



0(a) = - I dte~""— . (4.88) 

\ + t 



The integral in Eq. (4.88) is convergent and, moreover, well-defined since there are no 
singularities along the integration contour. So in this example the asymptotic series 
has been resummed unambiguously. However, sometimes this method gives rise to 
ambiguities. Take for example a„ = nl in which case B{t) = 1/(1 - t) and 

1 r°° 1 

0(a) = - dte-"" . (4.89) 

a Jo l-t 

In this case the integral is not well defined since there is a singularity at f = 1 on the 
integration contour. This singularity is called the renormalon pole. One can still give a 
meaning to this resummed series by introducing a prescription how to integrate around 
the pole 

1 r°° 1 

0(a)± = - dte~""- . (4.90) 



\ -t±i€ 



This resummation gives rise to ambiguities, since one does not know whether to take 
the -I-, - or a combination like the principal value prescription. But as long as the 
differences between the prescriptions are not that big, it does not really matter. The 
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difference is equal to the residue nit - 1, which is called the renormalon residue. In 
this case it is equal to 

0(a)+ - 0(a)_ = — exp (-l/or) . (4.91) 

a 

So in this example the series can be resummed up to ambiguities of order exp(- 1 1 a) I a. 
These ambiguities are tiny when a is small. However for larger values of a they can 
become important. These ambiguities indicate that the non-perturbative corrections to 
Eq. (4.85) could be in this case of order exp(-l /q')/q'. 

Renormalons in the nonlinear sigma model 

Renormalons appear in the vacuum expectation value of the a field. A calculation 
shows that at zero temperature one can write (David, 1986) 



Anm^ r d^p \ dn ^/ 1 \ 



where m^^ = A^ exp(-47r/g^). This equation is in agreement with the gap equa- 
tion (4.28) if m^ is split as follows m^ - m^^ + m^^^jN. The part of the integral 
in Eq. (4.92) that has the infrared renormalon pole in the Borel plane is in fact the 
contribution from the integrand in the limit P-^/m^ -^ oo (its "perturbative tail") 



/ 



(fP 1 dUoo 1 /A2\ 1 A2 ^ 



{2n)^ Iloo dm^ An \m^j An m^ 



where x - log(A^/m^) and Ei(x) the exponential integral which is defined below. In 
the limit x ^ oo, the logarithmic integral has the asymptotic expansion 



pQ I /~*oo —hx 

'Ei(x) - y -% + ine''' = \ Ab ^—- + ine'"^ , (4.94) 

f-^ x"+i Jo 1 - ^ 



where aig{b) = ±s. From Eq. (4.94), it is clear that there is a renormalon pole at 
b = \. This shows that when A ^ oo the value of {a) is inherently ambiguous at 
next-to-leading order, due to the freedom in the choice of prescription. David (1984) 
has shown that this ambiguity also arises in dimensional regularization. 

The same problem appears in the calculation of the effective potential even at its 
minimum, but not in the gap equation. The latter can be seen from the last term of Eq. 
(4.20) which contributes to the gap equation as follows 

\ d r <fP . \ r <fP I dU 



2dm^ 



r d^P If d^P 1 dU 
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The ambiguity that would arise from this term when removing its perturbative tail (of. 
Eq. (4.93)) cancels in the gap equation (4.28) against the one arising in m^ (cf. Eq. 
(4.92)). 

The perturbative tail of the effective potential, i.e. the first two terms of D2 de- 
fined in Eq. (4.54), corresponds to poles in the Borel plane at b = and b = I, 
respectively. Since m^ is only temperature independent at the minimum (at LO only, 
but that is sufficient because one is working at NLO), the subtraction of the perturba- 
tive tail will become temperature dependent, except at the minimum. Since subtract- 
ing temperature-dependent divergences renders the remaining temperature-dependent 
terms meaningless, it is impossible to define a finite effective potential at finite temper- 
ature. In order to avoid any renormalon ambiguity, it also not possible to obtain a finite 
effective potential or even an unambiguous minimum at zero temperature. However, 
as was shown in Sec. 4.5 the quantity f^ - ^^=" is free from renormalon ambiguities 
and is finite after temperature-independent coupling constant renormalization. 

Finally, using Eq. (4.92) one can investigate the finite temperature dependence 
of renormalon contributions to {a) and the effective potential. One can show that 
Eq. (4.92) at finite temperature has exactly the same renormalon contribution, i.e. 
neither the pole nor the residue become temperature dependent. Secondly, the per- 
turbative tail of the effective potential which is given by the first two terms of D2, 
corresponds to poles in the Borel plane at Zj = and b - \. The positions of the renor- 
malon poles are not affected by temperature. Only the residues become temperature 
dependent, except at the minimum of the potential, as we concluded earlier. The fact 
that renormalon pole positions are not affected by temperature, but residues are, is also 
the case for the thermal ultraviolet renormalons in 0"* in 3 -1- 1 dimensions studied by 
Loewe and Valenzuela (2000). 

4.8 Summary and Conclusions 

To summarize, the pressure in the nonlinear sigma model was calculated at finite 
temperature to next-to-leading order in the 1/A'^ expansion. The main result is that 
one can obtain an unambiguous, finite pressure, by subtracting the zero-temperature 
value of the pressure and renormalization of the coupling constant in a temperature- 
independent way. This procedure cannot be carried out away from the minimum of 
the effective potential and it was argued that defining a finite, effective potential by 
the subtraction of the so-called perturbative tail, leads to ambiguities associated with 
infrared renormalons. In general, these become temperature dependent, and this casts 
doubt on the usefulness of defining a finite effective potential outside the minimum. 
Since it turns out that, as it should, physical quantities are independent of the choice 
of prescription, one can apply any prescription for the logarithmic integral without 
worrying about the possible ambiguities. This is what will be done in the remaining 
chapters of this thesis. 
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The expression for the pressure was calculated numerically at finite temperature. 
This calculation shows that the 1 IN expansion is a meaningful expansion for all tem- 
peratures. The high-energy approximation that was originally applied to the nonlinear 
sigma model in 3-1-1 dimensions by Bochkarev and Kapusta (1996) was also studied. 
In 1-1-1 dimensions, where one can compare with exact numerical results, it was found 
that it underestimates the pressure for all temperatures. An improved approximation 
was suggested, the so-called high-temperature approximation. This approximation has 
the advantage that it is quite easy to produce numerical results and agrees better with 
the exact results. At asymptotically high temperatures the pressure approaches that of 
a gas of A^ - 1 free massless particles. It was found that the pressure divided by the 
pressure at infinite temperatures has a very weak dependence on N. Similarly behavior 
is found in lattice calculations as a function ofNf and Nc- 



Chapter 5 

Thermodynamics of the CP^~^ model 



Using the 1 /N expansion, the influence of quantum instantons to the thermodynamics 
of the CP^^"^ model in 1+1 dimensions is studied. The auxiliary field effective po- 
tential is calculated to next-to-leading order in 1 /N and turns out to have temperature- 
dependent ultraviolet divergences. These divergences can only be renormalized at its 
minimum just like in the nonlinear sigma model discussed in Chapter 4. By using that 
the CP^ model is equivalent to the 0(3) nonlinear sigma model it is argued that the 
pressure for intermediate temperatures is dominated by the effect of quantum instan- 
tons. This chapter is based on: The effect of quantum instantons on the thermodynam- 
ics of the CP^-^ model, J.O. Andersen, D. Boer and H.J. Warringa, hep-th/0602082. 

5.1 Introduction 

It was discovered by Belavin et al. (1975) that the classical equations of motion of 
Euclidean QCD have finite action solutions with nontrivial topology. These solutions 
are called instantons. They are stationary points of the classical action. In perturbation 
theory one most often only takes into account the trivial vacuum and small fluctuations 
around it. However, in the non-perturbative regime the instanton solutions and the 
fluctuations around it can really contribute to physical quantities as was first observed 
by 't Hooft (1976). He showed that instantons are giving rise to an additional source of 
U(1)a symmetry breaking in QCD, which is necessary to explain the relatively large 
mass of the ?/ meson. 

In this chapter it is investigated whether instantons can have an influence on the 
thermodynamical quantities. To answer this question, the CP^~^ model in l-i-l dimen- 
sions will be studied since it admits instanton solutions for all A^. The CP^~^ model 
was introduced by Eichenherr (1978). It is possible to investigate its non-perturbative 
regime using the l/N expansion (D'Adda et al, 1978). An attractive feature of the 
QpN-i model in l-i-l dimensions are its similarities with QCD in 3-1-1 dimensions. 
Like QCD it is an asymptotically free theory. Moreover, classically the CP^^^ model 
is scale-invariant. Due to renormalization of the quantum corrections a mass scale is 



80 Thermodynamics of the CP^ ' model 

introduced in the model. Furthermore, a mass gap for the scalar fields is generated 
dynamically by the interactions. The local U(l) symmetry of the CP^^^ model gener- 
ates a long range interaction between the scalar fields. Since the Coulomb potential in 
1 + 1 dimensions is proportional to the distance between the charges, this model has 
confinement (D'Adda et at, 1978). Finally as alluded to above, the classical equations 
of motion of the CP^'^ model admit instanton solutions (Golo and Perelomov, 1978). 
Lazarides (1979) investigated the CP^^"^ model at finite temperature for the first time 
and showed that at non-zero temperature the model is no longer confining. For a re- 
view on the CP^^^ model at finite temperature, see Actor (1985). Instantons in QCD 
at finite temperature are discussed by Gross et al. (1981) and in the large Nc limit by 
Schafer (2004). 

As is explained in Chapter 3, thermodynamical quantities like the pressure can be 
obtained by minimizing the finite temperature effective potential. In this chapter the 
finite temperature effective potential will be calculated to next-to-leading order in 1 jN. 
Its zero-temperature counterpart was calculated by Campostrini and Rossi (1992). It 
turns out that the effective potential contains temperature-dependent ultraviolet diver- 
gences. These divergences become temperature-independent at the minimum. Ther- 
modynamic quantities like the pressure are defined at the minimum of the effective 
potential and hence can be calculated without such renormalization problems. The 
same phenomenon is also found in the 0(A'^) nonlinear sigma model in l-i-l dimen- 
sions (Chapter 4) and the 0(A'^) linear sigma model in 3-1-1 dimensions (Chapter 6). A 
general explanation of how these temperature divergences arise is given in Sec. 3.4. 

Instantons are characterized by a quantized winding number Q, which is a topo- 
logical invariant. The instantons give a contribution proportional to e\p{-nNQ/g^) to 
the partition function, where the non-analyticity in gj^ shows that the contribution is 
non-perturbative. This contribution also indicates that instantons effects disappear in 
the limit A'^ ^ oo as was argued by Witten (1979). As a result one has to perform 
calculations to next-to-leading order in l/N to investigate the effects of instantons 
in the 1/A'^ expansion. It was argued by Witten (1979) using the 1/A'^ expansion at 
zero temperature that quantum corrections let the instanton configurations to disap- 
pear. However, Jevicki (1979) showed that classical instantons still are present in the 
quantum effective action, not longer as stationary solutions but as poles. However, at 
finite temperature, Affleck (1980a,b,c) found that the large-A'^ quantum effective action 
does contain stationary solutions with quantized topological charge, called quantum 
instantons. Hence, at finite temperature configurations with non-trivial topology can 
contribute to physical quantities. In this chapter the same conclusion will be drawn. 
The pressure is calculated for topological configurations with zero winding number. 
The contribution from other configurations will be ignored since they are difficult to 
calculate. It is found that for intermediate temperatures, where the leading order con- 
tribution to the pressure divided by T^ displays a sharp increase, the next-to-leading 
1 /N contributions gives rise to a negative pressure. This is unphysical because it re- 
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suits in a negative entropy. Therefore, the negative pressure indicates that an important 
contribution, most probably from the configurations with non-zero winding number, is 
left out. This conclusion is further strengthened by using a correspondence discovered 
by Eichenherr (1978) between the CP^ model and the 0(3) nonlinear sigma model. 
In the 1/A'^ approximation of the 0(3) nonlinear sigma model one implicitly includes 
all instanton configurations. Hence the difference between the pressure of the Q - 
sector of the CP^ model and the pressure of the 0(3) nonlinear sigma model gives 
the contribution of the topological configurations with Q i^ to the pressure of the 
CP^ model. In this way a strong hint is found that the topological configurations with 
Qi^O give a large contribution to the pressure for intermediate temperatures. 

This introduction will be ended by mentioning some related studies of the CP^^^ 
model. Schwab (1982), Miinster (1982) and Miinster (1983) have investigated the 
CP^^^ model on the sphere 5^. Miinster (1983) remarked that it is naive to expand 
only around the saddle point A^ = because it leads to an improper treatment of zero 
modes. Furthermore, the CP^^^ model has been studied on the lattice as well, see 
for example Campostrini et al. (1992a,b) and Vicari (1993). A way to investigate the 
relevance of the topological solutions is by adding a term i}Q to the action. Olejnik 
and Schierholz (1994) and Schierholz (1995) calculated iov N = A the free energy as a 
function of §. They found that depending on the size of the coupling constant there is 
a phase transition from the confined to the deconfined phase for ■& < n. Azcoiti et al. 
(2004) found that forN = 10 the CP symmetry is spontaneously broken for ■& = n. 
The CP^^^ model also has applications in condensed matter physics. Pruisken et al. 
(2003) and Pruisken and Burmistrov (2005) have applied this model to investigate the 
quantum Hall effect. Ichinose and Yamamoto (1990) studied antiferromagnetism using 
the CP^"^ model at finite temperature. 

This chapter is organized as follows. In Sec. 5.2 the CP'^^^ model is introduced. 
The calculation of the effective action and effective potential are performed in Sees. 
5.3 and 5.4. In Sec. 5.5 the results of the calculation of the pressure are discussed. 
Finally a summary and conclusions are given in Sec. 5.6. 

5.2 The CP^-^ model 

The CP^^^ model is described by the following Lagrangian density 



X = ^5^<^*5^0,- + Xint , <p]<Pi=Nlgl, i=\...N, (5.1) 



1 

where 0(x) is a complex scalar field and gh is the bare coupling constant. Under 
an U(l) transformation which is parametrized by cr(x), 0(x) -^ exp[/o"(x)]0(x). By 
requiring that the Lagrangian density is invariant under local U(l) transformations, the 
interaction term Xim can be determined from the transformation rule of the Lagrangian 
density which reads 

6£ - i{d^<T) (5^0*) cPi + ^Xint = . (5.2) 
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Hence the interaction term should contain derivatives. If the interaction Lagrangian 
density is chosen as follows 

Xint = ^(0;5^0,) (0*5^0,) , (5.3) 

then Eq. (5.2) is satisfied. As a result the model becomes invariant under global SU(A'^) 
transformations as well. 

The Lagrangian density can also be written in terms of gauge fields A^ which 
transform under a U(l) gauge transformation as 

A^{x) ^ A;(x) = A^{x) - d^,(T{x) . (5.4) 

The Lagrangian density expressed in gauge fields is given by 

Z^\ \dM'- , f.cPi = N/gl . (5.5) 



where the covariant derivative D^ = d^ + iA^. By solving the equations of motions of 
the A^ fields one finds 

A, = /|0;5,^, . (5.6) 

If this expression is inserted into Eq. (5.5) the original Lagrangian density Eq. (5.1) 
and Eq. (5.3) is recovered. 

Coset models 

The CP^^^ and the 0(A'^) nonlinear sigma model are examples of coset models. A 
coset model is a theory of fields 0(x) which take values in the coset space G/H. Here 
G denotes a compact Lie group and H a closed subgroup of G. The coset space G/H is 
a manifold though only a Lie group if // is an invariant subgroup of G. The coset space 
has dimension dimG/H = dimG - dim//. In 2-dimensional Euclidean space, is a 
map from R^ -^ G/H, whereas at finite temperature this map becomes S xR ^ G/H. 
In a coset model the fields transform as (p{x) -^ g(pix) where g e G. The Lagrangian 
density of a coset model is given by 

X = jlf'g'J{<(>)df,<f>*dy(f>j , (5.7) 

where rj^y is the space-time metric (d^y in Euclidean space) and gij((p) is a metric on 
the coset space G/H. The metric gij is to be chosen in such a way that the Lagrangian 
density is invariant under transformations G. 

Coset models have an interesting connection to the topology of space-time. If 
one identifies all points at space-time infinity, the fields (p{x) become a map from 
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S^ ^ G/H. These maps can be topologically nontrivial, which impHes that there are 
maps (f>{x) which cannot be deformed continuously into each other. Such maps can be 
classified according to the second homotopy group of G/H, n2iG/H). In the case that 
the second homotopy group is nontrivial, the action of the coset model has more than 
one stationary solution, which in Euclidean space-time are called instantons. At finite 
temperature all points at xi = ±00 can be identified, as a result space-time becomes 5^ 
as well. Hence finite temperature instantons (also called calorons) do exist, but unlike 
instantons at zero temperature they satisfy the boundary conditions for scalar fields in 
imaginary time (Lazarides, 1979; Affleck, 1980a; Gross e? a/., 1981; Bruckmanne^aZ., 
2005). 

The O(A^) nonlinear sigma model is an 0{N)/0{N - 1) = S^~^ coset model with 
real scalar fields. Since in a coset model, the fields take their value in the coset space, 
it follows that 0,-0j - N/gl, where j - 1 . . . A'^ and N/gl is the radius oi S^'^. If the 
metric is chosen to be gij = dij the Lagrangian density becomes invariant under G. 
Furthermore the Lagrangian density discussed in Chapter 4, Eq. (4.1), is recovered 
from Eq. (5.7) in this way. Since ;r2(5'^~^) i^ only for N = 3, the 0(A'^) nonlinear 
sigma model has instanton solutions for A'^ = 3 only. 

The CP^-i model is an S\J{N)/\J{N -\) = CP^-'^ coset model. The space CP^^'^ 
is the so-called complex projective space, it is the space of A/^-dimensional complex 
vectors z satisfying the equivalence relation z ~ z' if z = Az' where A e C and 
the relation z^z = c. Hence the A'^ scalar fields 0,(x) of the CP^^^ model satisfy 
0*(x)0,(x) = N/g^. The equivalence relation together with the constraint directly 
translates into the requirement of U(l) gauge invariance. The corresponding metric on 
the coset space can be read off from Eq. (5.3), gij = dij + g^4>i4>*;/^- It turns out that 
7T2{CP'^^^) = Z, hence the CP^^^ model has instanton solutions for any N. 

Instantons 

Classical instantons are solutions to the classical equations of motion in Euclidean 
space-time which have finite action. They can be characterized by an integer topolog- 
ical charge Q which is given by 

Q = ^.j jd^xe,yd, (0*5,0,) . (5.8) 

Using Eq. (5.6) the expression for Q can be written in terms of the gauge field as 
follows 

Q^^fd^xe^yd^A,. (5.9) 

To show that Q is an integer, consider a configuration with finite action. In that case the 
fields should go to a constant times a phase factor (arising from gauge invariance) at 
infinity, lim_j^_^oo (pix) = (pco exp[/o"(x)]. Since Q is the integral over a total divergence. 
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it can be written as a surface integral at infinity. As a result one finds 

1 r^"" do-{R = 00,,?) 1 
Q=2^j dz^ ..a = —Act , (5.10) 



d& In 

where Act denotes the change in cr by going along a circle at infinity. Since the fields 
with finite action should be continuous functions, cr is a multiple of 2n. Hence Q can 
only take integer values. 

Like any vector, C' - D/^cp' satisfies \C' ± ie^vC'y^ > 0. Working out this equation 
gives a lower bound on the Lagrangian density of the CP^^^ model 

X > ±jie^yd^ (<P*dy4>i ) - +--^e^y5^Ay . (5.11) 

Integrating this equation over x results in the classical action S which satisfies the 

following lower bound 

nN 
S>—\Q\. (5.12) 

gh 

For instantons this lower bound turns into an equality. The explicit form of the instan- 
tons can be found by solving the self-duality equation 

D^(Pi = ±ie^yDy<Pi . (5.13) 

The solutions of this equation, the instantons, are of the following form (Golo and 
Perelomov, 1978) 

<l>^(^)-{Nlglfe^-^^^P^^^. (5.14) 

'^ "^ \w{Xq ± lXi)\ 

where w{xo±ixi) is a smooth A^-dimensional complex vector function. The one instan- 
ton (with Q = 1) for example can be written as (Golo and Perelomov, 1978; D'Adda 

etal, 1978) 

Aui + [(xo - ao) - Kx[ -ai)] v, 
4>i(x) ^ ; — rjT , (5.15) 

where a^ is the position of the instanton in space and A the size of the instanton. The 
constants u and v obey the following relations, |mP = N/g^, |vp = N/g^ and u* -v = 0. 
These results can be generalized to finite temperature. If all points at xi = ±00 are 
identified the space changes from S xR. into 5^. To obtain a finite action at non-zero 
temperature, the fields cp still have to go to a constant, limjc^oo (p(x) = (poo exp[/o"(x)]. 
Due to the finite temperature boundary condition cr{xo, x\) = cr(xo +jS, xi) -1- 2nN. The 
topological charge is still quantized at finite temperature (Affleck, 1980a). The action 
has the same bound as at zero temperature. The explicit form of the instanton at finite 
temperature (also called calorons) follow from solving the self-duality equations using 
the boundary condition on xq, see Affleck (1980b). 
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A correspondence between CP^ and 0(3) 

The CP^ model is equivalent to the 0(3) nonlinear sigma model, which was shown 
by Eichenherr (1978), see also Banerjee (1994). This means that there is a one-to-one 
correspondence between those two models. Such an equivalence for example implies 
that one can obtain the pressure of the CP^ model by calculating the pressure of the 
0(3) nonlinear sigma model. The correspondence can be made explicit by writing the 
0(3) nonlinear sigma fields x{x) as follows 

Xaix) = (gllN)"^ (p*{x){cTa)ij 4>j(x) , a = l...3, (5. 16) 

where cTa are the three 2 x 2 Pauli matrices satisfying {o-a,cri,} = 26 ah- Using Eq. (5.16) 
the 0(3) nonlinear sigma Lagrangian density, Eq. (4.1), with constraint turns into the 
CP^ Lagrangian density, Eq. (5.1), with the corresponding constraint. 



5.3 Effective action 

To obtain the effective action, the constraint from Eq. (5.1) can be implemented by 
introducing an auxiliary Lagrange multiplier field a. This gives rise to the following 
Lagrangian density which is equivalent to the original Lagrangian density Eq. (5.1) 

£=^- \D^(l>if -'-a (4>*<Pi - N/gl) (5. 17) 

Since the Lagrangian density is now quadratic in the fields, the Gaussian integration 
over these fields can be performed. This results in the following effective action 

5eff - A^Tr log [-D^Hy - ia) + i—^ f <fxa{x) . (5.18) 

The covariant derivative D^ transforms as U^DfiU. Hence, as expected this effective 
action is invariant under local U(l) transformations. In order to obtain the effective 
potential, ^eff has to be expanded around the vacuum expectation values of the a 
field. The vacuum expectation value of the a field is purely imaginary as is proved 
in Eq. (4.12). Therefore a can be expressed in terms of the sum of its vacuum ex- 
pectation value inp- and a quantum fluctuating field a as follows, a = inp- + a/ ^/N. 
For convenience the gauge fields are rescaled by a factor Va^- This rescaling with 
factors of -0V does not have any effect on the final results since the a and A^ fields 
are integration variables of the path integral. It is just a convenient way to organize the 
1 /N expansion. Working out the square of the covariant derivative and inserting the 
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vacuum expectation value gives the following effective action 



S,, = NTrXog[-d^ + m^- — {d,,A^]^-A,A^-i— 






m 



y/N 



a(x) 



(5.19) 



Now as shown by Affleck (1980a) ^eff contains stationary solutions A^ at finite 
temperature, which vanish at zero temperature according to Witten (1979). These so- 
lutions have a quantized topological charge given by Eq. (5.9). Since these solutions 
are stationary points of an action in which quantum effects are incorporated, they are 
called 'quantum instantons'. In the limit of high temperature, the exact form of these 
quantum instantons can be found (Affleck, 1980a). Such instantons need to be consid- 
ered in a full calculation of the pressure. As a first step to investigate the relevance of 
the quantum instantons to the pressure, only the fluctuations around the trivial config- 
uration A^ = are taken into account in this chapter. The result for the pressure of the 
CP^ model with Q = will be compared to the pressure of the 0(3) nonlinear sigma 
model in which the contribution of all quantum instantons is take into account. Be- 
cause the CP^ model is equivalent to the 0(3) nonlinear sigma model this comparison 
should therefore give the contribution to the pressure of the quantum instantons. 

To obtain the pressure for the Q = configuration, S eft has to be expanded around 
the vacuum expectation value im^ and the trivial configuration A^ = (with Q = 0). 
Expanding S es gives 



5eff = A^Tr log (-52 + m^) + NTr ^^ 



(-1) 



k+i 



f^M^ + Tb--Tfe|^.'^il 



k=i 



2slJ 



^4 



-d"^ + m^ 



m — 



Va^ 



a{x) 



(5.20) 



With help of the following relations 






TrO 


= {x\0\x), 

Jx 


(5.21) 


(x 
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-52 + m2 


y) 


r 1 

(P(x-y) i 

jp P^ + m^' 


(5.22) 


-i{x A^5^ y) 


e'P^^-yp^Af'ix) , 
Jp 


(5.23) 




-i{x ^^A^' 


y) 


e'''^'~y^P,A^'(y) , 
Jp 


(5.24) 



the traces in 5 eft- can be evaluated. To next-to-leading order in l/N these relations lead 
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to 
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-/Tr 



Tr 



1-^2 + ^2 j 



-52 + 



m^ 



Jx Jp 



m^ 



P^ +W?- 



0, 



A^,A^' 



-52 + 



m^ 



-'^v 






(5.25) 
(5.26) 
(5.27) 



Tr 



a 



a 



-52 + rrP- -52 + m2 



- r e'^(->')a(x)a(y) r — -1 2 7^1^ ' (5.28) 



jTrl— r:^ T— r^ t + 



d^ + rrfi -52 + rrfi —d^ + m2 -<92 + ^2 ^ 



Jx. 



e"^'''-y>a{x)A^'(y) 



r ^ 

Jo [(P + 02 + , 



XXP jQ[{P + Q)^+m^][Q^+m^] 



-0, (5.29) 



Tr 



-52 + m2 -52 + m2 



Jxxp Jq [{P + 02 + m^] [Q^ + m^] 

Using the results above one gets the following effective action up to corrections of 
order 1/ Va^ (D'Adda et al, 1978) 

Seff-A^Trlog -52 + m2 --^/3y + /^ S(x) - - 

^ ^2^2 2 Jx \gl Jp P^ + m 

1 [ a{x)Y{x-y)a(y) + l- [ A^{x)A^y{x -y)A'(y) + ol-^], (5.31) 

2 Jxj 2 Jx,y \ ^/NJ 



+ 



■,Y ^ JX,Y 

where the inverse a and gauge field propagators are given by (D'Adda et al, 1978) 



r 1 i_ 



iff^ 



Vp2 + 4^2 + v^ ^ 

2^ V/'2(P2 + 4^2) "° (^ Vp2 + 4^2 _ Vp2^ 



:log 



, (5.32) 
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{Pi. + 2Q^){Py + 2Qy) 

QQ^+m^ Jq [{P + 02 + ^2] ^Q2 + ^2] 



JoQ^+ m^ Jo 



= U, 



'liv 



Pj£^ 

P2 



K{P), (5.33) 
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q(P + Q? + m^ 
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+ (P^ + Am^) 
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2n 



P^ + 4m2 



^2 



log 



e (^ + 2)^ + m2 ^2 + ^2 

V/'2 + 4m2 + VP2'1 



V/'2 + 4m2 - Vp2| 



(5.34) 



The trace of the inverse gauge field propagator was evaluated using dimensional reg- 
ularization in J = 2 - 6 dimensions, in that case (5|^ = 2 - e. Of course using a 
Pauli-Villars regulator one obtains the same result (D'Adda et al, 1978), but one has 
to be careful in applying this regularization technique in this case because it explicitly 
breaks gauge invariance. 

Equation (5.31) shows that although a kinetic term for the gauge fields was ab- 
sent from the classical action, it follows that such a term is generated by the quan- 
tum corrections. The Lorentz structure of the inverse gauge field propagator follows 
from gauge invariance which requires the inverse propagator to be transverse, that is 
P^A^"' - 0. 

The nonzero temperature results can be obtained by changing the integrals over 
momenta into sum-integrals. Due to gauge invariance, the inverse gauge field propa- 
gator is still transverse at nonzero temperature. Its tensor structure at finite temperature 
is the same as at zero temperature (this is typical for 1 -1- 1 dimensions, in more dimen- 
sions Lorentz symmetry breaking terms can appear at finite temperature). At finite 
temperature one finds (Lazarides, 1979) 



r(P) = 



1 



1 



Vp2 + 4rrfi + ^fp^] 



2^ V^2(pi74;;^ (^ ^/p2:^4^ - Vp2 ^ 



:log 



+ TriP) , 



(5.35) 



K(P) - :;- 






P2 + 4m2 r VP2T4^^WP2 
P^ °^ivp2+4m2- VP2, 
where TriP) is up to a factor 2 equal to Eq. (4.44) 

^2 + 2pq 



+ {P^+4m^)TT{P) , (5.36) 



TriP) 



1 rdq . -,, 



2pqY + Api.uj, 



(5.37) 



Q^q 



In the low-momentum limit at zero temperature the inverse gauge field propagator 
becomes 

(5.38) 



^^ * T^ (^'^^^ - P^P^) ■ 
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This is expression is up to a constant (which could be absorbed in the gauge fields 
by a redefinition) equal to the ordinary inverse photon propagator which gives rise 
to a Coulomb potential between two charges. In 1 + 1 dimensions, this potential is 
proportional to the distance between the charged particles. Therefore the CP^'^ model 
has confinement at zero temperature (D'Adda et al, 1978; Samuel, 1983). However, 
confinement disappears at non-zero temperature as was shown by Lazarides (1979), 
see also Actor (1985). 

5.4 Effective potential 

By integrating the effective action ^eff over the quantum fluctuations a and A^ one 
obtains the effective potential. The leading term of the effective potential can be read 
off directly from Eq. (5.31). In order to obtain the next-to-leading order corrections 
a Gaussian integration has to be performed. For the contribution arising from the a 
fluctuations this is straightforward. But to calculate the gauge field contribution to 
the effective potential one has to fix the gauge, because otherwise it is not possible to 
invert the gauge field inverse propagator. In the generalized Lorenz gauge this boils 
down to adding the gauge fixing term 

^£(^A,f, (5.39) 

to the effective action. As a consequence the ghost fields give a contribution to the 
effective potential which is £, log P^. After subtracting T-independent constants it is 
found that the gauge field and ghost contribution to the effective potential together is 
given by 

%auge(m2) - / log P2 - 1/ log det U^y + ^P^Py 

= \^^ogP^-^--^\ogM'^. (5.40) 

which is independent of ^ as expected. The same result could of course also be ob- 
tained in any another gauge, for example the Aq = gauge. In that gauge the ghost 
fields give a contribution of ^j,\og{P^^, whereas the Ai gauge field gives a contribu- 
tion of -^^log[(l - Pj//'^)Aj[^]. Together, this leads of course to the same result as in 
the generalized Lorenz gauge. 

From Eq. (5.31) and the results above, one obtains the complete finite temperature 
effective potential 

^(m^) = N^i^oim^) + ^NLo('«^) , (5.41) 

where the leading order contribution is given by 

^Lo(m2) - ^ - / log(p2 + n?) , (5.42) 
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and the next-to-leading order contribution by 

^NLo(m2) = -i^ log T{P) - i^ log ^^^{P) + i^ log P^ . (5.43) 

The effective potential is ultraviolet divergent. To regulate the divergences an ul- 
traviolet momentum cutoff A is introduced. This yields 



2 2 

m m 



2gl 4;r 






m^ 



+ —T^Jo(/3m) . (5.44) 

47r 



Here 

Jo(fim) = — -^^n(co,) . (5.45) 

T^ Jo t^^/ 

The minimum of the leading order effective potential obeys the following gap equation 

1 1 /A2\ 1 
- = —\og\ — \ + —J,(fim), (5.46) 



where 



Jo ^q 

The gap equation can be rendered finite by the substitution g-j^ -^ Zlg^(n), where 



Jx(J3m) = A \ -^n{ojq). (5.47) 



1 ?2 /^2X 

From this renormalization prescription it follows that the leading order beta-function 
is given by 

2>_ tis'W s* , ^/i 



The negative sign shows that the CP^^"^ model is asymptotically free. 

In order to calculate the next-to-leading order contribution to the effective poten- 
tial, "yNLO is written as a sum of divergent (D) and finite parts (F) in the following 
way 

^NLo(m2) - -J (Di +D2 + F1+F2 + F3+ F4) - ^T^ . (5.50) 

where the divergent and finite quantities are defined through the following relations 

Di + Fi - r log f (P) , F3 = r log f (P) , 

(5.51) 

D2 + F2 = J log a;;(p) , F4 - J log a;;(p) . 
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Here t{P) = 2n ^| P\P^ + 4m^)r{P) and A^(P) = 2n^|P^/{P^ + 4m^)A''^{P). 

The functions Di and D2 are functions which contain all divergences of the next- 
to-leading order effective potential. In order to obtain these ultraviolet divergences, 
the high-momentum limits of r{P) and A'^{P) are needed. In the high-momentum 
approximation {\p\ » T) one finds using Eq. (4.48) 






log hr + 



m^ 



2m^ Im^JiiPm) 



P^ 



P2 



log TT + 



m. 



6m^ Im^JiiPm) 



P^ 



P2 




(5.52) 
(5.53) 



where m^ = m^ exp[-/i(jSm)] and m^ - m^ exp[2 - Ji(J3m)]. 

The divergences Di and D2 can be obtained by integrating Eq. (5.52) and Eq. (5.53) 
over momenta. It is found that 
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, (5.54) 


2 
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A^ log log — T 


/A2\ 
- 2i ■ ^'^ 

\miJ 


-1- 6rrp- log log — r 
\m2/ 


. (5.55) 



From these two equations it can be seen that (through tfi^ and m^) the effective po- 
tential contains temperature dependent divergences. They cannot be renormalized in 
a temperature-independent way. However, by using the leading order gap equation 
(5.46) one can show that m^ and m^ become temperature-independent at the minimum 
of the effective potential. Therefore the quadratic divergence and the divergence pro- 
portional to the li function become temperature-independent at the minimum. As a 
result these divergences can be renormalized at the minimum, as will be done explic- 
itly in the following section. 

The finite functions Fi and F2 will be obtained numerically. In order to calculate 
these functions, the divergences will be written partly in terms of an integral as is 
discussed in Sec. 2.5. This prevents subtracting large quantities which can give rise 
to big numerical errors. The functions Fi and F2 are calculated using the following 
expressions 



Fy 



Fo - 



r I log 



•X' 
•X' 



r(P) 



P I log 



l0g(P2/m2) 

KiP) 



l0g(p2/^2) 



- Irrf" log log ( — 

A2 

- 6m^ log log ( -^ 

mi 



(5.56) 
(5.57) 



here P denotes the principal value integral. At zero temperature it is found that Fi ~ 

2 2 
y-ys and F2 ~ f-ci, where ci « 0.611671457 This is in agreement with the 



92 



Thermodynamics of the CP model 



2.5 

2 

1.5 

1 

0.5 



-0.5 

-1 



1 1 1 1 


1 





/ \ 






/ N 






/ N 






" 


Fi/T^ 


— 


/ 

/ \ 


F2/T^ 




- 




_ 


/ S 






/ s 






s 






' N 






/ ^. 






; 






/ --^ 






/ 






V '' ^^-—^ — — 






^v / ^^^^^ 






K^^L^^"""'''^ 






1 1 1 1 


1 





10 



12 



/3m 



14 



Figure 5.1: The finite functions Fi(JSm)/T^ and F2(fim)/T^. These functions con- 
tribute to the next-to-leading order effective potential. 



zero temperature calculations of Campostrini and Rossi (1992). For convenience the 
finite temperature parts of Fi and F2 are defined as Fi = Fi - m^yEl{2n) and f 2 - 
F2 - m^c\l{2n). These functions divided by T^ depend on jSm only and are displayed 
in Fig. 5.1. For large ySm these functions go to zero because that is equivalent to taking 
the limit of zero temperature. For small jSm they go to zero as well, for the same reason 
as discussed in Sec. 4.6 for the function F2. 

The finite functions F^jT^ and F/^IT^ were calculated using the method explained 
in Sec. 2.5. They are displayed in Fig. 5.2. The fim large limit of f 3 and F4 can 
be obtained by noting that for large jSm the temperature dependent part of the inverse 
propagator does not contribute to F3 and F4. Furthermore, the dominant contribution 
to the difference of a sum-integral and an integral arises from the low momentum 
modes. So the large m behavior of the zero temperature inverse propagators can be 
used to obtain a large fim approximation for F3 and F4. It is found that 



2jp 



n , 



logP^ ^--T^ 



2jp 



p ^ 



(5.58) 



As one can see in Fig. 5.2 this is in agreement with the numerical calculations. 

The small fim limit of F^IF^ and Fa.jF'^ can be obtained as well. For that one 
needs the small fim limit of T{P). That limit can be found by first performing the 
momentum integration and taking the zero modes as is explained around Eq. (4.78). 
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Figure 5.2: The functions F'i{pm)IT^ and F/[{pm)IT^ as a function of /3m. These 
functions contribute to the next-to-leading order effective potential. 



In this limit 



This results in 



np) 



1 



I5m P^ + Arn^p^ 



(5.59) 



(5.60) 



F3 « ^\ogP^ + ^J[\og{P^+4m^)-£\og{P^)^0, 

F4 « ^/log/'2_iriog(p2^4^2^ + Ilog|^-2J«0. (5.61) 

These limits are also in agreement with the numerical calculations as can be seen from 
Fig. 5.2. 



5.5 Contribution of quantum instantons to the pressure 

In the previous section the effective potential was evaluated. It was found that the ef- 
fective potential contains temperature-dependent ultraviolet divergences. At the min- 
imum these temperature-dependent divergences will disappear as will be discussed 
now. To calculate the effective potential at the minimum one only needs to solve the 
leading order gap equation (5.46) as was shown by Root (1974) and as discussed in 
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Sec. 4.5. As a result the leading and next-to-leading order contributions to the pressure 
are given by 



'Pi^o=^lo{ml)-^l^\ml), 



^nlo-^?;lo(4)-^nEo('«o)' 



(5.62) 



where m| is the solution of the leading-order gap equation (5.46) at temperature T. 
By using the leading order gap equation (5.46) it can be shown that at the minimum 
the divergent terms D\ and D2 become 
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(5.63) 



(5.64) 



eT' 



Hence the temperature-dependent quadratic divergence and the divergence propor- 
tional to the li function become temperature-independent at the minimum of the ef- 
fective potential. They cancel in the calculation of the pressure due to the subtraction 
of the zero temperature contribution. The divergences proportional to the log log func- 
tion can be absorbed into the coupling constant ^^. The renormalization factor Zl to 
next-to-leading order becomes in this way 



— = i + ^iog — +-^loglog — +<9 — 
Z] In ^\ii^ Nn ^ ^\fi^ \N^ 



(5.65) 



From this renormalization prescription it follows that the beta-function to next-to- 
leading order in 1 /N is given by 



dyU 



n \n^ 



(5.66) 



Using the results above it follows that the leading and next-to-leading order con- 
tributions to the pressure are given by 
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(5.67) 



(5.68) 



The results of the calculation of the pressure are displayed in Fig. 5.3 for the ar- 
bitrary choice g^ijj. = 500) = 10, for different values of N. As one can see, for low 
temperatures and all finite values of N the pressure first decreases for increasing values 
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Figure 5.3: Contribution of the zero winding number configurations to the pressure P 
of the CP^~^ model normalized to A'^r^, as a function of temperature, and for different 
values of N. 



of T. A decreasing pressure is unphysical because it implies that the entropy which is 
the derivative of the pressure with respect to temperature becomes negative. But this 
is in disagreement with the third law of thermodynamics that states that the entropy 
is minimal at zero temperature. If one believes that the 1 /N expansion is a good ex- 
pansion (which is for example shown in the 0{N) nonlinear sigma model discussed 
in Chapter 3) it is likely that the reason for this negative pressure is that the effective 
action, Eq. (5.20) was only expanded around the vacuum A^ = solution with zero 
winding number. The contribution of other vacua (quantum instantons) with nonzero 
winding number was left out of the calculations. As one can see from the figure, the 
problem of the negative pressure becomes less severe if A'^ becomes larger. This is in 
agreement with the fact that the instanton contribution vanishes in the N ^ oo limit 
(Witten, 1979) as is discussed in Sec. 5.1. Moreover, that this problem arises at low 
temperatures is also reasonable because instantons are non-perturbative in g. They 
become less important at small couplings, so they should vanish at high temperatures 
due to asymptotic freedom. 

Using the equivalence between the 0(3) nonlinear sigma model and the CP^ model 
(see Sec. 5.2) the exact contribution of the quantum instantons with nonzero winding 
number to the pressure can be found. Since there are no gauge fields in the 0(3) 
nonlinear sigma model and because the integration over the scalar fields can be done 
exactly the effects of all quantum instantons are automatically included its effective 
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Figure 5.4: Pressure f of the 0(A'^) nonlinear sigma model to NLO in \ IN iov N - ?> 
compared to the contribution to the pressure of the configuration with zero winding 
number Q of the CP^ model. The pressures are normalized to T^ and displayed as a 
function of temperature, for the g^(jj. - 500) = 10. 

potential (see Chapter 4). Due to the equivalence, the pressure of the CP^ model 
should be exactly equal to that of the 0(3) nonlinear sigma model. In Fig. 5.4 the result 
of the NLO calculation of the pressure of the 0(A'^) nonlinear sigma model iov N - 3 
is compared to the contribution to the pressure of the £.P^ model with zero winding 
number. It can be seen that for very low and high temperatures both pressures coincide. 
For intermediate temperatures these pressures differ. This difference is displayed in 
Fig. 5.5. This is a strong indication that quantum instantons give a sizable contribution 
to the pressure and other thermodynamical quantities at where the pressure (divided 
by r^) increases considerably. 



5.6 Summary and Conclusions 

In this chapter the effect of quantum instantons on the thermodynamical quantities was 
investigated. For that the effective potential of the CP^^^ model expanded around a 
background with zero winding number was calculated to next-to-leading order in 1 /N. 
It was shown that the effective potential contains temperature-dependent divergences 
which only can be renormalized at the minimum of the effective potential. Hence 
thermodynamical quantities (which are all defined at the minimum of this effective 
potential) can be rendered finite like in the (non)linear sigma models discussed in 
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Figure 5.5: Contribution of the configurations with nonzero winding number Q to the 
pressure P of the CP' model. The pressure is normalized to T^ and displayed as a 
function of temperature. 



Chapters 4 and 6. 

It was found that for all values of A'^ the contribution of the vacuum with zero 
winding number gives rise to a negative pressure, which is large for intermediate tem- 
peratures where the leading order pressure divided by T^ increases strongly. Since this 
is unphysical, it indicates that the quantum instantons must give a considerable contri- 
bution to the pressure for intermediate temperatures. In agreement with the vanishing 
of the instantons in the N ^ oo limit this problem of the negative pressure becomes 
less severe for large values of N. 

For N = 2 the exact contribution of the quantum instantons was found using an 
equivalence between the 0(3) nonlinear sigma model and the CP^ model. In the 1/A^ 
approximation to the 0(3) nonlinear sigma model one implicitly integrates over all 
quantum instantons and one finds a finite pressure in next-to-leading order in l/N. 
Comparing this result to the pressure of the CP^ model with zero winding number 
gives the exact contribution of the quantum instantons for N - 2. They give a large 
contribution for intermediate temperatures where the pressure divided by T^ raises 
quickly. 

A possible extension of this work would be to add a term {^Q to the effective action 
and investigate the dependence of the thermodynamical quantities on ??. Moreover one 
could try to take the fluctuations around the quantum instantons with nonzero winding 
number into account explicitly and see whether they will reduce the problem of the 



98 Thermodynamics of the CP^ ' model 

negative pressure, but this is not attempted in this thesis. 



Chapter 6 

Thermodynamics of 0(N) sigma models in d=3+l 



The O(A^) linear and nonlinear sigma model in 3 + 1 dimensions are low-energy effec- 
tive theories for QCD. The pressure of these two low-energy effective models will be 
computed up to next-to-leading order in l/N by minimizing the effective potential. It 
is found that the finite temperature effective potential contains temperature dependent 
divergences which become temperature-independent at the minimum. The calculated 
pressure can serve as a prediction for the pressure of QCD at low temperatures and will 
be compared to approximations discussed previously in the literature. Finally a mass 
bound on the sigma meson is presented. This chapter is based on: Thermodynamics 
of 0{N) sigma models: l/N corrections, J.O. Andersen, D. Boer and H.J. Warringa, 
Phys. Rev. D70 116007, (2004). 

6.1 Introduction 

Although the QCD Lagrangian possesses a chiral symmetry in the limit of zero quark 
masses, the true QCD ground state does not respect this symmetry. The chiral symme- 
try is spontaneously broken by quantum effects. To be specific, QCD with Nf massless 
quarks has a global S\J{Nf)L x S\J(Nf)R symmetry, which for the ground state at low 
temperatures is broken down to an S\J(Nf)v symmetry. According to Goldstone's 
theorem, there is a massless, spinless particle for each generator of a broken global 
continuous symmetry. In this case this implies the occurrence of Nl - 1 Goldstone 
bosons. In phenomenological applications Nf is either two or three, and one also has 
to take into account the explicit symmetry breaking due to the nonzero quark masses. 
Both the spontaneous and the explicit chiral symmetry breaking are apparent in the 
low-energy hadronic particle spectrum, where the expected number of relatively light 
mesons is observed (e.g. the three pions for Nf = 2). At sufficiently high temperatures 
one expects the chiral symmetry to be restored and lattice simulations of QCD sug- 
gest that this happens at a temperature of approximately 170 MeV depending on the 
number of quarks and their masses. 

In the case of two flavors, the situation is the simplest, since one can exploit the 
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fact that the SU{2)l x SU(2)« symmetry is locally isomorphic to S0(4). If baryons 
(nucleons in this case) and the p and a> meson are not included the simpler 0(4) linear 
sigma model can be used as a low-energy effective theory for describing the dynamics 
of three pion fields and one sigma field. These four fields form a four-dimensional 
vector (p in the fundamental representation of 0(4). At low temperature, the 0(4) sym- 
metry is spontaneously broken down to 0(3), where the sigma field acquires a vacuum 
expectation value and the three pions are interpreted as the Goldstone bosons. At high 
temperatures the vacuum expectation value of the sigma field vanishes. According to 
the calculations performed in this chapter, the symmetry is restored via a second-order 
phase transition in agreement with other calculations in the literature. If the symmetry 
is restored, the pion and the sigma field have the same mass. 

For Nf > 2 there is no connection between the S\J{Nf)L x S\J(Nf)R model and 
the 0(A'^) linear sigma model since the symmetries differ. The 0(A'^) sigma models 
have besides their relevance to low-energy QCD also other applications. For example 
the Higgs field in the Standard Model is described by an 0(4) linear sigma model. 
In that case the scalar fields are also coupled to the electroweak gauge fields and the 
fermions. The spontaneous symmetry breakdown of the 0(4) symmetry causes the 
Higgs field to acquire a vacuum expectation value. As a result the weak vector bosons 
and the fermions (except for the neutrinos) become massive. Another application is in 
cosmology. It is believed that during early times the universe went through a period 
of rapid expansion, called inflation. The field which controls the inflation, called the 
inflaton can be described by a linear sigma model (see for example Kolb and Turner 
(1990)). In condensed matter physics the O(A^) sigma models are used to model spin- 
spin interactions (see for example Itzykson and Drouffe (1995)) . 

In this chapter both the 0(A'^) linear sigma model and 0(A'^) nonlinear sigma model 
in 3-1-1 dimensions will be studied at finite temperature and to next- to-leading order in 
the 1 /N expansion. At zero temperature, the 1 /N expansion was applied to 0(A'^) sigma 
models a long time ago at leading order (LO) by Coleman et al. (1974) and at next-to- 
leading order (NLO) by Root (1974). At finite temperature the LO 1/A'^ contribution 
has been studied by Meyers-Ortmanns et al. (1993). In that case, the effective potential 
is that of an ideal gas and thus straightforward to compute. The NLO 1 IN corrections 
are less trivial, since they involve a momentum-dependent self-energy and cannot be 
evaluated analytically. A high-temperature expansion was performed by Jain (1993) to 
obtain purely analytical results for the linear sigma model. Similarly Bochkarev and 
Kapusta (1996) resorted to a "high-energy" approximation, which makes the calcula- 
tions manageable. However, this approximation is uncontrolled and it is difficult to 
assess how reasonable it is, unless one calculates the full NLO 1/A'^ corrections. These 
corrections will be calculated in this chapter and compared to the "high-energy" ap- 
proximation. 

The O(A^) sigma models have also been studied in detail at finite temperature using 
various other approaches. A systematic study has been carried out by Chiku and Hat- 
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suda (1998) using optimized perturbation theory. The method was used to calculate 
spectral functions, properties of the elfective potential, and dilepton emission rates. 
The 2PI formalism (Cornwall et at, 1974) (see Sec. 3.1) has also been used to exam- 
ine various properties of the 0(A'^) linear sigma models at finite temperature (Baym 
and Grinstein, 1977; Amelino-Camelia and Pi, 1993; Amelino-Camelia, 1997; Roh 
and Matsui, 1998; Petropoulos, 1999; Lenaghan and Rischke, 2000; Nemoto et ai, 
2000; Aarts and Martinez Resco, 2004; Arrizabalaga et ai, 2005; Roder, 2005), see 
Petropoulos (2004) for a recent review. For example, in several papers the temperature 
dependence of the pion and sigma masses, and of the vacuum expectation value of the 
sigma field, have been investigated. The calculations of the scalar field effective po- 
tential as a function of temperature have been carried out in the Hartree approximation 
and the large-A'^ limit (Amelino-Camelia and Pi, 1993; Amelino-Camelia, 1997; Roh 
and Matsui, 1998; Petropoulos, 1999; Lenaghan and Rischke, 2000; Nemoto et ai, 
2000). In these cases, the gap equations for the propagators are easy to solve since 
the self-energy reduces to a local mass term. In the Hartree approximation (which 
implies ignoring momentum-dependent self-energies), the result has been shown to 
be problematic (and a first-order phase transition occurs), which has been remedied 
by including more diagrams in the truncation (Verschelde and De Pessemier, 2002; 
Baacke and Michalski, 2003), resulting in a second-order phase transition. If one goes 
beyond the Hartree approximation or includes the full NLO contributions in the 1 /N 
expansion, the gap equations become nonlocal and very difficult to solve. 

The advantage of using the IPI formalism over the 2PI formalism is that one does 
not have to make a Hartree or other approximation to calculate the full 1 /N corrections 
to the thermodynamical quantities. However, to describe out of equilibrium phenom- 
ena, the 2PI formalism is favored as is explained in Sec. 3.1. 

The approach followed in this chapter is similar to the study of the nonlinear sigma 
model in 1-1-1 dimensions (see Chapter 4) and to that of the CP^^"^ model in l-i-l di- 
mensions (see Chapter 5). Similar conclusions about the renormalization of the effec- 
tive potential in 3-1-1 dimensions as in l-i-l dimensions will be drawn. It turns out that 
at NLO, temperature-independent renormalization is only possible at the minimum of 
the auxiliary field effective potential. This aspect of the 1 /N expansion was missed in 
previous work by Jain (1993) and Bochkarev and Kapusta (1996), since the renormal- 
ization is considerably simplified or even ignored in the various approximations. 

Since explicit chiral symmetry breaking plays a very important role in the actual 
hadron spectrum at low energy, the case of explicit symmetry breaking will also be 
considered in this chapter. The results change considerably and moreover, a critical 
temperature cannot be determined in that case, since the second-order phase transition 
turns into a smooth cross-over. 

The nonlinear sigma model in 3-1-1 dimensions is non-renormalizable and should 
be viewed as an effective theory, which is valid up to a certain energy scale where new 
physics enters. Strictly speaking the linear sigma model is renormalizable, but since it 
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becomes a trivial theory in the limit where the cutoff goes to infinity (see for example 
Amelino-Camelia (1997)), it will be treated as a theory with a finite cutoff. Given a 
finite cutoff, terms are called divergences when they are increasing in magnitude with- 
out bound as the cutoff is increased. The low-energy physics should be independent of 
such terms (decoupling) and they can be subtracted in the renormalization procedure 
in order to avoid increasing sensitivity to the ultraviolet cutoff as it grows. On general 
grounds, one expects the temperature dependence to be insensitive to an increasing 
cutoff due to the exponential suppression provided by the Bose-Einstein distribution 
function. Therefore one expects the renormalization to be possible in a temperature- 
independent way. 

This chapter is organized as follows. In Sec. 6.2, the effective action of the linear 
and nonlinear sigma model in the l/N expansion are discussed. In Sec. 6.3, the effec- 
tive potential and gap equations to next-to-leading order are calculated. In Sec. 6.4, 
the results for the pressure at next-to-leading order for general A'^ and for the special 
case of N = 4 are presented. Also, the so-called high-energy approximation is dis- 
cussed and compared with exact numerical results. Sec. 6.5, is devoted to the choice 
of parameters for N = 4,m order to make contact with low-energy QCD phenomenol- 
ogy. A bound on the scalar sigma meson is derived too in this section. In Sec. 6.6, a 
summary and conclusions are given. 

6.2 Effective actions 

The Euclidean Lagrangian density of the 0(A^)-symmetric linear sigma model with a 
symmetry breaking term proportional to H is given by 

£=\ (d^4>if + ^ (<P4i -Nfl,f - ^Hcf>N , (6.1) 

where / = 1 . . .N. Summation over repeated indices is implicitly understood. The 
subscript b denotes a bare quantity. The coupling constants are rescaled with factors 
of A'^ in such a way that for large N the action naturally scales as N as is explained in 
Sec. 3.2. 

It is possible to eliminate the quartic interaction term from Eq. (6. 1) by introducing 
an auxiliary field which is denoted by a, in order to allow for Gaussian integration. To 
this end one adds to the Lagrangian density Eq. (6.1) the term 



£a = — 

2Ah 



iAh 



-iS(^.-^.-A^/.y 



(6.2) 



such that one has 



£=\ (d^,4>f - \a (<P4i - Nfl,) + ^a' - <^mn ■ (6.3) 
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If one integrates over the a fields, the original Lagrangian in Eq. (6.1) is recovered. 
Hence Eq. (6.3) is equivalent to the original Lagrangian density of the linear sigma 
model. In the limit Ah ^ oo the term quadratic in a vanishes and one obtains the 
Lagrangian of the nonlinear sigma model Eq. (4.1) with a symmetry breaking term. 

If explicit symmetry breaking is absent (H = 0), the field acquires a vacuum 
expectation value by spontaneously breaking the symmetry. Because of the residual 
0(A'^- 1) symmetry, it is possible to write 4> = i^i,^2,- ■ ■ , ^N-i, o"), such that only cpf^ = 
(T has a nonzero expectation value. For H > the same argument applies, because the 
action is minimal when the cr field is the only one that acquires an expectation value. 

Integrating over the ;r's gives the following effective action 

Seff - ^{N - l)Trlog {-d^ - ia) + f dr fd\ ^ {d^af - '-acr^ 



2 •^'^'^ lAb 



. (6.4) 



The scalar fields cr and a can be written as a sum of space-time independent vacuum 
expectation values irn^ and o", and quantum fluctuating fields a and o" 



a = im + 



a 



Va^ 



cr 



Va^, 



cr + cr . 



(6.5) 



Equation (6.2) can be used to show that the vacuum expectation value of a is purely 
imaginary (for a proof, see Sec. 4.4). The vacuum expectation value of cr is propor- 
tional to VA^, which follows from Eq. (6.1). Substituting Eqs. (6.5) into Eq. (6.4), the 
effective action ^eff can be written as 



>eff 



-(A^ - l)Tr log f -(9^ -F m^ - -^ I - /3VNHd- 



-^f2 l2_Jl 
I'^'^y ^/N) 2Ah\- ^/N 



-—im^ — ^1 - V/7//(T 



(6.6) 



Expanding Eq. (6.6) in powers of 1/ V^ up to corrections of order 1/ VA^, one finds 



Nm'^ _ I- 

- — NHcr + yN X terms Imear m a and cr 

lAb 

2 Jp \ -icr P^ + m'^ I 



(6.7) 
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where ;^'^ = {a{P), d-{P)) is a vector containing the Fourier transforms of a and o", and 
the function Il{P, m) is given by 



n(P, m) = { ^ " . (6.8) 

*T„ (2 + m < P ->- 11'^ -I- "1^ 



■m^{P+ 2)2 + m2 



6.3 Effective potential and gap equations 

One can obtain the effective potential to next-to-leading order in the l/N expansion 
from Eq. (6.7) by performing the Gaussian integral over the fluctuating fields a and o". 
Up to corrections of order 1 /N, the effective potential can be written as 



^(m\ d-) = N^i^oim^, (T) + ^NLoC/w^ cr) , (6.9) 



where 



^Lo(m\a-) = 'l-[fl^-a-^) + ^^+Hcr-^^\ogiP^ + m^), (6.10) 
^nlo('«^'t) - -j{ log I{P,m). (6.11) 

Here, 

I{P, m) = \6n^n(P, m) + + — . (6.12) 

Ah P^ -I- m^ 

To derive the effective potential, divergent constants which are independent of a, m 
and the temperature were subtracted. Equivalently, these terms can be removed by 
adding a vacuum counterterm to the effective potential. In the following, such terms 
are simply dropped. 

In thermodynamic equilibrium, the system will be in the state that extremizes the 
effective potential with respect to m^ and o". This extremum can be found by differen- 
tiating the effective potential with respect to m^ and a, which gives 

dn(P,m) 2d-^ 



S 1 2m2 1 r dm2 (p2+m2)^ / , ,x 

2 r 1 

N%T^m 



2 -J ^ 1 

\m + 



1 

^n(/',m) + f + J^ 

^ ' ^ h P^+rr, 



2 

\ 
o- = H . (6.14) 



These equations are often referred to as gap equations. Solving the gap equations 
gives m and o" as a function of the parameters f„, H and A, and of the temperature. The 
solution of the gap equation is needed to calculate thermodynamical quantities Uke the 
pressure. 
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The inverse a and a propagators can also be obtained from Eq. (6.7). For this the 
2x2 matrix in Eq. (6.7) has to be inverted. As a result one finds 

Dl\P,m) = ln(P,m) + l+ "^ (6.15) 

2 Ah P^ + m^ 

2o-2 
Dg}{P,m) - P' + m^+ ■ (6.16) 

ll{P, m) + 2/ Ah 

In these equations the values for m^ and o" are determined by solving the gap equations. 
The 71 propagator is equal to (see Eq. (4.11)) 

Expanding a around its vacuum expectation value im^ in the previous equation, and 
using the expression for a propagator, Eq. (6. 15), it follows that the inverse n propa- 
gator in momentum space is given by 



-,-1/ 



K 



1 1 



D,\Pm) = P' + m' + -2 -^-T^, T ; tt^ ■ (6-18) 

Gap equation (6.14) can be used to show that in the broken phase where a i^ (for 
H = 0) the pion propagator has a pole at P^ = 0. This implies that also to next-to- 
leading order in l/N the pions are massless, in accordance with Goldstone's theorem. 

From Eq. (6.16) it follows that in the unbroken phase, the cr mass becomes equal 
to the leading order mass of the n field, which is m^. It may appear therefore that 
the cr and n masses are not equal at next-to-leading order in the unbroken phase, but 
this is not a correct conclusion. This is because the cr field only starts to propagate 
at next-to-leading order, so its 1 /N mass corrections require a next-to-next-to-leading 
order calculation. In the calculation of the pressure to next-to-leading order one only 
needs the leading order masses as will be explicitly shown below in Eq. (6.51). 

The leading-order and next-to-leading order contributions to the effective potential 
in 3 -I- 1 dimensions will be explicitly calculated in the next subsections. The integrals 
over momentum are regulated using an ultraviolet momentum cutoff A. Furthermore 
it is assumed that the cutoff is large compared to the other scales in the problem, that 
is A » m, InT. 

Leading-order contribution 

The leading-order contribution to the effective potential is 

m2 / 2 A2 ,\ r4 



-I- 



4 

m 



64;r2 



327r2 /a2\ 1 



+ Hd-, (6.19) 
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where the function Joifim) is 



Jo(fim) = 



32 r°° p^ 

3r4 Jo ^~p 



n{cOp) . 



(6.20) 



Here, n(a>p) = [exp(JSa>p) - 1]"^ is the Bose-Einstein distribution function. Since a 
finite effective momentum cutoff A was introduced in the theory, one makes an error 
in the evaluation of Jq by summing over all Matsubara modes and integrating up to 
infinite momenta instead of up to A. However, this error is negligible as long as A » 
m, 2nT. This remark also applies to the functions /i , K^ and Kf defined below. For an 
investigation of how one could apply a finite cut-off in the calculation of sum-integrals 
see Amte and Rosenzweig (1993). 

Equation (6.19) contains ultraviolet divergences in the sense explained in Sec. 6.1. 
These divergences can be dealt with by defining the renormalized parameters /^ and 
A as 



/2 = fl,-A^/l6n\ 
32n^ , /A2\ 327r2 



(6.21) 
(6.22) 



where A = A(ji). The renormalization group equation for the running coupling A that 
follows from Eq. (6.22) is 



_ dA _ A"- 
dju \6n^ 



(6.23) 



The y6-function is exact to all orders in A^ in the large-N limit, but differs from the 
perturbative one obtained at one loop, see Peskin and Schroeder (1995). However, at 
next-to-leading order they agree as will be shown in the next subsection. After this 
renormalization, the leading-order effective potential becomes 



^LO = ^(/.^-C.V 



m 



647r2 



327r2 ///2\ 1 

—^Hm^ri 



+ 



64n^ 



Jo(fim) + Hd-. (6.24) 



Since the potential term in the Lagrangian should always have a minimum in order 
to have a stable theory. Ah must be positive (cf. Amelino-Camelia (1997) for a detailed 
discussion). From Eq. (6.22) it immediately follows that there is a maximal value for 
the cutoff given by 

/I6;r2^ 

Amax = /^ exp — — 



(6.25) 



Therefore the linear sigma model should be viewed as an effective theory, which is 
at most valid up to the cutoff given by Eq. (6.25). Taking the cutoff to infinity is 
equivalent to taking A to zero, which implies that the theory is trivial. One should 
keep in mind that the renormalized leading-order effective potential does not depend 
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Figure 6. 1 : Leading order sigma condensate iov H -0 (solid line) and H i^Q (dashed 
line). The parameter choice is discussed in Sec. 6.5. 



explicitly on A, but is only valid for m and T much smaller than Amax- When A = Amax 
the linear sigma model reduces to the nonlinear sigma model, since in this case A}, = oo. 
The leading-order renormalized gap equations follow from differentiating Eq. (6.24) 
with respect to m^ and a and are given by 



G = 167rV^ 



H 



m cr , 



where 



2,^,2 



G = T'Ji(J3m) + \6n'(T' - m" log I ^) - ^^" "" 



A 



Here, the function /i (J3m) is defined as 

o noa 2 

7i06m) = — dp^ 

T^ Jo Wn 



n{Dp) . 



(6.26) 
(6.27) 



(6.28) 



(6.29) 



In Figs. 6.1 and 6.2 the vacuum expectation value of the sigma field and the pion 
mass are respectively displayed as a function of temperature. The cases with {H + 0) 
and without (// = 0) explicit symmetry breaking are considered. Both figures were 
calculated using the parameter set discussed in Sec. 6.5. 
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Figure 6.2: Leading order pion mass for // = (solid line) and H i^ (dashed line). 
The parameter choice is discussed in Sec. 6.5 



If H - 0, one can show by using the gap equation (6.27) that either m = or 
o" = 0. From the gap equation (6.26) it follows that for m = the expectation value of 
(T has the temperature dependence 



'^=a//^-^- 



(6.30) 



The order parameter for symmetry breaking, o- vanishes continuously. Hence at T = 
Tc = Vl2/;r there is a second-order phase transition. Below Tc the 0(A'^) symmetry 
is broken spontaneously to 0(A^ - 1) since & i^ 0. Above T^ the 0(A'^) symmetry is 
restored and one has o- = and m i^ 0. 



Next-to-leading order contribution 

In this subsection it will be shown that it is not possible to renormalize the next-to- 
leading order effective potential in a temperature-independent way. It turns out that one 
can only renormalize the effective potential at the minimum, since the temperature- 
dependent divergences become temperature independent by using the leading order 
gap equations. Hence physical quantities like the pressure, which can be obtained 
from the minimum of the effective potential can be renormalized consistently. To 
show this, the divergent parts of the effective potential will be extracted. This can be 
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done analytically. A perhaps more familiar example in which the effective potential 
is meaningless outside the minimum arises at zero temperature in gauge theories. In 
this case, the effective potential depends on the gauge-fixing condition except at the 
minimum (Jackiw, 1974; Nielsen, 1975; Kobes et ai, 1991). 

In order to isolate all divergences, one in principle needs to evaluate Il{P, m) in- 
cluding corrections of order m^/A^, since such terms can also give rise to divergences 
in the effective potential. However, since the linear sigma model is an effective theory, 
Eq. (6.1) should be viewed as the part containing only the so-called relevant operators, 
see Polchinski (1984). For instance, irrelevant operators of mass-dimension six are 
not included in the Lagrangian density. Such operators also contribute to Ii{P,m) at 
order 1/A^. Therefore, for consistency with Eq. (6.1), order 1/A^ terms in IV{P, m) are 
not considered. Although it is possible to obtain an exact analytic expression for the 
zero-temperature part of IV{P, m), here it will only be given up to order 1/A^ since this 
expression is much less complicated. As a result one finds 



n(P,m) 



1 



16;r2 



,A2\ //'2 + 4^2 fV/'2 + 4m2- VP2^ 

lQg| — l + l+V P2 ^°^ 



m^ 



VP^ + 4m2 + VP2| 

+ I{T{P,m), (6.31) 



where the temperature-dependent part of Ii{P,m) equals (Bochkarev and Kapusta, 
1996) 



1 /^co 

nr(P,m)--— d^^log 
8;r^7? Jo ^q 



q^ + pq + A^ 



q^ - pq + A2 



n{a)q) . 



Here 



P4 + 4m^pl 
4P2 



In the limit P ^ m,T, IIt{P, m) can be approximated by 



nriPm)^.±, 



r2 Am'T'p' (3P^-4p^)T^ 



(6.32) 



(6.33) 



(6.34) 



The next-to-leading-order effective potential has only ultraviolet divergences. Us- 
ing the leading order renormalization of Ai„ it is easily seen that I{P,m) (defined in 
Eq. (6.12)) becomes finite. Also, the difference 



log I{P, m) , 



(6.35) 



is finite (cf. Sec. 2.5). Therefore, all possible divergences of "yNLO can be isolated by 
calculating 

-^floglnuiPm), (6.36) 
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where Ium{P, ni) is the high-momentum (HM) approximation to I{P, m). It gives the 
large-P behavior of I{P, m). After averaging over angles, one finds 



log/HM ^logC 



, i£i __L (^f '£i, 



P^Ci 2P^\Cij P^Ci 



where 



Ci 
C2 



1 + iogK 



-2m^ 



m^ 



C3 = +m^ 



21og|^|-l 



m^ 



+ 327rV^ + 2r2/iOSm), 



m^ [327r2o-2 + 2T^JiQ3m)] 



(6.37) 

(6.38) 
(6.39) 
(6.40) 



By integrating the function log /hm over P, all the divergences of the NLO eff"ective po- 
tential can be found. The logarithmic and power divergences are given by the quantity 
D, which is 



D 



1 



167r2 



2\+32n'-/Ah 



A^e 



-m^A^ 



1 



,l+32ff2/ii, 



1 + 2ei+32;rVA li 



A^ 
G -F 2m'' log I — 

m^ 



1 



,l+32;r2/ift 



(6.41) 



while the terms that have a small cutoff dependence through their dependence on Ab, 
are given by the quantity E, which is 



1 



167r2 



3„r(-G.|M>og(l.^).(G-2„,f-J^ 



(6.42) 



Since G depends explicitly on the temperature, it is impossible to renormalize the next- 
to-leading-order effective potential in a temperature-independent way. However, at the 
minimum, the leading-order gap equation (6.26) can be used to show that G = ISn^f^. 
Hence, the divergences become independent of the temperature at the minimum and 
renormalization can be carried out in a temperature-independent manner. This will be 
discussed in more detail next. 

The divergence proportional to G in Eq. (6.41) is independent of m in the minimum 
even though it depends on f,r- This divergence can be removed by vacuum renormal- 
ization. The divergent terms which are proportional to m^ can be removed by defining 
the renormalized parameter f,r as 



fn - fn,h \^'^ KJ 



A2 



1 A^ 



N I6n^ N4n^ 



,\+32n^/At 



1 



gl+32;r2/^6 



(6.43) 
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The remaining divergence is proportional to m^ and is removed by renormalizing A}, 
as follows 

From Eq. (6.44), the /3-function governing the running of A can be obtained 

which coincides with the standard one-loop jS-function in perturbation theory, see Pe- 
skin and Schroeder (1995). One could argue from this renormalization that one can 
only trust the 1/A'^ expansion for A/^ » 8. Although the 1/A'^ correction to the y6 function 
indeed has a large coefficient, this is not the case for the effective potential itself as will 
be shown below. 

As mentioned, the terms in E have a small cutoff dependence through their depen- 
dence on X},. These terms will not be renormalized, since they do not grow without 
bound with increasing cutoff and are not strictly speaking divergences. The effective 
potential does not become increasingly sensitive to them with increasing cutoff. The 
term in the first line of Eq. (6.42) becomes smaller if A is increased and the absolute 
value of the other term from E increases as a function of A, but is bounded by a fi- 
nite number which is independent of A}, and A. Moreover, renormalizing these terms 
would invalidate the 1 jN expansion, because of their magnitude. This is similar to 
ordinary perturbation theory, where one is only allowed to do finite renormalizations 
that do not invalidate the perturbative expansion. A final reason for not renormalizing 
these terms is the connection with the nonlinear sigma model {Ah = oo). In that case, 
the terms from E are not divergent and fj^ will be renormalized just as in Eq. (6.43) 
with Ai, = oo. 

The next-to-leading order correction changes the critical temperature T^. Since 
the next-to-leading order gap equations are complicated, it is not possible to obtain an 
analytical expression for the critical temperature at next-to-leading order. However, in 
the limit of small Ah and H = the gap equations simplify to 



i 



From the gap equations it follows that the critical temperature at NLO is 



Tc = a/ — i^A . (6.48) 

^' 1 + 2/N 

This result is the same as obtained by Jain (1993) and Bochkarev and Kapusta (1996). 
It is probably only correct in the weak-coupling limit and Tc may depend on A at NLO 
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in 1 /N. In the next section it can be seen that the transition at NLO remains second 
order. 



6.4 Pressure 

Like in Chapters 4 and 5, the pressure P{T) is defined as the value of the effective 
potential at the minimum at temperature T minus its value at the minimum at zero 
temperature. As was shown in the previous section, it is possible to renormalize the 
effective potential at the minimum. The pressure is therefore a well-defined quantity. 
In order to determine the value of the next-to-leading order effective potential in the 
minimum, only the gap equation to leading order Root (1974) is needed. Writing the 
solutions to the gap equations as 

'"^ ^ '"lo + '"nlo/^ ' (6-49) 

o- = o-Lo + o-NLo/N , (6.50) 

and Taylor expanding the effective potential (6.9), one obtains (up to 0{l/N) correc- 
tions) 

^{m^,a-) - A^n^L0('«L0'0"L0) + ^NL0('«L0''5"L0) 



2 d^LoJm') 



+ 0"NLO 



do- 



0"=0-L0 



The last two lines of Eq. (6.51) vanish by using the leading-order gap equations. In 
the following, the pressure P will be written as 

^ = A^^LO + ^NLO . (6.51) 

From the discussion above, it follows that 

^LO - ^lo{ml,cTT)-^l^\ml,cTo), (6.52) 

^NLO - %^o{m^T'CrT)-'Vll^{ml,cro), (6.53) 

where m^ and d-j are the solutions of the leading-order gap equations (6.26) and (6.27), 
at temperature T. 

In the following, the results for the numerical evaluation of the leading and the 
next-to-leading order contributions to the pressure for general N will be presented. 
Subsequently the N = 4 case which is of relevance for QCD phenomenology will be 
discussed. 

As will be motivated in Sec. 6.5, the following values for the parameters will be 
used: AQu = 100 MeV) = 30, fjr = 47 MeV (note that fj^ as defined in Eq. (6.1) is 
1 /2 times the more conventional definition) and if there is explicit symmetry breaking 



6.4. Pressure 



113 



0.12 



0.1 - 

0.08 - 

jA 0.06 - 

0.04 - 



"1 r 



0.02 - / 



-' 



_L 



"1 r 



J L 



J L 



_L 



"1 r 



J L 



H = Q 
//t^O _| 

I 



50 100 150 200 250 300 350 400 450 500 

T (MeV) 

Figure 6.3: Leading-order pressure !Plo normalized to T^ as a function of temperature 
without and with explicit symmetry breaking. 



H - (104 MeV)^ in order to reproduce the physical value of the pion mass at T = 0. 
A realistic choice of parameters would allow to compare to lattice QCD simulations 
of the Nf = 2 case for T <Tc, although this would require extrapolation of the lattice 
results down to the actual pion masses. 

Leading-order contribution to the pressure 

The leading order pressure can be obtained from Eqs. (6.24) and (6.52). In Fig. 6.3, the 
leading-order pressure normalized by T'^ is shown. If // = 0, the pions are massless 
below Tc. Then it can be see from the leading-order effective potential (6.24) that the 
pressure becomes equal to the pressure of an ideal gas of massless particles: !Plo - 
n^T^ /90. If H i^ 0, the pions are massive. In the limit of zero temperature the !Plo/^^ 
goes to zero because of the Boltzmann suppressing factor exp(-m/r). 

Next-to-leading order contribution to the pressure 

To calculate the next-to-leading order contribution to the pressure, !Pnlo is decom- 
posed as follows 

!Pnlo - Dimr) - D(mo) + Fi + F2 , (6.54) 

where D(m) is the term containing logarithmic and power ultraviolet divergences given 
in Eq. (6.41), and Fi and F2 are finite terms defined below. 
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The term Fi is defined by 



2Jp 



log 



U{P, mT) + — + 



IcTj 



h P^ + 



mt. 



-log 



2 2d-l 

n(/',mo) + — + 



h P^ + 



mr, 



- Dimj) + D{mo) . (6.55) 



Since the term F\ contains the finite cutoff-dependent term E (which is defined in 
Eq. (6.42)) it has a small dependence on the cutoff as well. The function Fi was 
calculated numerically by rewriting the terms involving D as an integral like it was 
done in Sec. 4.4. Then it is possible to subtract the integrands, instead of the large 
values of the integral. In this way, it is easier to avoid large numerical errors. 
The function F2 is defined by 



-i£i„g 



2 2a-2 

U{P, mT) + — + 



h P^ + 



m; 



(6.56) 



In order to calculate the function F2 a modified Abel-Plana formula was used, see 
Sec. (2.5). It turns out that due the suppression of high momentum modes at finite 
temperature, low momentum modes gives the main contribution to F2. This shows that 
the "high-energy approximation" of Bochkarev and Kapusta (1996) (implicitly defined 
as such a difference and directly related to f 2) is invalid, since in that approximation 
the high-momentum modes are assumed to give the main contribution. 

After renormalization, it turns out that the next-to-leading order contribution to the 
pressure is 



m; 



8;r2 



^^ 



^NLo-T^logFr -^T^logFr 



mt 



m7, 



^ 



2\ 



8;r2 



m: 



+ F1+F2 , 



(6.57) 



()/ 



which is shown in Fig. 6.4. At T = 0, fm^olT^ can be calculated exactly. This can 
be used as a check of the numerical calculations. At T = 0, clearly F\ - 0, and hence 
Pj^Lo/T'^ - F2/T'^. At r - 0, for low P, I{P,m) is dominated by 327rW(^^ + m^). 
This gives 



F2 



^^log{P^+m^) 



(6.58) 



For H ^0, the mass m vanishes and so 'Pnlo/T'^ ^ -;t^/90. For H ^ (104 MeV)^ it 
follows that Pt^ho/T'^ = 0, since T is in that case much smaller than m, such that the 
pressure is exponentially suppressed. 

The pressure ior H = is approaching the H i^ pressure at high temperatures, 
indicating that the effects of the explicit symmetry -breaking terms become smaller at 
higher temperatures. This is because // is a temperature-independent constant. 
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Figure 6.4: Next-to-leading order contribution to the pressure normalized to T^, as 
function of temperature ior H - Q and H = (104 MeV)^, for different values of the 
cutoff A. 



Pressure of the 0(4) linear sigma model 

In order to make contact with two-flavor low-energy QCD the 0(4) model is studied 
in this subsection. In Fig. 6.5, the pressure normalized by T"^ ior N = 4 and H = 
to next-to-leading order as function of T is shown. The LO pressure below T^ equals 
the pressure of a gas of four massless non-interacting scalars. This follows imme- 
diately from Eqs. (6.24) and (6.20). At NLO the sigma field becomes massive. For 
temperatures much lower than nia-, the contribution to the pressure from the sigma is 
Boltzmann suppressed and (to good approximation) it holds that P = 7t^T^/30, which 
is the pressure of a gas of three massless non-interacting scalars. From the calcula- 
tions presented here it can be concluded that the transition to next-to-leading order is 
of second (or higher) order since the derivative of the pressure is not diverging. 

In Fig. 6.6, the pressure normalized by T^ ior N - 4 and H = (104 MeV) to 
next-to-leading order as function of T is shown. 

In Figs. 6.5 and 6.6 the cutoff A = 5.0 GeV. To make contact with low-energy 
QCD, a few comments on this choice are in order. For the low-energy chiral La- 
grangian, the cutoff is usually taken to be Snf^ (using our definition of fj^), which is 
around 1.2 GeV. However, for the present purpose this value would be at the limit 
of applicability, since the critical temperature at which chiral symmetry is (approxi- 
mately) restored is only about a factor of 8 smaller and the requirement that 2nT «: A 
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Figure 6.5: LO and NLO pressure normalized to T^, for A'^ = 4 as a function of 
temperature, for // = and A = 5.0 GeV. 



should be satisfied. In this way one ensures that one sums over sufficient Matsubara 
modes. Therefore, the cutoff is taken to be considerably larger to reduce the sensitivity 
of the results to the cutoff. However, only for temperatures considerably lower than T^ 
one expects that the result is of actual relevance for the QCD pressure. 



Pressure of the 0(4) nonlinear sigma model 

In the limit X\, = oo, the Lagrangian of the nonlinear sigma model is obtained. In 
the nonlinear sigma model there are no counterterms available for logarithmic diver- 
gences. Therefore only f„^\, will be renormalized as in Eq. (6.43) with X\, = oo. This 
implies that F2 has a weak cutoff dependence. 

In Fig. 6.7 the pressure of the 0(4) nonlinear sigma model without explicit sym- 
metry breaking (H = 0), through next-to-leading order in l/N is shown. The pressure 
has been calculated for different values of the cutoff. The LO result for A = 20 GeV 
is included. For comparison, also the pressure resulting from the "high-energy ap- 
proximation" employed by Bochkarev and Kapusta (1996) is shown. A considerable 
difference between our results and those of the "high-energy approximation" is ob- 
served. 

In the approximations made by Bochkarev and Kapusta (1996), the T = part of 
sum-integrals are omitted such that every f is replaced by ^. Hence Il{P, m) is sim- 
ply replaced by IVjiP, ni). In the high-energy approximation the latter is furthermore 
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approximated by 



Ut{P, m) 



T^Jiipm) 



(6.59) 



32;r2A2 ' 

where /i and A^ are given by Eqs. (6.29) and (6.33), respectively. The term involv- 
ing (T^ in Eq. (6.12) is omitted because Bochkarev and Kapusta (1996) incorrectly 
assumed that this term is IjN suppressed with respect to the other contributions in 
Eq. (6.12). As a result the pressure in the "high-energy approximation" reduces to 



r 



Nm' 



.(/^<^)-f|;iog(/.^+M=)-i£iog— 



4m'^pl 



Defining the functions 



32 r°° p^ 

3r'* Jo c^p 



(6.60) 



(6.61) 



where u}± = yjp^ + nfi ± m, the pressure becomes 

'^4 i,jQ^\ „27^4 'Y'X 



p = ^(f;-^') 



NT'^Joifim) n^T 



6An^ 



90 64;r2 



[K^^ipm) + K-iPm)] . (6.62) 



Expanding Eq. (6.62) in powers of m/T and rescaling with factors of A'^, Eq. (52) of 
Bochkarev and Kapusta (1996) is obtained. 
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For completeness, the gap equations in this approximation are given by 

\(iT?fl = T^Jiifim) + 16n^d-^ - — [^^(jSm) + K:[(/3m)\ , (6.63) 



m cr 



0, 



where the functions Kf are 



Kf(J3m) = ±-;t dp ^ ^ n{co^) . 

T^ Jo (iJp m 



(6.64) 



(6.65) 



There are several problems with the approach of Bochkarev and Kapusta (1996). 
Firstly, it is incorrect to ignore zero-temperature contributions to the pressure and it 
also obscures renormalization issues. Then, as was argued below Eq. (6.56), the ar- 
guments for applying the high-energy approximation are not valid. Furthermore, one 
cannot neglect the term proportional to o" in I{P, m) for T < Tc- Fourthly, as the so- 
lutions to the gap equation (6.63) indicate, m/T becomes significantly larger than one 
for T > Tc, hence the m/T expansion breaks down. If one were to use Eq. (6.62) 
instead, one finds that the pressure even becomes negative above T k 300 MeV. An- 
other problem is that for T <Tc their pressure is equal to that of a massless gas, which 
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is incorrect since the sigma meson is massive and included at NLO. Hence, one ex- 
pects a deviation from the ideal-gas pressure niT <Tc. Finally, at high temperatures 
the next-to-leading order pressure should become approximately equal to the leading 
order pressure because chiral symmetry will be restored. This is not the case for the 
pressure calculated by Bochkarev and Kapusta (1996) as can be seen from Fig. 6.7. 

Jain (1993) has calculated the thermodynamic potential to NLO in the 0(A'^) linear 
sigma model using a high-temperature expansion. This approximation breaks down 
at low temperatures. Therefore it is not useful to compare those results to the results 
obtained in this chapter since they are strictly speaking only valid at low temperatures 
(much lower than the cutoff A). 



6.5 Choice of parameters 

The plots in the preceding sections are made using particular choices of the parameters, 
namely, /!(// = 100 MeV) = 30, f„ = 47 MeV (note that this f„ differs from the more 
conventional definition by a factor of 1/2) and if there is explicit symmetry breaking, 
H = (104 MeV)^. In this section these choices are motivated. For simplicity partly 
leading-order calculations are used for fixing the parameters. 

The values for fj^ and nij^ are chosen to be roughly equal to their measured values: 
f„ = 47 MeV and m,r - 138 MeV (the average of the measured masses of the ;r",7r^ 
and 7T~). These values are used for determining the parameter H as follows. Given 
a choice of A at some scale jj the LO renormalized gap equations (6.26) and (6.27) 
are solved for & and m^, such that m^ = m^ (which is the correct identification at 
leading order, see Eq. (6.16)). For the choice of /l(// = 100 MeV) = 30, this results in 
H = (104 MeV)l 

The choice of A is motivated by considerations on the maximal value of the cutoff 
and the sigma mass. As explained below, the sigma mass turns out to be maximal 
if /l(ju = 100 MeV) = 80. The problem with this choice of A is that in that case 
the maximal value of the cutoff is 720 MeV. This is very low and allows us only 
to do calculations up to around T = 50 MeV. Therefore, a lower value is chosen: 
A(ju = 100 MeV) = 30. Using that parameter choice A^ax = 19 GeV and the sigma 
mass is equal to 256 MeV and 350 MeV in the case oi H = md H = (104 MeV)^ 
respectively. 

To obtain the mass of the sigma field, one has to find the poles of the propagator 
in Minkowski space. The physical mass nipb is often defined by the solution to the 
equation 

-m , +m + ReZ{pQ - /mph + e, p - 0, m) - , (6.66) 

where S is the self-energy. Using the expression for the inverse o- propagator in the 
1/A^ expansion, Eq. (6.16), and choosing fi = nia- (because all results are independent 
of jU, one can choose fi as one likes, the choice fi = nia- is just convenient), one finds 
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that at r ^ and for // = 

"^ = 1 , 32.^ , .^ • ^6.67) 

^ "•" /l(m^) "'' l+32ff2/i(nv) 

Equation (6.67) can be maximized with respect to /l(mo-). This implies that mo- < 
Vl67r/;r w 333 MeV, which is lower than the averaged measured value of 400 - 800 
MeV. An analogous bound was found using the laige-N expansion by Einhorn (1984) 
on the mass of the Higgs boson. Patkos et al. (2002) show a figure in which the bound 
on the sigma mass can be seen indirectly, although they do not comment on this rather 
interesting fact. 

For H i^ a. similar bound applies. In that case the maximal value of the sigma 
mass can be found by maximizing 



32;r2/.'+2m2log(g) + 



647T^n 



ml = ml + Re , ^ ' , (6.68) 



">ste)+>+V 



AM , 1 , lml~4ml 



iogi— fPii+/^ 






with respect to A},. Finding the maximum of the previous equation requires solving the 
following equation for nzo- 



ml = 



V^+A^O 



l+^l+A\milmi) 



ml , (6.69) 



where 



A(x) ^ r^ - - , - - log , . (6.70) 

' ml n) Vl-4/x n * \ 1 + yi - 4/x/ 



By solving this equation it turns out that the maximal value of m^ is equal to 433 MeV. 
This is on the low side of the experimental measured values between 400 and 800 MeV. 
The reason that a rather low bound (which might be unrealistic) on the sigma meson 
mass is found could be because possible essential three-flavor physics was missed out 
in the calculation. 

The bounds in the previous section were obtained using non-perturbative large-A'^ 
expansion. So these bounds are valid for any value of A. In perturbation theory in 
small A, the mass of the sigma field is m^ = m^ + Ahcr^. This result can be found by 
expanding the expression for the inverse cr propagator, Eq. (6.16). One might think 
that this indicates that the sigma mass can grow to any value by increasing A. However, 
by increasing A, the sigma field will be coupled stronger to the pions. Since the sigma 
field is not stable, but can decay into two pions, the decay width is increased when A 
becomes bigger. As a result this decay shifts the mass of the sigma meson to lower 
values. 
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6.6 Summary and Conclusions 

In this chapter, the thermodynamics of the O(A^) linear and nonlinear sigma models to 
next-to-leading order in the 1 IN expansion was studied. 

At next-to-leading order it was shown that one can renormalize the effective poten- 
tial in a temperature-independent manner only at the minimum of the effective poten- 
tial. By renormalizing the next-to-leading order effective potential in the minimum the 
beta function for A to next-to-leading order was found. This beta function is consistent 
with the perturbative result. 

The pressure for the linear and nonlinear sigma model to next-to-leading order as 
a function of temperature was calculated numerically. The results show that for the 
calculation of the pressure IjN im good expansion, even ifN-4. With a relatively 
realistic choice of the parameters a prediction for the pressure of QCD for tempera- 
tures below Tc was made. The results for the pressure disagree significantly with the 
calculations of those by Bochkarev and Kapusta (1996). This is due to the fact that in 
the calculations performed in this chapter the zero-temperature contributions are not 
neglected and that the next-to-leading order contribution is treated without resorting to 
a "high-energy approximation". 

It was also found that in the linear sigma model the sigma mass has an upper 
bound at zero temperature. This bound depends only on the parameters f,^ and m„. For 
a realistic choice of these parameters, this implies that the mass of the sigma meson 
is smaller than 433 MeV. This does not necessarily have consequences for the real 
sigma meson, since the full three-flavor physics was not taken into account. It would 
be interesting to investigate the temperature dependence of this bound. 

An interesting extension of these calculations would be the the calculation of spec- 
tral functions at finite temperature. It would be interesting to see how the bound on 
the sigma meson mass depends on temperature. The methods developed in this chapter 
could also be useful for the study of more complicated models incorporating additional 
features of low-energy QCD. 



Chapter 7 

The phase diagram of the N JL model 



In order to obtain information about the structure of the QCD phase diagram, one can 
investigate the phase diagram of low-energy effective theories for QCD. In this chap- 
ter phase diagrams of such an effective theory, namely the Nambu-Jona-Lasinio (NIL) 
model, will be presented. These phase diagrams are obtained numerically as a func- 
tion of temperature and of the up, down and strange quark chemical potentials. Phases 
with broken chiral symmetry, color superconducting phases and phases in which the 
pseudoscalar mesons condense can be found in the calculated diagrams. It is shown 
numerically that color superconducting and pseudoscalar condensed phases are sepa- 
rated by a first order phase transition. This chapter is based on: Color superconduc- 
tivity versus pseudoscalar condensation in a three-flavor NJL model, H.J. Warringa, 
D. Boer and J.O. Andersen, Phys. Rev. D72 014015, (2005). 

7.1 Introduction 

The NJL model (Nambu and Jona-Lasinio, 1961), the instanton liquid model (Shuryak, 
1982) and random matrix models (see for example Verbaarschot and Wettig (2000)) 
are examples of effective models which can be used to study the QCD phase diagram 
at finite temperature and densities. Despite their shortcomings, it is expected that these 
models do describe the qualitative features of the QCD phase diagram in regions not 
accessible to perturbative or lattice QCD. For example, a critical point at finite baryon 
chemical potential has been predicted using effective models by Halasz et al. (1998). 
Furthermore, it is expected that due to the existence of an attractive interaction, quarks 
can form Cooper pairs. This is just as electrons do in ordinary superconductivity due to 
the attractive phonon interaction between electrons. As a result, at high baryon chem- 
ical potentials and low temperatures, quark matter can for instance be in a two-flavor 
color-superconducting phase (2SC) (Bailin and Love, 1984; Alford et al., 1998; Rapp 
et al, 1998) in which two flavors pair, or in a color-flavor locked (CFL) phase (Alford 
et al., I999a,b) in which three flavors pair. In these color-superconducting phases gaps 
in the quasi-particle excitation spectrum on the order of 100 MeV arise. These gaps 
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are large compared to the mass gaps of free quarks (the up and the down current quark 
masses are about 5 MeV), indicating that the physics in the color superconducting 
phase is very different from that in the deconfined phase. 

Many results on high density phase diagrams presented in the literature have been 
obtained for equal up, down and strange quark chemical potentials. However, this 
may not be directly relevant for heavy-ion collisions or compact stars. For example 
to model compact stars, one has to enforce electric and color neutrality and weak 
equilibrium. For this reason, different flavor and color chemical potentials have been 
introduced in the NJL model (Alford and Rajagopal, 2002; Steiner et at, 2002). This 
gives rise to a more complicated phase diagram in which one also finds new 2SC- 
like (Neumann et at, 2003) and gapless superconducting phases (Alford et al, 2000; 
Shovkovy and Huang, 2003; Alford et al, 2004; Riister et al, 2004; Abuki et al, 
2004; Riister et al, 2005; Blaschke et al, 2005). It is still an open issue whether these 
phases are really stable. Similarly, in heavy-ion collisions, a difference between the 
quark chemical potentials arises if the number densities of the different quark flavors 
are not the same. This difference can cause interesting observable effects such as two 
critical endpoints (Klein et al, 2003; Toublan and Kogut, 2003). However, instanton 
induced interactions tend to suppress this effect as was shown by Frank et al. (2003). 

In addition to a more complicated structure of the superconducting phases, differ- 
ent chemical potentials can also trigger pseudoscalar condensation (Son and Stephanov, 
2001 ; Kogut and Toublan, 2001). This has been confirmed on the lattice at zero baryon 
chemical potential (jub = fiu+l^d+l^s) but finite isospin chemical potential (juj = iXu-jJ-d) 
by Kogut and Sinclair (2002). In a pseudoscalar condensed phase, depending on the 
flavors involved, the charged pion, neutral or charged kaon field acquires a vacuum 
expectation value. As a result parity is broken spontaneously. Pseudoscalar condensa- 
tion in the two-flavor NJL model has been studied by Toublan and Kogut (2003) and 
Barducci et al. (2004), as a function of the different chemical potentials at zero and 
finite temperature. An extension to three flavors was carried out by Barducci et al. 
(2005) as wen. 

The phase diagram of the three-flavor NJL model as a function of the different 
chemical potentials including both pseudoscalar condensation in the quark-antiquark 
channel and color superconductivity had not yet been addressed. This is the sub- 
ject of this chapter. At zero temperature, pseudoscalar condensation is possible if 
l/^M - IJ-d\ > nijj, as was shown by Son and Stephanov (2001), or if |//„ (/ - //,! > mx 
(Kogut and Toublan, 2001). On the other hand, color superconducting phases occur if 
the chemical potentials are large and approximately equal. Therefore, one can imag- 
ine scenarios where for example //„ w p^i (the Fermi surfaces of the u and d quark 
should be sufficiently close for Cooper pairing to occur) and //„ w -//., (the Fermi 
surfaces of the u and s should be sufficiently close for kaon condensation to occur), 
with \^u4 - IJ-s\ > f^K. In such a case a 2SC phase is competing against a phase in 
which kaons condense. Hence the phase diagrams presented in this chapter are not 
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a superposition of the phase diagrams with only pseudoscalar condensation and with 
only color superconductivity, which were calculated before in the literature. From the 
calculations it follows that a coexistence phase of pseudoscalar condensation and color 
superconductivity does not occur for the parameters chosen and that these phases are 
separated by a first-order transition. However, it is not excluded that other choices of 
parameters may lead to such a coexistence phase, just as coexistence of color super- 
conductivity and chiral symmetry breaking may occur in the NJL model for specific 
ranges of parameters as was found by Blaschke et al. (2003). Here a coexistence phase 
is to be understood as a phase in which two condensates are nonzero simultaneously 
to be distinguished from a mixed phase in which phases coexist. 

In this chapter pseudoscalar condensation in the quark-antiquark channel is stud- 
ied. The pseudoscalar diquark interaction was not taken into account. This interaction 
is suppressed relative to the scalar diquark interaction due to instantons (Rapp et al., 
1998). However, in absence of instanton interactions, if one neutralizes the bulk matter 
with respect to color and electric charges it is possible to have pseudoscalar diquark 
condensation with rather large gaps as was shown by Buballa (2005b). According 
to Buballa (2005b) pseudoscalar diquark condensation in the NJL model is similar 
to pseudoscalar condensation in the CFL phase studied with eff"ective chiral models 
studied by Casalbuoni and Gatto (1999), Son and Stephanov (2000), Casalbuoni et al. 
(2000), Schafer (2000), Bedaque (2002), Kaplan and Reddy (2002) and Forbes (2004). 

Overall charge neutrality conditions are not applied in this chapter, this allows us to 
compare to the previous studies of Toublan and Kogut (2003), Barducci et al. (2004) 
and Barducci et al. (2005). Imposing neutrality conditions is necessary to describe 
situations in nature where the quark matter is realized at high densities such as in 
the core of neutron stars. Neutrality conditions would qualitatively affect the phase 
structure, leading for example to the observation that in a macroscopic volume of 
electric neutral quark matter in equilibrium with the weak interactions the 2SC phase 
is energetically disfavored (Alford and Rajagopal, 2002; Steiner et al., 2002). 

To obtain the phase diagrams an auxiliary field effective potential will be mini- 
mized with respect to the vacuum expectation values. This effective potential will be 
calculated in leading order in 1/A^c where Nc is the number of colors. For simplicity 
next-to-leading order 1 jNc corrections are not taken into account. However, it is pos- 
sible to calculate these corrections as was shown by Hiifner et al. (1994). The leading 
order 1 INc expansion amounts to the so-called mean field approximation. 

This chapter is organized as follows. In Sec. 2, it will be argued how the NJL 
model qualitatively arises from the QCD Lagrangian density. The NJL model itself 
and the choice of parameters will be treated in Sec. 3. In Sec. 4, the calculations and 
some of its technical aspects are discussed. In Sec. 5, the phase diagrams are presented 
and a summary and conclusions are given in Sec. 6. 
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7.2 From QCD to the NJL model 

The NJL model can serve as a low-energy effective theory of QCD. In this section 
it will be discussed qualitatively which approximations one has to make in order to 
obtain the NJL model from QCD. A related discussion about the connection between 
QCD and the NJL model can be found in Bijnens et al. (1993) and Bijnens (1996). 

The QCD Lagrangian density in Minkowki space, Eq. (1.1), can be written as 
Xqcd - Xf + Xym + gJfi^^a^ where the kinetic term of the fermions is given by 

Xf = 'A('/5^-Mo+//ro)^, (7.1) 

and the terms involving the gluon fields by 

Xym = —:Fa F% ■ (7-2) 

The interaction between the gluons and the quarks is written in terms of a fermion 
color current which is equal to 7fJ - if/T"y^if/, where T" is a generator of SU(3). The 
partition function of QCD can be factored in the following way 

Z - f D^DipDAl exp (/5qcd) - ( D^Dijj exp (iSf + iSiAJ^i) , (7.3) 

where the quark interaction action S'intL/] is given by 

S int [/] - -/ log r DA^ exp lis YM + ig f /^A^ j . (7.4) 

By adding a gauge fixing term and the corresponding ghost fields one can perform the 
integration over the gauge fields in such a way that the quark interaction action can be 
written as 

SintlJ] - 5i„,[7 = 0] + ^ f fr^{x)rfl{x - y)4(y) + ... . (7.5) 

where for simplicity the 6-point and higher fermion interaction terms will be neglected 
from now on (if g is small this is certainly allowed). In a more realistic low-energy 
effective theory the effect of the 6-point and higher fermion interactions should be 
taken into account. In Eq. (7.5) I^^ix - y) is the exact gluon propagator of pure glue 
QCD in a specific gauge 

lfl{x-y)^{A'',{x)Al(y)) . (7.6) 

In the limit of high momenta the pure gluon propagator becomes (in the Feynman 

(^ = 1) gauge) 

-if) j,?^^ 
KKp) = ^^ • (7-7) 
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At low momenta non-perturbative effects will modify this propagator. Since the exact 
gluon propagator cannot be calculated analytically in the low-momentum regime, one 
can try to make an ansatz for the exact gluon propagator. Such an ansatz should reduce 
to Eq. (7.7) in the high-momentum limit. One of the non-perturbative effects is that 
gluons can form massive bound states called glueballs (see for example Morningstar 
and Peardon (1999)). The effect of these glueballs should somehow be reflected in the 
gluon propagator. As an ansatz one could take 

where M^ mimics the effect of a massive glue-ball. This mass should be on the order 
of the scale at which non-perturbative physics starts to play a role, which is roughly 
1 GeV^. Comparing this ansatz to lattice calculations of the exact propagator for pure 
gluon QCD in the Landau gauge (^ - 0) (Bernard et at, 1994; Bonnet et ah, 2000; 
Silva and Oliveira, 2004) shows it to be fairly realistic. 

Since the objective is to obtain a low-energy effective theory for QCD, the ansatz 
can at low-momenta be crudely approximated by 

^at^P^-[Q p^>A\ ^'-^^ 

where A is an ultraviolet momentum cutoff. This ansatz, Eq. (7.9) will eventually lead 
to the NJL model. The term 1/M^ acts as a sort of coupling constant and should be in 
the order of 1 GeV~^. Because one has to fix a gauge to obtain the gluon propagator, 
one should realize that after a gauge transformation in principle the ansatz for the 
propagator should transform as well. Since this is usually not done in the NJL model, 
gauge invariance is lost, but a global SU(3)c symmetry is retained. 

To obtain a more realistic low-energy effective theory for QCD one should use 
a more general ansatz for the exact gluon propagator. Such a better ansatz could be 
found by fitting it to lattice calculations of the gluon propagator which were performed 
by Bernard et al. (1994), Bonnet et al. (2000) and Silva and Oliveira (2004). However, 
in most cases this will lead to non-local interactions between the quarks which will be 
more difficult to solve analytically in the end, see for example Hashimoto et al. (2005). 
In a first approximation the local interaction could be sufficient. 

Now at finite temperatures and densities, the term ^intL/ - 0] is the contribution of 
the pure gluon part of QCD to the total pressure of QCD. This term does not depend on 
chemical potentials but only on temperature. If one investigates phase diagrams as a 
function of chemical potential at a fixed temperature this term is just a constant added 
to the effective potential. Hence to first approximation the effects of 5int[/ - 0] can be 
neglected. At finite temperature, the exact gluon propagator depends on temperature, 
this could be taken into account in a more realistic ansatz. 
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As a result the QCD partition function can be approximated by 

Z^C f TliffDilj exp (/ r Xnjl J , (7.10) 

where the Lagrangian density of the NJL model is given by 

£nJL = £i - gJp^a , (V.ll) 

here g - g^/2M^ is a coupling constant. The term C contains the gluonic energy 
density and does not depend on chemical potentials but only on temperature. 

The qualitative arguments to obtain the NJL model will not play a role in the fol- 
lowing sections. Below the NJL model will be taken as a starting point. One must 
therefore be careful to extend the conclusions obtained below to QCD without cau- 
tionary remarks, especially at chemical potentials and temperatures close to the cutoff. 
Nevertheless, there are certainly reasons to believe that the qualitative features of the 
NJL model also extend to QCD. For instance the meson spectrum is well described by 
the NJL model, see for example Klevansky (1992). In the following section the NJL 
model will be discussed in more detail. 

7.3 The NJL model 

The NJL model (Nambu and Jona-Lasinio, 1961) was originally intended to describe 
the interactions between nucleons. Nowadays, however, it is mainly used as an effec- 
tive model for QCD. In QCD the interaction between the quarks is caused by gluon 
exchange. In the NJL model this gluon exchange is approximated by a four-point 
color current-current interaction. This is like the four-Fermi theory of the weak inter- 
actions in which the low-energy exchange of a massive vector boson is modeled by a 
four-point interaction. At energies much lower than the mass of the vector bosons, the 
weak interactions are almost independent of momentum. Therefore in that case the 
four-Fermi theory works well. In the previous chapter it was argued that by making 
an ansatz for the gluon propagator in the non-perturbative regime a fermion four-point 
interaction arises too. 

There are many NJL-like models. In this chapter a three-flavor NJL model with 
three colors is discussed. The different NJL models are based on a Lagrangian density 
which contains a kinetic term and a four-point color current-current interaction 

X = «A {iy^d^ - Mo + fiyo) ^-g {^yTaifff . (7. 12) 

The color, flavor and Dirac indices of the fermion fields tfr are suppressed in this equa- 
tion for notational simplicity. The diagonal mass matrix Mq contains the bare quark 
masses mo„, moj and mo.v. The matrix fi is also diagonal and contains the quark chem- 
ical potentials //„, jj^ and fig. The matrices Ta are the 8 generators of SU(3) and act in 
color space. They are normalized as Tr TaTh = 26 ab- 
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By applying several Fierz transformations to the current-current interaction (see 
for example Buballa (2005a)) and including only terms which give rise to attractive qq 
and qq channels, one obtains the following Lagrangian density 

X = ^ {iy^d, - Mo + //ro) ^ + Z-qq+ Z\q + £^q, , (7.13) 

where the quark-antiquark, the scalar diquark and the pseudoscalar diquark interaction 
term are respectively given by 

£g, - G [(^i,^)' + (^^/rs-A)'] , (V.14) 

-Cjq = jG{ij/tA^BCiy54f^)[V tA^BCiysilj) , (V.15) 

£^1 - ^G(^fA^BC^T(«A^fA/lBC^) , (7.16) 

where A,Be |2, 5, 7), since only the interaction in the antisymmetric color and flavor 
triplet channel is attractive. The matrices Aa are the 9 generators of U(3) and act 
in flavor space. They are normalized as Tr AaAh = 26 at- The matrices ta are the 
generators of U(3) and act in color space. Their normalization is Tr tJi, = 25 ab- To 
remind the reader, the antisymmetric flavor matrices A2, As and A-] couple up to down, 
up to strange and down to strange quarks, respectively. The charge conjugate of a 
field ip is denoted by i/'c - CiJ/^ where C = iyoyi- The coupling strength 3G/4 of the 
diquark interaction is fixed by the Fierz transformation (see Buballa (2005a)) and will 
be used in this chapter. However, some authors discuss the NJL model with a different 
diquark coupling constant (see for example Riister et al. (2005) for a comparison). 

The NJL model has a symmetry structure similar to QCD, except for the color 
symmetry. The NJL model is only invariant under global SU(3) color transformations 
since it does not contain gauge fields. In absence of quark masses and chemical poten- 
tials, the Lagrangian density of the NJL model has a global SU(3)c x U(3)£ x U(3)s 
symmetry. Due to the non-vanishing quark masses, the symmetry is like in QCD 
broken down to SU(3)c x U(3)y. Since the masses of the quark and the chemical po- 
tentials will be different, the symmetry of the Lagrangian density is further reduced to 
SU(3)c X U(l)g X U(l)/ X U(l)y, where B, I, Y stand respectively for baryon, isospin 
(the z-component) and hypercharge number. An additional U(1)a symmetry break- 
ing term is necessary to explain the large mass of the 77 particle. As was shown by 
't Hooft (1976) this is due to the non-trivial topological structure of the vacuum, the 
instantons. The effect of these instantons can be modeled by the following effective 
six-point interaction 

Xinst - K [det <A(1 - r5)«A + det <A(1 + ys)^] , (7. 17) 

where the determinant acts only in flavor space. In the two-flavor case one can argue 
by applying a Fierz transformation on this interaction, that the scalar diquark interac- 
tion is more attractive than the pseudoscalar diquark interaction (Rapp et al., 1998). 
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Therefore condensation in the scalar diquark channel is favored over pseudoscalar di- 
quark condensation. Likewise this is also expected to hold in the three-flavor case. 
Therefore from now on the pseudoscalar diquark interaction will be neglected. 

Unfortunately the six-point instanton-induced interaction gives a very complicated 
contribution to the eff"ective potential when allowing for pseudoscalar quark-antiquark 
condensates in the mean field approximation or beyond (Osipov et at, 2005). There- 
fore in this chapter this instanton-induced interaction is left out but it would be inter- 
esting to take this term into account in a further study. 

The results that will be presented in this chapter are obtained with the following 
choice of parameters which were also used by Buballa (2005a) 

mou - 5.5 MeV , mo^ = 5.5 MeV , mo, = 112 MeV , 
G-2.319/A2, K = 0, A = 602.3 MeV. *• ^^ 

These are the precise values used in the calculations, but clearly not all digits are 
significant implying that a small change of parameters will not have a big influence 
on the results. This choice of parameters gives rise to constituent quark masses M„ = 
Md = 368 MeV and M., = 550 MeV. In that case the pion and kaon have a mass of 
respectively 138 MeV and 450 MeV. 

7.4 Effective potential 

To obtain a phase diagram one has to investigate the behavior of order parameters. A 
good order parameter vanishes in one phase and is non- vanishing in another and related 
to symmetry breaking and restoration. In massless QCD good order parameters indi- 
cating chiral symmetry breaking are the quark-antiquark condensates (uu), (dd) and 
(ss). Since quarks in reality have a mass, chiral symmetry is already explicitly broken 
in the QCD and NJL Lagrangian. Even if one uses massless quarks, different quark 
chemical potentials also break the chiral symmetry. Therefore the quark-antiquark 
condensates will be nonzero in every phase. So the quark-antiquark condensates are 
no good order parameters in a strict sense, but since their values in the broken phase are 
much larger than in the approximately restored phase, one can use these condensates 
to distinguish different phases. So the quark-antiquark condensates are "approximate 
order parameters". The quark-antiquark condensates are not the only possible order 
parameters, others follow naturally from the discussion below in which it is explained 
how to calculate the condensates. 

As is discussed in Chapter 3, the quantity from which one can derive thermody- 
namical quantities is the effective potential. One also can find the values of certain 
condensates from this effective potential and hence determine the phase diagram. To 
calculate the effective potential of the NJL model it appears to be useful to introduce 1 8 
auxiliary real scalar fields aa and yS^ and 9 complex scalar fields Aab in the following 
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way 

r ^ r - - (1 1Q^ 

■^ -^ AG 3G ■ ^ ' 

By shifting these auxiliary fields as follows 

3 J 
cua^ a:a + IGif/Aai// , Aab -^ ^AB " tGiA tAABCysi// , (7.20) 

fia^l3a + 2GiJ^AJy5iff , A\g ^ A^^ + -G«A?A^BCr5«A^ , (7.21) 

the Lagrangian density Eqs. (7.1), (7.14) and (7.15) becomes quadratic in the fermion 
fields 



X = (A [iydft - Mo - aaAa - iysPa^a + fiyo) (A ■ 



AG 

^AABiptAABCyiip'^ + \A\j^iI/ tAAsCynp -— ^ 



(7.22) 



which allows integration over the quark fields. These quark fields \fj have 3 (color) 
x3 (flavor) x4 (Dirac) - 36 components. As is explained in appendix 7.A the in- 
tegration is most easily done by introducing a two-component Nambu-Gorkov field 
^^ - up^, ijA I Vz which has 72 components in this case. After integration over the 
fermion fields one obtains an effective action which depends on the auxiliary fields. 
To obtain the effective potential to leading order in 1/A^c, these auxiliary fields are re- 
placed by their vacuum expectation values. The leading order 1 INc approximation is 
equivalent to the so-called mean field approximation. From now on a a, Pa and Aab 
stand for the vacuum expectation values of the auxiliary fields. The quantum fluc- 
tuations around the vacuum expectation values of the auxiliary fields are not taken 
into account in the mean field approximation. The vacuum expectation values of the 
auxiliary fields are directly related to the quark condensates 

aa = -2G {ipAaip) , Pa = -2G (ipAaiysi/y) , 
Aab ^ tG ^(A tAABCysi//) . 

The condensate Aab is complex. To calculate its vacuum expectation value one should 
in principle average over all its phases. Because of the U(1)b symmetry of the action, 
this averaging will let the vacuum expectation vanish. However, in a real physical sit- 
uation the system will just pick one direction giving rise to a nonzero vacuum expecta- 
tion value of Aab- The phase will however be undetermined unless a small asymmetry 
is brought in. This is situation is similar to the magnetization below the phase transi- 
tion temperature in the Ising model. In that case averaging over all possibilities will 
also give rise to a zero net magnetization. However the real system should choose to 
have either a finite positive or negative magnetization. It turns out that Aab is a good 
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order parameter, giving rise to a BCS gap in the quasi-particle excitation spectrum 
indicating color superconductivity. 

In this chapter it is assumed that all condensates are space-time independent. When 
this restriction is dropped it is possible to find the crystalline Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) phase (Alford et al. , 2001). This is a phase in which quarks with 
different momentum magnitudes pair. One can find this phase between the phase in 
which chiral symmetry is broken and the color superconducting phases. 

After going to imaginary time, the thermal effective potential "V to leading order 
in 1 /Nc (which is equivalent to the mean-field approximation) reads 

T^^i^.^^-I y fi^^logdeti., (7.24) 

where ^ is a 72 x 72 matrix 

and 

Di - 1/ ® (iyoPo+yiPi) - yu ® To - (Mo+aaAa) ® Id -pa^a ® iys , (7.26) 
D2 - 1/ ® {ijopo +yiPi) +^i®yo- (Mq + aa^l) » 1^ - jia^l ® iys . (7.27) 

The matrix 1 is the identity matrix in color (c), flavor (/), or Dirac {d) space. 

The values of the condensates and the phase diagram are determined by mini- 
mizing the effective potential 'V with respect to the condensates. To make the min- 
imization procedure easier, one can take advantage of the fact that certain conden- 
sates must vanish. Firstly, application of QCD inequalities derived by Weingarten 
(1983) shows that in QCD at zero temperature, global SU(3)y symmetry breaking 
cannot be driven by condensates of the type {uiy=,u) (Vafa and Witten, 1984b; Son 
and Stephanov, 2001). Therefore, Po, P3 and/Jg are zero. Outside the phase in which 
diquarks condense, it was checked numerically that this is correct at non-zero temper- 
atures and chemical potentials as well. At zero chemical potential and temperatures a 
theorem by Vafa and Witten (f 984a) states that all parity violating condensates should 
vanish in QCD, this is also in agreement with the calculations. Moreover, although 
perturbative one-gluon exchange cannot distinguish between p^ and a^ condensation 
with k € {1,2,4,5,6,7}, pseudoscalar condensation is favored due to the instanton 
interaction (Son and Stephanov, 200f). One can therefore set all a^s again with 
k e {1,2,4,5,6,7} to zero. Numerically it was found that this is correct, even in 
the absence of instanton interactions. 

One can further simplify the minimization procedure by exploiting the symme- 
tries of the NJL model. The free energy (which is the energy in the minimum of the 
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effective potential) is invariant under the same transformations as the Lagrangian den- 
sity, so applying a U(l)-flavor transformation to all condensates leaves the free energy 
invariant. The (3 condensates transform under a U(l)-flavor transformation which is 
parametrized by an angle a as follows 



/3l,4,6 -^ COStKi,4,6/3l,4,6 + sin CC 1,4,6 y62,5, 7 , (7.28) 

/32,5,7 -^ - sin Q'l,4,6jSl,4,6 + COS Q'l,4,6jS2,5,7 ■ (7.29) 



Therefore, using the U(l)-flavor transformations one can choose the pseudoscalars to 
condense in the P2, Ps, and yS? channels, and set fix, JS4 and ySg to zero. The phase in 
which P2, 1^5 and/or yS? is non- vanishing is called the ii^ jn~ , K^ jK^ and/or K^ IK^ 
condensed phase, respectively. The condensates ^2, Ps andyS? are good order param- 
eters since they break the U(l)j, x U(l)j X 11(1)^ symmetry which is present in the 
Lagrangian density into U(1)b x U(1) if only one pseudoscalar condensate arises or 
to U(1)b if there are more. Moreover, these condensates break parity symmetry. To 
summarize for all quark-antiquark condensates, only the aQ, ai„ o-g and/32»/35>/57 con- 
densates will be considered below. 

Because of the global SU(3)c symmetry, one can also rotate away several diquark 
condensates. Without loss of generality one can minimize with respect to A22, A25, 
A55, A27, A57 and A77. In principle, all six diquark condensates can have a phase. It is 
always possible to remove two of them by using the two diagonal SU(3)c transforma- 
tions. As long as there is no pseudoscalar condensation, one can use the U(l)-flavor 
symmetries to rotate away three other phases. As a result either A25, A55, A27, or A57 
has a phase (Buballa, 2005a). However, this reduction is not completely possible if 
pseudoscalar condensation occurs. By choosing the pseudoscalars to condense in the 
P2, f^5 and/37 channels, one breaks the U(l)-flavor symmetry. Hence if pseudoscalar 
condensation arises in one channel, one can in general rotate away one phase less in 
the diquark sector. If it occurs in more channels, two phases less can be rotated away. 
However, numerically it is found that allowing for such a complex phase leads to di- 
quark condensation only in the A22, A55 and the A77 channels. The A25, A27, and A57 
diquark condensates do not arise or can be rotated away. Moreover it is found that 
pseudoscalar condensation in the quark-antiquark channel does not coexist with color 
superconductivity, such that one can always take the diquark condensates to be real. 

The different possible color-superconducting phases are named as follows (Riister 
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etal, 2005) 

A22 ^ 0, A55 ?t 0, A77 i= CFL , 

A77 - 0, A22 i= 0, A55 7t uSC , 

A55 - 0, A22 i= 0, A77 i^O dSC , 

A22 - 0, A55 7^ 0, A77 7^ sSC , 

A22^0, A55=0, A77=0 2SC, 

A55 * 0, A22 = 0, A77 = 2SCus , 

A77 7t 0, A22 = 0, A55 = 2SCds . (7.30) 

The abbreviation CFL stands for color-flavor locked phase. If there is exact SU(3)v 
flavor symmetry, the three diquark condensates in this phase have equal size and the 
vacuum is invariant under a combined rotation in color and flavor space (Alford et al. , 
1999a,b). In the uSC, dSC or sSC phase the up, down or strange quark always takes 
part in the diquark condensate, respectively. In the 2SC phase an up and a down 
quark form a diquark condensate, in the 2SCus and the 2SCds phase this condensate is 
formed by an up and and strange quark and a down and a strange quark, respectively. 

To calculate the effective potential one needs to evaluate a determinant of a 72 x 72 
matrix. Only in special cases such as when all masses and chemical potentials are 
equal and in absence of pseudoscalar condensation one can perform over the sum over 
Matsubara frequencies analytically and hence simplify the effective potential some- 
what further. But in the more general cases which will be discussed in this chapter 
this either is not possible or gives rise to very complicated equations. One is therefore 
restricted to do numerical calculations. 

To calculate the determinant of the effective potential in an efficient way, one can 
multiply the matrix K with diag(lc® l/^To^ Ic® l/®yo) which leaves the determinant 
invariant. In this way, one obtains a new matrix K' with ipo's on the diagonal. By 
determining the eigenvalues Aj of the matrix K' with po = 0, one can reconstruct the 
determinant of K for all values of po which is namely Y[]=ii'^i + ^Po)- After summing 
over Matsubara frequencies, one finds using Eq. (2.55) 



72 , 

T Y, log det^ - Z [ Y + ^ ^°§ (1 + '''''''') 

PQ=(2n+l)nT i=l 



(7.31) 



All that remains in order to determine the effective potential is to integrate over three- 
momentum p up to an ultraviolet cutoff A. 

The speed of the calculation of the effective potential depends heavily on how fast 
one can compute the eigenvalues. There are several ways to speed up the calculation. 
Firstly, the determinant of K does not depend on the direction of p. Therefore, one 
can choose p to lie in the z-direction. Together with the choice of the non- vanishing 
condensates mentioned above, this implies that K'(pQ = 0) becomes a real symmetric 
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matrix, which simpUfies the calculation of the eigenvalues. Secondly, one can inter- 
change rows and columns of K' without changing its determinant. By doing so, one 
can bring K' in a block-diagonal form. One can then determine the eigenvalues of the 
blocks separately which is significantly faster since the time needed to compute eigen- 
values numerically scales cubically with the dimension of the matrix. In the most 
general case with diquark condensation, one can always reduce the problem to two 
36 X 36 matrices. Moreover, if there is no diquark condensation, but only pseudoscalar 
condensation (Barducci et al. , 2005), the problem can be further reduced to computing 
the eigenvalues of two 6x6 matrices 

The eigenvalues were obtained using LAPACK routines (Anderson et al, 1999). 
After numerical integration over three-momentum p up to the cutoff, the condensates 
were determined by minimizing the eff"ective potential using MINUIT (James and 
Roos, 1975). To be certain that the minimization procedure did not end up in a local 
minimum the continuity of the minimized eff"ective potential as a function of chemical 
potentials and/or temperature was always checked. 

7.5 Phase diagrams 

In this section results for the phase diagrams of the NJL model with u, d, and s quarks 
are presented. The phase diagrams are plotted as a function of the chemical potentials 
and temperature. To determine the locations of the phase boundaries, the behavior of 
the condensates was examined. The different possibilities are displayed in Fig. 7.1. If a 
condensate jumps discontinuously the transition is first-order (a), and this is indicated 
by a solid line in all phase diagrams. If its derivative has a discontinuity, the transition 
is second order (b), and this is indicated by a dotted line. If a condensate changes 
rapidly in a narrow range without vanishing, in other words no derivative has a discon- 
tinuity, there is a smooth cross-over (c), and this is indicated in the phase diagrams by 
a dashed-dotted line at the point were the condensate varies maximally. 




Figure 7.1: Typical behavior of condensates near a transition. Fig. (a) shows a conden- 
sate undergoing a first order transition. Fig. (b) a second order transition and Fig. (c) a 
cross-over. 
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Figure 7.2: Baryon number density as a function of baryon chemical potential for 
r = and fJ.u = l^d = l^s- The vertical pieces indicate the presence of a mixed phase. 
Using Fig. 7.4 one can conclude that upon increasing ns, quark matter is subsequently 
in a phase with broken chiral symmetry, in a mixed phase of broken chiral symmetry 
and 2SC, in a 2SC phase, in a mixed phase of 2SC and CFL, and finally in a CFL 
phase. For comparison, the nuclear matter density is about 0.17 fm"^. 
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Figure 7.3: The chiral condensates (solid line) and pion condensate (dotted line) as a 
function of /// - fiu - f^d with //„ = -/z^. 
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To calculate the phase diagrams, the effective potential was minimized on a grid of 
chemical potentials and or temperatures. The numerical uncertainties in the presented 
phase diagrams are in the order of the distance between the points in the grid which 
was about 3 MeV, this coincides with the thickness of the lines in the phase diagrams. 

One should keep in mind that the relation between chemical potential and number 
density is not linear. In Fig. 7.2 the baryon number density as a function of the baryon 
chemical potential is displayed. In that figure it can be seen that at a first-order phase 
boundary the number density increases discontinuously. At these particular densities, 
quark matter for example can be in a mixed state of normal and superconducting matter 
(see e.g. Bedaque (2002) and Lawley et al. (2005)). 

A typical outcome of the minimization procedure is displayed in Fig. 7.3, in which 
the chiral and pion condensates are plotted as a function of/// for //^ - 0. It can be seen 
from the figure that pion condensation arises via a second order transition for jij > m,^. 
It can also be seen that the chiral condensates ilu and dd are small but non-vanishing 
in the phase with pion condensation. 

Phase diagram with //g vs. T 

In Fig. 7.4 the phase diagram of the NJL model as a function of baryon chemical 
potential jjb = f^u + l^d + l^s and temperature is displayed for fiu = t^d = l^s- This 
phase diagram is similar to the qualitative QCD phase diagram displayed in Fig. 1.3. 
At low temperatures and baryon chemical potential the chiral symmetry is broken 
in the phases denoted by (a) and (g). By increasing the chemical potential at low 
temperatures one enters the color superconducting phases, first the 2SC phase (n/q) and 
thereafter again via a first order transition the CFL phase (t). Three critical endpoints 
can be seen in this diagram. 
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Phase diagram with /Ug = and T -0 

In Fig. 7.5 the phase diagram of the NJL model for /iy = and T = as a function of 
jUu and jUd is displayed. Outside the 2SC phase (q), the results agree qualitatively with 
the two-flavor calculations by Barducci et al. (2004) (see their Fig. 1), where color- 
superconducting phases were not taken into account. Moreover, Barducci et al. (2004) 
use different parameters, in addition to a form factor. One can clearly see that the phase 
diagram is symmetric under reflection in the origin. This is because the free energy 
is invariant under the transformation (jUu,lUd,l^s) -^ {-jJ-u, -jJ-d, -j^s), that stems from 
the symmetry between particles and antiparticles. Fig. 7.5 is also symmetric under 
interchange of u and d, because of the choice of equal up and down quark masses. 
This gives rise to the symmetry of the phase boundaries with respect to the diagonals. 

In general, horizontal and vertical lines in the phase diagrams arise if the pairing 
of one type of quark is not changed after a transition. In this case, the location of the 
phase boundary is determined by the properties of other quarks. Therefore, changing 
the chemical potential of the unchanged quark species cannot have a big influence on 
the location of the phase boundary. This results in the horizontal and vertical lines. For 
r = 0, one always finds these lines near the values of the constituent quark masses, 
i.e. Hu ~ Mu, fid ~ Md and fi^ * Mg (see for example Buballa (2005a)). The diagonal 
lines arise because at T = pion condensation can occur if |//„ -/Udl > ruj, - 138 MeV 
(Son and Stephanov, 2001). 

The diagram shows that if the chemical potentials are different, the transition to 
the 2SC phase (q) remains first order as was concluded by Bedaque (2002). Moreover, 
one can see from Fig. 7.5 that if //„ + jud it is possible to go through two first-order 
transitions before entering the 2SC phase (q) (similar to the situation discussed by 
Toublan and Kogut (2003) without color superconductivity). To have such a scenario 
at zero temperature, a minimum difference between jU„ and fid is required. In the 
present case this is about 35 MeV. Pion condensation (i) and the 2SC phase (q) are 
in this diagram separated by two phase transitions in contrast to the estimated (jU^,///) 
phase diagram of He et al. (2005) which correctness is therefore questionable. 
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Phase diagram with i^d = and T = 

In Fig. 7.6 the phase diagram for /Ud = and T = as a function of //„ and /u^ is 
displayed. Since the up and down quark masses are much smaller than the strange 
quark mass, this diagram is very different from Fig. 7.5. Besides the possibility of 
pion condensation in (h) and (i) also phases in which the charged kaon (k) and the 
neutral kaon condense (l)/(m) arise. The lines separating the charged kaon phase (k) 
from the chirally broken phase (a) are diagonal because at T = kaon condensation 
can occur if I//., - fiu,d\ > niK - 450 MeV (Kogut and Toublan, 2001) (the chosen 
parameter set gives rise to a somewhat low kaon mass, but this is not relevant for 
the qualitative features of the phase diagram). The 2SCus phase appears in (r). This 
phase is surrounded by phases in which the pions (h)/(i) and the neutral kaons (1) 
condense. One passes a first-order transition when going from the pion and neutral 
kaon condensed to the 2SCus phase. 
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Figure 7.6: Phase diagram for fi^i = and T = as a function of ju„ and /i.,. First 
and second-order transitions are indicated by solid and dotted lines, respectively. The 
letters denote the different phases, where a: uu + dd + ss, b: uu + dd, d: dd + ss, h: 
n^/n-, i: n^/n- + ss, k: K^/K' + dd, 1: K^^/K^\ m: K^/K^^ + uu and r: 2SCus + dd. 
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Phase diagrams with //„ ~ ^i^ 

In Fig. 7.7 the phase diagram for T = as a function of the up and down quark 
chemical potential and the strange quark chemical potential is displayed. In this phase 
diagram ju„ = //^ + e where 6 is a very small positive number. This e is necessary 
because when e = one is just at a first-order phase boundary between the phase in 
which the charged kaons condense (k) and the one in which the neutral kaons con- 
dense (m), as can be seen from Fig. 7.6. This nonzero value of e gives rise to a small 
asymmetry in the phase diagram. If e is chosen to be negative, the phases in which the 
neutral (l)/(m) and the charged (j)/(k) kaon condenses are interchanged. 
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+ dd, 1: K*-^/K^\ m: is:"/^" + uu, n: 2SC, q: 2SC + ss and t: CFL. 



Apart from the additional 2SC (n)/(q) and CFL (t) phases, the results presented 
here agree qualitatively with the three-flavor calculations by Barducci et al. (2005) 
(see their Fig. 7). Barducci et al. (2005) have used different quark masses and a dif- 
ferent coupling constant, and in addition employed a form factor to mimic asymptotic 
freedom. Therefore, one may conclude that the use of such a form factor does not 
affect the phase diagram qualitatively. The phase diagram Fig. 7.7 cannot simply be 
obtained by a superposition of phase diagrams obtained from a calculation with pseu- 
doscalar condensation, but without superconductivity (such as done by Barducci et al. 
(2005)), and one with superconductivity, but without pseudoscalar condensation (such 
as done by Gastineau et al. (2002)). Despite the fact that the two types of phases do 
not coexist, there is nevertheless competition between them. Figure 7.7 shows that the 
K^IK^ (l)/(m) and the K^ jK' (j)/(k) phase are separated from the 2SC phase (q) by 
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a first-order transition. This remains the case at finite temperature as is illustrated in 
Fig. 7.8. This figure displays the phase diagram as a function of /i,, and temperature, 
for fixed fiu = Md = 550 MeV. 
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Figure 7.8: Phase diagram as a function of /i.^ and T, for fixed fiu = l^d = 550 MeV. 
Here a cross-overs is indicated by a dashed-dotted line. The letters denote the different 
phases, where g: ss, 1: K°/K°, n: 2SC, q: 2SC + ss and t: CFL. 



Returning to the discussion of Fig. 7.7; the line jJ-u - jJ^d - l^s goes through the 
phase (a) in which chiral symmetry is spontaneously broken. At some point it enters 
via a first-order transition the 2SC -i- 5^5' phase (q), and finally goes into the CFL phase 
(t), again via a first-order transition. If there is a difference between jU„ = fid and 
jUg, one can see in Fig. 7.7 that as the densities increase, quark matter can go directly 
from a phase with chiral symmetry breaking (a) to a CFL phase (t) without passing the 
2SC phase (q) first. This can also occur in compact stars as was shown by Alford and 
Rajagopal (2002). 

It is also interesting to note that the phases (l)/(j) of kaon condensation can also 
occur outside the region of spontaneous chiral symmetry breaking. Assuming the 
phase transition towards chiral symmetry restoration coincides with the deconfinement 
transition (as appears to be the case in lattice studies at small baryon chemical potential 
and in some models), this would imply that condensation of a state with quantum 
numbers of the kaon may persist in the deconfined phase. This was first observed by 
Son and Stephanov (2001), based on a perturbative calculation at very high isospin 
chemical potential. 
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Figure 7.9: Phase diagram for T = as a function of //„ = /i^ and ju^. The letters 
denote the different phases, where a: uu + dd + ss, c: uu + ss, d: dd + ss, g: ss, i: 
;r+/7r- -F ^.? , I: K'^/K'^, o: 2SCus, q: 2SC + ss, r: 2SCus -F dd and t: CFL. 



Phase diagrams with /i„ = /i , 

In Fig. 7.9 the phase diagram at zero temperature as a function of //„ = fig and ju^ is 
shown. This diagram is similar to Fig. 7.5 for small strange quark chemical potentials 
(below the kaon mass). At larger strange quark chemical potentials the diagrams differ, 
exhibiting kaon condensation (1) and diquark condensation involving strange quarks 
(r)/(o)/(t). 

In Fig. 7.10 the phase diagram as a function of ju^ and T, for fixed ju„ = fig = 
550 MeV is displayed. In this figure one can find five critical points. Also, one can 
see in this figure that the K^^K^^ phase (1) is separated from the 2SCus phase (r) by 
a first-order transition for all temperatures. Furthermore, it is interesting that at finite 
temperature there is a first-order transition from the phase in which the neutral kaons 
condense (1) to the pion condensed phase (h). In Fig. 7.11 the lower-right corner of 
Fig. 7. 10 is enlarged for clarity. In this figure one can find all possible superconducting 
phases, including the more exotic uSC (u), dSC (v) and sSC (w) phases. For //„ = fid = 
jUs one goes from the CFL phase (t) via the 2SC phase (n) to the chirally restored phase 
when raising the temperature. However, small differences between //„ = jUg and fid can 
cause one to go through completely different phases. 
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Figure 7.10: Phase diagram as a function of ju^ and T, for fixed //„ - iig - 550 MeV. 
The letters denote the different phases, where f: dd, h: n^ jn" , 1: K^ jK^, o: 2SCus, r: 
2SCus -I- dd, t: CFL and w: sSC. The upper phase is phase with chiral symmetry in 
which all condensates vanish. The lower right corner of this figure is enlarged in Fig. 
7.11. 
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Figure 7.11: Same as Fig. 7.10. The phases that occur are f: dd, n: 2SC, o: 2SCus, p: 
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7.6 Summary and Conclusions 

In this chapter the phase diagram of the three-flavor NJL model including pseudoscalar 
quark-antiquark condensation and color superconductivity was studied as a function 
of the different quark chemical potentials and temperature. The NJL model has a rich 
and interesting phase structure. The pseudoscalar condensed and color superconduct- 
ing phases are competing and are separated by a first-order phase transition. As was 
discussed, this need not be the case for other (less conventional) choices of the param- 
eters of the model. 

Furthermore, it was concluded that at zero temperature and zero strange quark 
chemical potential, a minimum asymmetry of about 35 MeV between the up and the 
down quark chemical potentials is required in order to have two first-order transitions, 
when going from the phase with spontaneous chiral symmetry breaking to the 2SC 
phase. 

The results provide a qualitative check and extension of several earlier calculations 
that appeared in the literature. The calculations in this chapter simultaneously allowed 
for nine different nonzero condensates, uu, dd, ss, n^ In^ , K^ /K^, K^^ IK^, A22, A55 
and A77. The new aspects of the phase diagrams are often located in regions, where 
the quark chemical potentials are large and very different in magnitude for the differ- 
ent flavors. Although such situations are not necessarily realized in compact stars or in 
heavy-ion collisions. The study of such unusual situations is nevertheless interesting 
for a fundamental understanding of the theory. Moreover a comparison with future 
lattice data may provide interesting information. This is especially relevant for pseu- 
doscalar condensation in the phase where chiral symmetry is restored and also for the 
complicated superconductivity phase structure close to the cutoff of the model. 

This work can be extended in several ways. For example, one can take into ac- 
count the instanton-induced interaction ('t Hooft, 1976), Eq. (7.17) in the mean field 
approximation or using the method of Osipov et al. (2005). If one has pseudoscalar 
condensation, this is much more difficult than in the normal case though. Another 
very useful extension would be the inclusion of the neutrality conditions (Alford and 
Rajagopal, 2002; Steiner et al., 2002), in which case the phase structure changes and 
for instance gapless phases will occur. It would also be interesting to see how the 
results depend on the strength of the diquark coupling and on the quark masses. Fur- 
thermore, one could add the LOFF phase (Alford et al., 2001). In this crystalline 
phase, quarks of different momentum magnitudes can pair. One could also include 
vector interactions. In this case spin-1 diquark condensation (see for example Pisarski 
and Rischke (2000); Buballa et al. (2003)) and an induced Lorentz-symmetry break- 
ing (ISB) phase (Langfeld and Rho, 1999) are among the possibilities. It would also 
be worthwhile to take pseudoscalar diquark condensation (Buballa, 2005b) into ac- 
count. Finally, one could try to go beyond the mean-field approximation as was done 
by Hufner et al. (1994). 
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7.A The Nambu-Gorkov formalism 

In case a Lagrangian density contains terms which are proportional to ip^Tip, the inte- 
gration over the quark fields is non-trivial. In that case it useful to apply the Nambu- 
Gorkov formalism which is discussed for the Lagrangian density used in this chapter, 
Eq. (7.22). 

Taking the fermionic part of the Lagrangian Eq. (7.22) and transposing half of the 
quark-antiquark term gives 






(7.32) 



where M = Mq - aa^-a - iysl^a^a- By introducing a two-component Nambu-Gorkov 
field 

one can rewrite the Lagrangian in Eq. (7.32) as Xf = ^^5^ with 

^ ^ / ^\stAABCys iCrCd^ + M^ - fijo \ .^ 34. 

\ iy^'df, - M + fiyo -AABtA^BCys j' 

In the path integral the measure 'Dif/'Dip will be changed into D*F. Integration over the 
Nambu-Gorkov fields yields that S ^s = ^logDet ( 5 -5^1/2 = ^logDet5 since 
S is antisymmetric in this case. The determinant is invariant under interchange of an 
even number of rows. Using that one can put the part containing the kinetic term on 
the diagonal. The eff'ective action is simplified further if one multiplies from the left 
with diag (1, C) and from the right with diag (1, C) where C is the charge conjugation 
matrix, resulting in 

5eff = logDet ( '•r'^^. - ^//^ro ^ABtAABy. \ 

\ -A^j^tA^-Bys ly'^df, -M' - idyo j 

This expression was evaluated further in Eq. (7.24), using a basis for a complete set of 
functions which satisfy the anti-periodicity conditions for fermions in imaginary time 
like in Eq. (2.43). 
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In this thesis, entitled "Thermodynamics of QCD-inspired theories" the behavior of 
(quark) matter under extreme high temperatures and densities is investigated. Such 
extreme conditions can be reached for example in the early universe shortly after the 
big bang, during a heavy-ion collision and inside the core of a neutron star. They can 
be so extreme that normal nuclear matter does longer exist, such that matter is in a 
different phase than usual, like for example the so-called quark-gluon plasma (at high 
temperatures and/or densities) or the color-superconducting phase (at high densities). 

A general introduction to this thesis is given in Chapter 1 . Firstly, it will be dis- 
cussed what could happen to matter under such extreme conditions. Then, the relevant 
aspects of quantum chromodynamics (QCD) (the theory which describes the interac- 
tions between the quarks) will be explained shortly. Thereafter, an overview of the 
results and problems of previous calculations of the pressure and the phase diagram of 
quark matter are given. The relevance of these kind of calculations for the understand- 
ing of a compact star, heavy ion collisions, and the state of the universe shortly after the 
big bang, will be examined as well. At the end of the chapter the term "QCD-inspired 
theories" will be explained. 

In Chapter 2 it is discussed how one can perform calculations at finite temperature 
and densities in a field theory like QCD. For that, results from statistical physics will be 
applied. Firstly, a short review of classical statistical physics is given, which thereafter 
will be generalized to a quantum mechanical approach. Using path-integrals it will be 
discussed how this approach can be used in order to investigate quantum field theories 
at finite temperature and densities. In calculations at finite temperature and densities 
one often has to perform summations over an infinite number of terms. A method to 
perform these sums exactly will be discussed. Also it will be explained in detail how 
these summations can be computed numerically for more difficult cases. 

The third chapter treats the so-called effective potential and the IjN expansion. 
The effective potential is an important quantity, since the pressure and the phase dia- 
grams can be obtained by minimizing this potential. The 1 jN expansion is a method 
which can be used to investigate the non-perturbative (in the coupling constant) behav- 
ior of some theories. It will be explained how this expansion can be obtained by the 
introduction of auxiliary fields. At the end of this chapter an important result of this 
thesis, which is the occurrence of temperature-dependent divergences in the effective 
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potential, will be discussed. Such divergences cannot be renormalized properly. These 
divergences become independent of the temperature at the minimum of the effective 
potential. That is just at the point at which one has to evaluate the effective potential in 
order to obtain the physical quantities like the pressure and the phase diagram. At the 
minimum, these divergences are no problem because there they can be renormalized 
in a systematic way. 

The subject of Chapter 4 is the thermodynamics of the nonlinear sigma model in 
two dimensions. This model has some aspects with QCD in common such as asymp- 
totic freedom. Using the 1 /N expansion it is possible to investigate this model in the 
domain where it is non-perturbative in the coupling constant. By studying this model 
one can learn something about QCD at finite temperature, for example how the behav- 
ior of an asymptotically free theory changes when going from low to high tempera- 
tures. In this chapter, the effective potential will be calculated to next-to-leading order 
(NLO) in 1 /N. It turns out that this effective potential contains temperature-dependent 
divergences which can only be renormalized properly at the minimum. By minimiz- 
ing the effective potential, the pressure as a function of temperature is obtained. It is 
found that 1 /A'^ is a good expansion method. Furthermore it turns out that the pressure 
divided by the pressure in the limit of infinite temperature is almost independent of A'^. 
Similar behavior is found in lattice calculations of SU(A^) Yang-Mills theory. 

Chapter 5 discusses the thermodynamics of the CP'^"^ model in two dimensions. 
This model is an extension of the nonlinear sigma model in which also electromag- 
netism is taken into account. This model shares also some features with QCD, it 
contains for example like QCD so-called topological non-trivial vacua. The influence 
of these vacua on the pressure is investigated using the same methods as in Chapter 4. 
It is found that the non-trivial vacua give a large contribution in the region in which 
the pressure as a function of the temperature raises quickly. This result is a possible 
indication that topological non-trivial vacua are important for the thermodynamics of 
QCD in the neighborhood of the phase transition. 

In Chapter 6 the thermodynamics of the linear and nonlinear sigma model in four 
dimensions is investigated. These models are low-energy effective theories of QCD. 
The pressure of these models is calculated to NLO in 1 /N. The most important re- 
sult of these calculations is a prediction for the pressure of two-flavor QCD at low 
temperatures. The NLO correction is compared to predictions of this correction in the 
literature. These predictions turn out to be wrong. Finally a relatively low upper bound 
on the mass of the sigma meson is found. The mass of the sigma meson has not yet 
been determined very accurately experimentally. If this upper bound is a result of the 
approximations or occurs really in QCD should be investigated further. 

The phase diagram of the NIL model is examined in Chapter 7. The NJL model 
is a low-energy effective theory of QCD as well. The phase diagrams are calculated 
by minimizing the effective potential. The calculations have been carried out for sit- 
uations in which the densities of the up, down and the strange quark are different. 
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such as in a compact star or in a heavy-ion collision. Next to the phase in which one 
can find ordinary matter, different color-superconducting phases and phases in which 
the pseudoscalar mesons condense are found. The results of this chapter are an ex- 
tension of earlier calculations in which either only the color-superconducting phases 
or the phases in which the pseudoscalar mesons condense were taken into account. 
The main result is that these previous calculations are extended, such that one gets a 
more complete picture of the phase diagram. Furthermore the competition between 
the color-superconducting phase and the phase in which pseudoscalar mesons con- 
dense is studied. It is found that they are separated by a first-order phase transition. A 
possible extension of this work would be the inclusion of color and electric neutrality 
constraints, such that obtain the phase diagram of a compact star. 



Samenvatting 



Het standaardmodel is de natuurkundige theorie die de deeltjes waaruit materie is 
opgebouwd en hun onderlinge wisselwerkingen (met uitzondering van de zwaarte- 
kracht) beschrijft. Dit model heeft tot nu toe, op het aantonen dat neutrino 's een hele 
kleine massa hebben na, alle experimented testen glansrijk doorstaan. Hoewel de 
wisselwerkingen tussen de deeltjes onderling dus goed begrepen zijn, is het gedrag 
in situaties waarbij een heleboel deeltjes met elkaar wisselwerken vaak beperkt be- 
kend. Een voorbeeld van zo'n situatie is materie onder extreem hoge temperaturen en 
dichtheden, zoals tijdens de oerknal en in hele compacte sterren. In dit proefschrift 
zijn enkele berekeningen uitgevoerd die tot een beter begrip van materie onder zulke 
extreme omstandigheden kunnen leiden. 

Van atomen naar quarks 

De materie om ons been is opgebouwd uit atomen. Een atoom bestaat uit een posi- 
tief geladen kern waar elektronen, die negatief geladen zijn, zich omheen bevinden. 
De elektronen blijven in de buurt van de kern omdat positief en negatief geladen 
deeltjes elkaar via de elektromagnetische wisselwerking aantrekken. De elektromag- 
netische kracht tussen twee ladingen wordt veroorzaakt door de uitwisseling van/o- 
tonen (hetzelfde soort deeltjes waaruit licht bestaat). Deze kracht wordt beschreven 
door quantumelektrodynamica (QED), de quantumtheorie van het elektromagnetisme. 
QED is het best geteste onderdeel van het standaardmodel. 

De kern bevat bijna alle massa van het atoom en is ook weer opgebouwd uit deel- 
tjes: positief geladen protonen en neutrale neutronen. De kern van het waterstof atoom 
bestaat uit een proton, terwijl bij voorbeeld het stabiele goudatoom 79 protonen en 118 
neutronen bevat. 

Omdat protonen elektrisch geladen zijn, zorgt de elektromagnetische wisselwer- 
king voor een afstotende kracht tussen de protonen. Desondanks zijn er veel stabiele 
atoomkernen waarin de protonen en neutronen toch bij elkaar blijven. Daarom moet er 
nog een aantrekkende kracht tussen de protonen en de neutronen bestaan. Deze kracht 
wordt de sterke wisselwerking genoemd. Het blijkt dat de aantrekkende kracht tussen 
protonen en neutronen een gevolg is van de wisselwerkingen tussen de quarks, dat zijn 
de deeltjes waaruit protonen en neutronen zijn opgebouwd. 
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Een proton is opgebouwd uit drie quarks, twee zogenaamde up quarks en een down 
quark. Een neutron bestaat uit een up quark en twee down quarks. Naast de elektrische 
lading hebben de quarks nog een zogenaamde kleurlading. Quarks oefenen een kracht 
op elkaar uit door uitwisseling van gluonen. Deze kracht is op hele kleine afstand 
veel sterker dan de elektromagnetische kracht, waardoor een proton en ook sommige 
atoomkernen stabiel zijn. 

Naast het up en het down quark zijn er nog vier andere zwaardere soorten quarks 
bekend, te weten het strange, charm, bottom en top quark. Met al deze quarks kunnen 
allerlei (instabiele) deeltjes gevormd worden. Maar vanwege de onderlinge interacties 
kunnen lang niet alle combinaties van quarks voorkomen. Bijvoorbeeld, een quark als 
vrij deeltje is onmogelijk. Onder normale omstandigheden vormen quarks altijd met 
z'n drieen een deeltje (bijvoorbeeld een proton of een neutron) of met z'n tweeen. In 
het laatste geval vormt het quark tezamen met een antiquark een deeltje, een zoge- 
naamd meson. Een voorbeeld hiervan is het positief geladen pion, dat een gebonden 
toestand van een up quark en een anti-down quark is. 

De wisselwerkingen tussen de quarks worden beschreven door quantumchromo- 
dynamica (QCD) (zie formule 1.1). QCD is net als QED een onderdeel van het stan- 
daardmodel. Met behulp van deze theorie kunnen bepaalde uitkomsten van botsings- 
experimenten tussen bijvoorbeeld twee protonen worden voorspeld. De resultaten van 
deze experimenten stemmen overeen met de voorspellingen, wat dus betekent dat de 
fundamentele wisselwerkingen tussen de quarks in principe goed begrepen zijn. 

Materie onder extreme omstandigheden 

In dit proefschrift is theoretisch onderzocht wat met materie opgebouwd uit quarks 
gebeurt onder extreem hoge temperaturen, zoals tijdens de oerknal en bij hele grote 
dichtheden, zoals in een compacte ster. In beide gevallen hebben we natuurlijk te 
maken met een heleboel quarks, veel meer dan in een botsingsexperiment tussen twee 
protonen. Net zoals het ondoenlijk is om te voorspellen hoe een individueel water- 
molekuul zich door water beweegt, kunnen we in dat geval geen uitspraken doen over 
het gedrag van ieder individueel quark. Maar het is wel mogelijk om uitspraken te 
doen over het gedrag van alle quarks tezamen, net zoals water te karakteriseren is 
door collectieve grootheden, als de temperatuur, de dichtheid en de druk. De relaties 
tussen deze grootheden kunnen met behulp van thermodynamica (warmteleer) worden 
uitgerekend en zijn ondere andere van belang voor het begrijpen van het gedrag van 
compacte sterren en het ontstaan van het heelal. 

Water kan in verschillende toestanden (ook wel fasen genoemd) voorkomen, te 
weten ijs, vloeibaar en damp. Al deze fasen zijn natuurlijk opgebouwd uit dezelfde 
watermolekulen, alleen gedraagt de verzameling watermolekulen zich in iedere fase 
complect anders. In ijs zijn de watermolekulen relatief sterk aan elkaar gebonden. 
Als de temperatuur verhoogd wordt, gaan de watermolekulen sneller trillen, op een 
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gegeven moment trillen ze zo hard dat ze niet meer in een soort roostervorm bij elkaar 
kunnen worden gehouden, water gaat dan over naar de vloeibare fase. Bij de overgang 
naar waterdamp gebeurt ongeveer hetzelfde. 

Algemeen wordt aangenomen dat iets dergelijks ook met quarks moet gebeuren. 
Zoals eerder is aangegeven, vormen quarks met z'n tweeen of met z'n drieen een 
deeltje zoals het neutron. Neem nu eens in gedachten een groot aantal neutronen in 
een doos. Als de temperatuur laag is dan blijven de quarks bij elkaar, maar als de 
temperatuur verhoogd wordt kan het gebeuren dat de quarks zo hard gaan trillen, dat 
ze niet meer bij elkaar kunnen worden gehouden in een neutron. De materie komt dan 
in een nieuwe fase die het quark- gluonplasma wordt genoemd. In deze fase hoeven 
de quarks niet langer met z'n tweeen of drieen bij elkaar te blijven. De temperatuur 
waarbij de faseovergang optreedt is ongeveer twee biljoen graden Celsius, dat is maar 
liefst honderdduizend keer de temperatuur in het binnenste van de zon. Vlak na de 
oerknal was het heelal waarschijnlijk eventjes een quark-gluonplasma. Er zijn sterke 
aanwijzingen dat het quark-gluonplasma onlangs voor een hele korte tijd in het la- 
boratorium is gemaakt tijdens een experiment in het Brookhaven National Laboratory, 
New York, waar goudkernen met enorme snelheden tegen elkaar werden geschoten. 

De neutronen in de doos kunnen natuurlijk ook op elkaar gedrukt worden. Op een 
gegeven moment wordt de dichtheid zo groot dat alle neutronen tegen elkaar aanzit- 
ten. Als er dan nog harder gedrukt wordt, dan gaan de neutronen overlappen, zodat ze 
hun identiteit als individuele deeltjes verliezen. In dit geval kan er kan een faseover- 
gang optreden naar een zogenaamde kleur-supergeleidende fase. Dit gebeurt bij enorm 
grote dichtheden van ongeveer honderd biljoen kilogram per liter. In theorie zou die 
faseovergang bereikt kunnen worden in het binnenste van een neutronenster. Een neu- 
tronenster is een uitgebrande ster met een enorm hoge dichtheid. Tot nu zijn er echter 
nog geen neutronensterren gevonden die zo'n hoge dichtheid hebben dat ze zich met 
zekerheid in de kleur-supergeleidende fase bevinden. 

De bovenstaande beweringen over het gedrag van deeltjes onder extreme om- 
standigheden kunnen in principe wat preciezer gedaan worden door berekeningen met 
QCD uit te voeren. Het blijkt dat de krachten tussen de quarks minder sterk worden 
naarmate de energie van de quarks groter wordt. Dit fenomeen heet asymptotische vrij- 
heid. Als gevolg daarvan is het mogelijk analytische berekeningen (berekeningen met 
pen en papier) met QCD bij hele hoge energieen uit te voeren. Maar deze energieen 
zijn nog veel hoger dan de energieen die de quarks hebben rond de faseovergangen. 
Een manier om QCD bij lagere energieen te bekijken is met behulp van computersimu- 
laties {lattice QCD). Het probleem van deze computersimulaties is dat ze niet goed 
werken bij grote dichtheden. Om toch iets te voorspellen en te begrijpen over het 
gedrag van QCD bij temperaturen en dichtheden rond de faseovergangen kun je een- 
voudiger theorieen bekijken die QCD tot op zekere hoogte benaderen (zoals in hoofd- 
stuk 6 en 7 van dit proefschrift) of die bepaalde aspecten met QCD gemeen hebben 
(zoals in hoofdstuk 4 en 5) (vandaar de uitdrukking QCD-inspired theories in de titel). 
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Een belangrijke uit te rekenen grootheid is de druk (aangeduid met het symbool 
P) van materie die uit quarks is opgebouwd, als functie van de temperatuur (T). Het 
resultaat van de exacte berekeningen is afgebeeld in figuur 1.1. De grijze banden 
geven de foutmarges aan. Duidelijk is te zien dat voor lagere temperaturen (in het 
bijzonder rond de faseovergang) de resultaten van de berekeningen onbetrouwbaar 
zijn (omdat de krachten tussen de quarks dan te groot worden). Het resultaat van de 
computerberekeningen is voor een verschillend aantal quarks afgebeeld in figuur 1.2. 
De resultaten van deze computerberekeningen zijn betrouwbaar voor temperaturen in 
de buurt van de faseovergang (ongeveer 170 MeV, dat correspondeert met 2 biljoen 
graden Celsius) en hoger. Om de druk bij lagere temperaturen te voorspellen kunnen 
QCD-achtige theorieen gebruikt worden, zoals in hoofdstuk 6 van dit proefschrift is 
gedaan. In figuur 6.6 is het resultaat van de berekening van de druk afgebeeld. De 
onderbroken lijn in deze figuur is een voorspelling is van de druk van QCD bij lage 
temperatuur. 

Een ander interessant voorbeeld is het fasediagram van QCD. Het fasediagram 
geeft, afhankelijk van bij voorbeeld de dichtheid en de temperatuur, aan in welke fase 
materie zich bevindt. In figuur 1.3 is het fasediagram van QCD geschetst. Op de hori- 
zontale as van dit diagram staat hb afgebeeld, een grootheid die een directe relatie heeft 
met de dichtheid (hoe groter ^b, hoe groter de dichtheid). Op de verticale as staat de 
temperatuur. Het quark-gluonplasma en de kleur-supergeleidende fase zijn in dit dia- 
gram terug te vinden. Het is belangrijk te weten dat dit diagram een schets is, het is 
namelijk niet precies bekend waar de lijnen die de verschillende fasen van elkaar schei- 
den in het diagram liggen. Dat komt doordat dit diagram gebaseerd is op berekeningen 
aan theorieen die QCD benaderen in een bepaald gebied. In het fasediagram afgebeeld 
in figuur 1.3 zijn de dichtheden van de verschillende soorten quarks allemaal gelijk. In 
hoofdstuk 7 van dit proefschrift zijn, met behulp van een QCD-achtige theorie, fasedia- 
grammen uitgerekend waarbij de dichtheden van de quarks allemaal verschillend zijn 
gekozen. 

Overzicht en resultaten van dit proefschrift 

Dit proefschrift getiteld, "Thermodynamics of QCD-inspired theories" (Thermody- 
namica van QCD-achtige theorieen), gaat dus over het gedrag van materie onder ex- 
treme omstandigheden. 

In hoofdstuk 1 wordt een algemene inleiding gegeven. Eerst wordt ingegaan op 
de vraag wat er zou kunnen gebeuren met materie onder extreme omstandigheden. 
Daama worden in het kort de relevante aspecten van quantumchromodynamica (QCD) 
uitgelegd. Vervolgens worden eerdere berekeningen van de druk en het fasediagram 
en de problemen die daar bij komen kijken behandeld. Ook wordt in hoofdstuk 1 aan- 
dacht besteed aan situaties waarbij zulke extreme omstandigheden kunnen voorkomen, 
te weten in het heelal vlak na de oerknal, tijdens zware-ionenbotsingen en in neutro- 
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nensterren. Tot slot wordt uitgelegd wat er bedoeld wordt met QCD-achtige theorieen. 

Hoofdstuk 2 gaat over hoe met een zogenaamde quantumveldentheorie als QCD 
berekeningen kunnen worden uitgevoerd bij eindige temperatuur en dichtheid. Daar- 
voor wordt een statistische analyse gebruikt. Eerst wordt er ingegaan op hoe dat in zijn 
werk gaat voor een klassieke theorie, daarna worden de resultaten veralgemeniseerd 
voor een quantummechanische aanpak. Met behulp van het zogenaamde pad-integraal 
formalisme wordt uitgelegd hoe dit kan worden gebruikt voor een quantumveldentheo- 
rie. In berekeningen bij eindige temperatuur moeten vaak sommaties over een oneindig 
aantal termen worden uitgerekend. Een methode om dat exact te doen wordt uitgelegd, 
en er wordt ingegaan op hoe er met oneindigheden, die in deze berekeningen kunnen 
optreden, moet worden omgegaan. Ook wordt ruim aandacht besteed aan hoe deze 
sommen voor ingewikkelder gevallen met de computer kunnen worden berekend. 

Het derde hoofdstuk behandelt de zogenaamde effectieve potentiaal en de IjN 
benadering. De effectieve potentiaal is een belangrijke grootheid, aangezien door 
het minimaliseren van die potentiaal de druk en de fasediagrammen kunnen worden 
verkregen. De 1 IN benadering is een methode waarbij inzicht kan worden verkregen 
in het lage energie gedrag van de theorieen die in dit proefschrift worden behandeld. 
Er wordt uitgelegd hoe deze benadering kan worden afgeleid door het introduceren van 
extra velden. Tot slot wordt ingegaan op het feit dat in berekeningen van de effectieve 
potentiaal er oneindigheden kunnen ontstaan die van de temperatuur afhangen. Zulke 
oneindigheden vormen een grote bedreiging voor de betekenis en bruikbaarheid van 
de antwoorden. Een belangrijk resultaat uit dit proefschrift is dat deze oneindigheden 
temperatuursonafhankelijk worden op het minimum van de effectieve potentiaal. Dat 
is juist het punt waar de druk en het fasediagram worden bepaald. In het minimum 
zijn deze oneindigheden geen probleem omdat ze dan door subtracties en herdefinities 
van bepaalde constanten op een systematische manier kunnen worden verwijderd zon- 
der de betekenis van de antwoorden te beinvloeden. Deze procedure wordt ook wel 
renormalisatie genoemd. 

Het onderwerp van hoofdstuk 4 is de thermodynamica van het niet-lineaire sigma 
model in twee dimensies (een ruimte- en een tijddimensie). Dit model heeft bepaalde 
aspecten met QCD gemeen, zoals asymptotische vrijheid. Maar met behulp van de 
\IN benadering kan hier het lage energie gedrag wel analytisch worden berekend. 
Door dit model te bestuderen kunnen we lets leren van het gedrag van QCD bij eindige 
temperatuur, bijvoorbeeld hoe een asymptotische vrije theorie zich gedraagt als er van 
lage naar hoge energie wordt gegaan. In dit hoofdstuk is de effectieve potentiaal uit- 
gerekend met behulp van de 1 IN benadering tot en met de eerste niet-triviale correctie 
(ook wel next-to-leading order (NLO) genoemd). Er worden temperatuursafhanke- 
lijke oneindigheden gevonden, die blijken te kunnen worden gerenormaliseerd op het 
minimum. Vervolgens wordt de druk als functie van de temperatuur uitgerekend. Met 
dit resultaat kunnen we ondere andere concluderen dat de 1/A^ benadering een goede 
rekenmethode is. Verder lijkt het gedrag van de druk gedeeld door de druk in the limiet 
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van oneindige temperatuur, onafhankelijk van A'^. Een soortgelijk resultaat is gevonden 
in computerberekeningen aan QCD. 

Hoofdstuk 5 gaat over de thermodynamica van het CP^'^ model in twee dimen- 
sies. Dit model is een uitbreiding van het niet-lineaire sigma model waarbij ook elek- 
tromagnetisme wordt meegenomen. Ook dit model heeft weer een aantal aspecten 
met QCD gemeen, het bevat bijvoorbeeld net als QCD zogenaamde topologisch niet- 
triviale vacua. Het effect van deze vacua op de druk is onderzocht. Daarbij wordt 
gebruik gemaakt van dezelfde methoden als in hoofdstuk 4. Uit de berekeningen volgt 
dat de topologisch niet-triviale vacua voor temperaturen waarbij de druk relatief snel 
stijgt een grote bijdrage leveren. Dit resultaat is een mogelijke aanwijzing dat deze 
niet-triviale vacua rond de faseovergangen in QCD belangrijk zouden kunnen zijn. 

In hoofdstuk 6 wordt de thermodynamica van het lineaire en het niet-lineaire sigma 
model in een tijd- en drie ruimtedimensies onderzocht. Algemene wordt aangenomen 
dat deze modellen een goede benadering van QCD bij lage energie zijn. Opnieuw 
is met behulp van de 1 IN benadering de druk uitgerekend tot en met de eerste niet- 
triviale correctie. Het belangrijkste resultaat van deze berekening is een voorspelling 
van de druk van QCD met twee soorten quarks bij lage temperaturen. Dit resultaat is 
vergeleken met eerdere schattingen uit de literatuur, waarvan geconcludeerd wordt dat 
ze foutief zijn. Tot slot wordt er een vrij lage bovengrens gevonden op de massa van 
het zogenaamde sigma meson. Deze massa is experimented nog niet goed vastgesteld. 
Of deze bovengrens een gevolg is van de benaderingsmethode of ook echt in QCD 
voorkomt zou verder moeten worden onderzocht. 

Hoofdstuk 7 gaat tenslotte over de fasediagrammen van het NJL model. Dit model 
geldt ook als een theorie die QCD kan beschrijven bij lage energieen. Door het mini- 
maliseren van de effectieve potentiaal worden de fasediagrammen uitgerekend. De 
berekeningen zijn uitgevoerd voor situaties waarbij de dichtheden van de up, down 
en strange quark verschillend kunnen zijn, zoals bijvoorbeeld in een ster of in een 
zware-ionenbotsing. Naast de fase waarin materie zich bij normale omstandighe- 
den bevindt, zijn verschillende kleur-supergeleidende fasen en fasen waarin pseu- 
doscalaire mesonen condenseren gevonden. De resultaten uit dit hoofdstuk zijn een 
uitbreiding op eerdere berekeningen uit de literatuur waarin of alleen rekening met de 
kleur-supergeleidende fase of met de fase waarin pseudoscalaire mesonen condenseren 
wordt gehouden. In dit hoofdstuk worden beide mogelijkheden meegenomen. Het 
belangrijkste resultaat is dat de eerdere berekeningen van fasediagrammen zijn uitge- 
breid. Verder is de competitie tussen de twee typen fasen onderzocht, een belangrijk 
resultaat is dat ze gescheiden zijn door een zogenaamde eerste-orde faseovergang. 

Dit onderzoek heeft er toe geleid dat het begrip van het gedrag van materie onder ex- 
treme omstandigheden lets is toegenomen. Toch blijven er nog een heleboel zaken 
onbegrepen, zoals bijvoorbeeld de precieze vorm van het QCD fasediagram. Hopelijk 
zullen verbeterde (computer )berekeningen en nieuwe experimenten (bijvoorbeeld bij 
CERN in Geneve) ons in de toekomst meer leren. 
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